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Abstract 



The wave function of a composite system is defined in relativity on a space-time 
surface. In the exphcitly covariant hght-front dynamics, reviewed in the present article, 
the wave functions are defined on the plane uj-x = 0, where u is an arbitrary four- vector 
with cu^ — 0. The standard non-covariant approach is recovered as a particular case for 
CO — (1,0,0, —1). Using the light-front plane is of crucial importance, while the explicit 
covariance gives strong advantages emphasized through all the review. 

The properties of the relativistic few-body wave functions are discussed in detail and 
are illustrated by examples in a solvable model. The three-dimensional graph technique for 
the calculation of amplitudes in the covariant light-front perturbation theory is presented. 

The structure of the electromagnetic amplitudes is studied. We investigate the 
ambiguities which arise in any approximate light-front calculations, and which lead to a 
non-physical dependence of the electromagnetic amplitude on the orientation of the light- 
front plane. The elastic and transition form factors free from these ambiguities are found 
for spin 0, 1/2 and 1 systems. 

The formahsm is applied to the calculation of the relativistic wave functions of two- 
nucleon systems (deuteron, scattering state), with particular attention to the role of their 
new components in the deuteron elastic and electrodisintegration form factors and to their 
connection with meson exchange currents. Straigthforward applications to the pion and 
nucleon form factors and the p — tt transition are also made. 
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Chapter 1 
Introduction 



1.1 The relevance of relativity in few-body systems 

The relevance of a coherent relativistic description of few-body systems, both for bound 
and scattering states, is now well recognized in nuclear as well as in particle physics. 

This is already clear in particle physics for the understanding of the wave functions 
of the valence quarks in the nucleon or in the pion, as revealed for instance in exclusive 
reactions at very high momentum transfer. The typical example of such a reaction is 
elastic scattering, and the extraction of the electromagnetic form factors, with their well 
known behavior at very high momentum transfer, the so-called scaling laws expected from 
QCD. 

The need for a coherent relativistic description of few-body systems has also be- 
come clear in nuclear physics in order for instance to check the validity of the standard 
description of the microscopic structure of nuclei in terms of mesons exchanged between 
nucleons. In this case, electromagnetic interactions play also a central role in "seeing" 
meson exchanges in nuclei. The forthcoming experiments at Thomas Jefferson National 
Accelerator Facility (former CEBAF) at momentum transfer of a few (GeV/c)^ are here 
of particular importance. 

In both domains, it is obligatory to have, first, a relativistic description of the 
bound and scattering states. It is also necessary to have a consistent description of the 
electromagnetic current operator needed to probe the system. This is mandatory in order 
to have meaningful predictions for the various cross-sections. 

A few relativistic approaches have been developed in the past ten years in order to 
meet these goals. Among them, two have received particular attention in the last few 
years. 

The first one is based on the Bethe-Salpeter formalism [Q] or its three-dimensional 
reductions. The Bethe-Salpeter formalism is four- dimensional and explicitly Lorentz co- 
variant. The calculational technique to evaluate electromagnetic amplitudes is based on 
Feynman diagrams and associated rules. Three-dimensional reductions result in equations 
of the quasi-potential type. 

The second one is Light-Front Dynamics (LFD) In this case, the state vector 
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describing the system is expanded in Fock components with increasing number of particles. 
The state vector is defined on a surface in four-dimensional space-time which should 
be indicated explicitly. The Fock components - the relativistic wave functions in this 
formalism - are the direct generalization of the non-relativistic wave functions. 

In the non-relativistic limit, (c — >■ oo), the wave function is defined at time t = 0, 
and the time evolution is governed by the Schrodinger equation, once the Hamiltonian of 
the system is known. This is the "instant" form of dynamics. Physical processes are thus 
calculated according to old fashioned (time ordered) perturbation theory. This form of 
dynamics is however not very well suited for relativistic systems, since the interaction of 
a probe (the electron, for instance) with the constituents of the system is not separated 
from its interaction with the vacuum fiuctuations. Moreover, the plane t = cte is not 
conserved by a Lorentz boost. 

In the standard version of LFD, the wave function is defined on the plane ^ + § = cte 
0. It is also equivalent to the usual equal time formalism in the infinite momentum frame. 
From a qualitative point of view, all the physical processes become as slow as possible 
because of time dilation in this system of reference. This greatly simplifies the description 
of the system. The investigation of the wave function is equivalent to make a snapshot 
of a system not spoiled by vacuum fiuctuations. It is thus very natural in the description 
of high energy experiments like deep inelastic scattering. The calculational technique is 
here based on the Weinberg rules [|]. 

This formulation has however a serious drawback since the equation of the plane 
t + z = cte (we take here and in the following c = 1) is still not invariant by an arbitrary 
rotation in space and time (Lorentz boost). This plane breaks rotational invariance. As 
we shall see later on, this fact has many important consequences as far as the construction 
of bound states (or scattering states) of definite angular momentum is concerned, or in 
the calculation of electromagnetic amplitudes. 

1.2 Why light- front dynamics? 

We shall present in this review a covariant formulation of LFD which provides a simple, 
practical and very powerfull tool to describe relativistic few-body systems, their bound and 
scattering states, as well as their physical electromagnetic amplitudes. In this formulation, 
the state vector is defined on the plane characterized by the invariant equation u-x = 0, 
where u is an arbitrary light-like four vector uj = (ci;o,ci;), with tu^ = 0, ^. With the 
particular choice u = (1,0,0,-1), we recover the standard LFD defined on the plane 
t + z = 0. The covariance of our approach is realized by the invariance of the light- 
front plane u-x = under any Lorentz transformation of both u and x. This implies in 
particular that u cannot be kept the same in any system of reference, as it is the case in 
the usual formulation of LFD with u = (1, 0, 0, —1). 

There is of course equivalence, in principle, between these approaches. Any exact 
calculation of a given electromagnetic process should give the same physical cross-section, 
regardless of the relativistic formalism which is used. In practical calculations however. 
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one can hardly hope to carry out an exact calculation of physical observables starting 
from "first principles" , although some simple theoretical systems can of course be used as 
test cases (the Wick-Cutkosky model for instance). 

From a practical point of view, one may thus be led to choose a particular formalism 
depending on the system one is interested in. This choice is of course dependent on various 
criterions, as well as function of personal taste. Among these criterions, let us mention 
those we think are the most important: 

i) The formalism should enable us to have physical insights into the various processes 
under consideration, at each step of the calculation. 

a) It should have a direct and transparent non-relativistic limit in order to gain 
more physical intuition from our present knowledge in the non-relativistic domain. This 
is particularly important in nuclear physics where a lot is already known from the non- 
relativistic description of the microscopic structure of nuclei and their electromagnetic 
interactions at the scale of 1 GeV or less 0. 

in) It should provide a as simple as possible calculational procedure to evaluate 
physical processes and compare them with experimental results. 

We would like to show in this review that the covariant formulation of LFD is, 
according to these criterions, of particular interest. As we shall see extensively in the 
following, it has definite advantages as compared to both the Bethe-Salpeter formalism, 
and the usual formulation of LFD. Let us recall the most important ones below. 

z) The calculational formalism - time ordered graph technique (in the light-front 
time) in contrast to old fashioned perturbation theory - does not involve vacuum fluctua- 
tion contributions. This strongly simplifies a priori the physical picture and calculations. 

a) In this approach, the wave functions - the Fock components of the state vector 
- satisfy a three-dimensional equation, and have the same physical meaning (probability 
amplitudes) as the non-relativistic wave functions. Their non-relativistic limit is thus 
explicit from the general structure of the wave function. This enables a transparent link 
with non-relativistic approaches in first 1/m order and in particular with the contribution 
of the dominant meson exchange current (the so-called pair term) in the two-nucleon 
systems. 

in) The explicit covariance of LFD is very important in practical applications: it 
allows to construct states with definite total angular momentum to separate contributions 
of relativistic origin from non-relativistic ones, and simplifies very much the calculations 
in the framework of a special graph technique. In these respects, the advantages of the 
covariant formulation in comparison to ordinary LFD are the same as the advantages of 
the Feynman graph technique in comparison to old fashioned perturbation theory. 

iv) A very important property of relativistic wave functions and off-shell amplitudes 
is their dependence on the orientation of the light-front plane. It takes place both in 
non-covariant and covariant approaches. In the covariant approach this dependence is 
parametrized explicitly in terms of the four-vector uj. On the contrary, exact on-shell 
physical amplitudes should not depend on the orientation of light-front plane. However, 
in practice, this dependence survives due to approximations. The covariant representation 
of the electromagnetic amplitudes allows to separate the physical form factors from the 
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unphysical contributions. 

As already mentioned, the present approach differs from the standard ones by the 
parametrization of the current choice of the z direction by an arbitrary one n. A second 
major difference is the reference to a particular field theory inspired dynamics, which 
implies that our approach describes states with an unfixed number of particles. This is 
a source of many problems, especially related to the mass dependence of the interaction. 
These problems, which will receive a particular attention, have a strong relationship with 
those encountered in nuclear physics when using an energy dependent nucleon-nucleon 
interaction (Bonn-E potential 0), or with the contribution of recoil and norm corrections 
to meson-exchange currents [||. 

Our aim in this review is to present the various facets of the covariant formulation 
of LFD, and compared them to the standard formulation. For completeness, we shall 
also make contacts with the BS formalism. The relevance of Light-Front Quantization 
in quantum field theory is the subject of an intense present activity. Among others, the 
non-perturbative problems of how the condensates present in many theories appear on 
the light-front, or how renormalization should be implemented, are not yet completely 
solved ip. However, directions of research are clearly identified (zero modes in LFD in 
particular). These problems are outside the scope of this review, although the explicit 
covariance of our formalism may be of particular interest in solving these problems. 

We believe that the real advantages of the covariant formulation of LFD can be 
apprehended in practice, i.e., in applications to relativistic nuclear and particle physics 
and field theoretical problems. The first applications to relativistic nuclear and particle 
physics are reviewed in the present paper, with particular attention to two-nucleon systems 
for which many experimental data exist at present or are expected in the near future. 

Many extensive or review papers have been devoted to the description of few- 
body systems in relativistic approaches. Some of them are based on different versions 
of the quasi-potential equations 0, |l^ . Many applications were made using the three- 
dimensional Gross equation O, O, O, |15|, |I6[. Two- and three- nucleon systems were 



investigated in the framework of the Bethe-Salpeter equation in refs. [|17| , [18], [T^, |20| (see 
for a review |21|). The solution of the Bethe-Salpeter equation was also found and applied 
to deep inelastic scattering on deuteron and to the deuteron electrodisintegration in refs. 



The standard version of the LFD and its applications to few-body systems was 
reviewed in refs. ^ ^ ^ . The first results on the covariant formulation 

of LFD are given in ref. 0. This approach has been investigated and developed also in 



refs. f32 



The content of this review is the following. We develop in chapter |^ the general 
properties of the covariant formulation of LFD, and derive the graph technique associ- 
ated to it. Chapter ^ is devoted to the properties of the two-body wave function: spin 
structure, equation for the wave function, as well as its connection to the Bethe-Salpeter 
amplitude. We pay particular attention to the construction of the angular momentum 
operator. We apply our formalism to the two-nucleon system in chapter ^ by first de- 
riving the nucleon-nucleon potential in this formalism. The two-nucleon wave function 



7 



(deuteron and ^5*0 wave function) is then constructed in chapter |^. The general structure 
of the electromagnetic amplitude is detailed in chapter where we show how to extract 
the physical form factors for spin 0, 1/2 and 1 systems, as well as transition form factors. 
We discuss in chapter ^ the electromagnetic form factor of the simplest states with J = 
and 1 of the Wick-Cutkosky model, as well as a few hadronic systems (pion, nucleon, 
...). We apply our formalism to the two- nucleon system (deuteron form factors and elec- 
trodisintegration cross-section). We also discuss, in a 1/m expansion, the relationship 
between our formalism and first order relativistic corrections taken into account as meson 
exchange currents (the so-called pair term) in non-relativistic approaches. 

The chapters 2, 3 and 6 are quite general and apply to any system. Chapters 4, 
5 and 7 are direct applications to the two- nucleon systems, as well as three quark and 
quark- ant iquark systems. 
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Chapter 2 

Covariant formulation of light-front 
dynamics 

We detail in this chapter the general properties of the covariant formulation of LFD. 
In contrast to the standard approach, the transformations of the coordinate system and 
of the light-front plane can be done independently from each other. These two types of 
transformations entail the corresponding transformation properties of the state vector and 
its Fock components. Particular attention is paid to the angular momentum operator. 
This formal field-theoretical introduction (sect. |2.1| ) can be omitted by a reader not 
interested in these details. 

We then present the graph technique associated to this formulation, for particles 
of spin 0, 1/2 and 1. We illustrate peculiarities of this graph technique by a few simple 
examples. 

2.1 Transformation properties of the state vector 

The state vector is defined in general on a hypersurface in space-time and therefore de- 
pends dynamically on its position. For example, the non-relativistic wave function iplt) 
depends dynamically on translations of the plane t = cte, i.e. on time t (by the trivial 
phase factor exp{—iEt) for a bound state, where the binding energy E is determined by 
dynamics). In ordinary LFD, the wave function is defined on the plane t + z = = cte. 
It depends dynamically on any translation of this plane. However, some Lorentz transfor- 
mations and rotations change the orientation of this plane |0. These transformations are 
thus also dynamical ones. That means lack of explicit covariance (i.e., the impossibility to 
transform the state vector from one reference system to another one without knowledge of 
the dynamics). This does not mean the absence of Lorentz covariance at all, since there 
exist anyhow a closed system of generators of the Poincare group, and the observable 
amplitudes calculated exactly would be covariant. In practical approximate calculations, 
however, the covariance is lost. 

In the covariant formulation of LFD, the wave function is defined on the general 
plane u-x = a, where oo is an arbitrary four vector restricted to the condition cj^ = 0,and 
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a is the "light-front time". When there is no need to refer to the a evolution, we shall 
take (7 = 0. The kinematical transformations of the system of reference are thus separated 
from the dynamical transformations of the plane uj-x — a. The dynamical dependence of 
the wave function on the light-front plane results in that case in their dependence on uj. 

This separation of kinematical and dynamical transformations has a definite ad- 
vantage in the sense that it provides a definite prescription for constructing bound and 
scattering states of definite angular momentum. Since the total angular momentum of a 
composite system is determined by the transformation properties of its wave function un- 
der rotation of the coordinate system, the construction of systems with definite J is now 
purely kinematical, in contrast to the formulation on the plane t + z = . The dynamical 
part of this problem, resulting from the presence of the interaction in the generators of 
the Poincare group which change the position of the light-front plane, is separated out 
and replaced by the so-called angular condition. 



2.1.1 Kinematical transformations 

Let us first specify the transformation properties of the state vector with respect to 
transformations of the coordinate system. We will use for this purpose a field-theoretical 
language. The operators associated to the four-momentum and four-dimensional angular 
momentum are expressed in terms of integrals of the energy-momentum Ty^y and the 
angular momentum M^^ tensors over the light-front plane lo-x = a, according to: 

P,^l T^uU^'^K^-x - a)d'x = P; + , (2.1) 

J,. = / M>,5(a;-x - a)d^x - + , (2.2) 

where the and int superscripts indicate the free and interacting parts of the operators 
respectively. For generality, we consider here the light-front time cr 7^ 0. 

The description of the evolution along the light-front time a implies a fixed value 
of the length of cJ, or, equivalent ly, of ujq. This is necessary in order to have a scale of 
(J. However, the most important properties of the physical amplitudes following from 
covariance do not require to fix the scale of uj and will be invariant relative to its change. 

We work in the interaction representation in which the operators are expressed in 
terms of free fields. For example, for a scalar field ^{x), the free operators P° have the 
form: 

^ J a\k)a{k)ky d^k , (2.3) 
Ji. = j<^\MkY{k,^^-k.^^d^k, (2.4) 

where a* and a are the usual creation and destruction operators, with [a{k) , a* {k')] — 
5^{k — k'). The operators P*"* and J'"* contain the interaction Hamiltonian i7'"*(x): 



P^"* = / W\x)6{uj-x - a) d^x , (2.5) 
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J. 



int 



H^^ {x){x^uj^ — Xi,Ufj)5{uj-x — a) d X 



(2.6) 



For the particular applications considered in this review, we do not need to develop the 
field theory in its full form. We therefore do not pay attention here to the fact that 
the field-theoretical Hamiltonian if*"* (a;) is usually singular and requires a regularization. 
The regularization of amplitudes in our formulation will be illustrated by the example of 
a typical self-energy contribution at the end of this chapter. 

Eq.( |2.5| ) is consistent with the expectation that only the component along the 
direction uj of the light- front "time" has a dynamical character in the light-front formalism 
P]. Under translation x ^ x' = x + a oi the coordinate system A —>■ A', the equation 
uj-x = a takes the form uj-x' = a', where a' = a + uj-a. The state vector is transformed in 
accordance with the law: 



<i^.{a)^<i)l{a') = Upo{a)<i)^{a) , (2.7) 

where Upo{a) contains only the operator of the four- momentum ( p.3|) of the free field: 

f/po(a) = exp(iP°-a) . (2.8) 

The "prime" at (f)'{cr) indicates that (p'icr) is defined in the system A' on the plane u-x' = a 
in contrast to 0(cr) defined in the system A on the plane u-x = a (the value of a being 
the same). The state vector (p'lcr') is defined in A' on the plane u-x' = a', which coincides 
with uj-x = a. Therefore no dynamics enters into the transformation ( |2.7| ). This is rather 
natural, since under translation of the coordinate system the plane u-x = a occupies the 
same position in space while it occupies a new position with respect to the axes of the 
new coordinate system, as indicated in fig. [I|. The formal proof of (^.T]), ( |2.8| ) can be 



found in [33 . 



(ox=a 




Fig.l. Translation of the reference system along the light-front time. 
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the result is similar 33 



In the case of infinitesimal four-dimensional rotations ^ x'^ = gx^ = + e^^x'^ ^ 

ci>u.{(y)^(t>'A(y) = UAg)ct>^{a) , (2.9) 



where u;^ = u;^ + e^^uj'^ and 

1 



2 

The operator is given by 



UAg) = 1 + -Jl,e>^\ (2.10) 



2.1.2 Dynamical transformations 

The properties of the state vector under transformations of the hypersurface are deter- 
mined by the dynamics and follow from the Tomonaga-Schwinger equation 

i6(j)/6a{x) = H"'\x) . (2.11) 



/^From the definition of the variational derivative in (|2.11|) we obtain: 

i5(P = W\x) 5V{x) , 

where 5V{x) is the volume between the initial surface and the surface obtained from the 
original one by the variation 5a{x) around the point x. 

Under the translation a ^ a + 5a oi the plane, the total increment of the state 
vector is obtained through the increment at each point of the surface: 



j H'''\x)5{u-x - a)d^x 5(7 . (2.12) 



This relation gives the Schrodinger equation. In the interaction representation in the 
light-front time, we have: 

id(l)/d(J = H{a)(j){a) , (2.13) 

where: 

H{a) = J Wf{x)5{uj-x - a)(fx, (2.14) 

and if™*(x) may differ from if*"*(x) because of singularities of the field commutators on 
the light cone. This point is explained below in sect. p.2.1| . 



Similarly, in the case of a rotation of the light-front plane, u ^ ^ uo'^ = u ^ + 5lj^, 
5uj^ = e^^oj^ , we find: 

I ( d d \ 

The increment of the volume over the point x is: 

5V = e^^ x^ J" 5{u-x - o) d^x , (2.16) 
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and it follows from ( p.l2| ) that |P3|| : 



where: 



j;:*0^(a)= VM</)^(a) , (2.17) 



VM = M^M^-^.^) , (2.18) 



and J*^* is given by ( |2.6| ). 

Equation (|2.17|) is called the angular condition. It plays an important role in the 
construction of relativistic bound states, as we shall explain below. 

The transformation of the coordinate system and the simultaneous transformation 
of the light-front plane, which is rigidly related to the coordinate axes, correspond to the 
successive application of the two types of transformations considered above (kinematical 
and dynamical). Thus, under the infinitesimal translation x x' = x+a of the coordinate 
system, A A', and of the plane, we have: 

Mct) ^ (P'M = (1 + iP-a)(l)^{a) . (2.19) 

Note that for the state with definite total four- momentum p (i.e., for an eigenstate of the 
four- momentum operator), the equations (|2.7|) and ( p.l9|) give: 



exp(iP°-a)0(cr) = exp(ip-a)0(cr + cu-a) . (2.20) 

This equation determines the conservation law for the four-momenta of the constituents, 
given in chapter |^. 

2.1.3 Role of the angular condition 

We are interested here in the state vector of a bound system. It corresponds to a definite 
mass M, four-momentum p, total angular momentum J with projection A onto the z axis 
in the rest frame, i.e., the state vector forms a representation of the Poincare group. This 
means that it satisfies the following eigenvalue equations: 

P,^'\p) = P,^'\p), (2.21) 

P2 0JA^p) _ ^2 ^JA(p) ^ (2.22) 

^3 <p'\p) = M X<p'\p) , (2.24) 



where S*^ is the Pauli-Lubanski vector: 



Sfi — ~^fiup'y P ■ (2.25) 



The state vector |p, A)^^ for a given J is normalized as follows: 

^(A',p'b, A)^ = 2po 5^^\p-p') S''^ . (2.26) 
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For convenience we introduce here another notation for the same state vector: 

\pA)^^<Pi\p). 

We omit in (p'^^ip) the argument cr, but show exphcitly the argument p. For simphcity, 
we have left out the u underscript when not absolutely necessary. 

We can now use the angular condition (|2.17|) and replace the operator j^"*, which 

is contained in J^^, in ( p.25| ), by L^^ioj). Introducing the notations: 



M, 



we obtain instead of eqs.( |2.23| ) and ( ^^.241 ): 

W^(j)-^^{p) - 



lJX 



ip) 



-mV(j + 1) (j)-^^{p) 

M X (f)-^^{p) . 



(2.27) 
(2.28) 



(2.29) 
(2.30) 



These equations do not contain the interaction Hamiltonian, once cj) satisfies ^2. 21\ ) and 
( 2.2i\ ). The construction of states with definite angular momentum becomes therefore a 
purely kinematical problem. 

At the same time, the state vector must satisfy the dynamical equation ( |2.17| ). In 
terms of the operators J^y and M^^, the angular condition ( |2.17| ) can be rewritten as: 



The commutation relations between the operators P^, J^j^, M^j/ have the form: 



~[P^J.J J up] 



Qp'yJi 



Vp 1 



1 



P7J 



, 

QppPn 9p.K.Pp 1 
QppPn 9p,K.Pp 1 

gppMy-y - gupM^y + Qu-yMpp - gp-yMyp , 

9ppJy^ QupJ p'y ~t~ 9y^J pp Qp'yJup • 



(2.31) 

(2.32) 
(2.33) 

(2.34) 

(2.35) 

(2.36) 

(2.37) 



Equations ( p.32| ), (|2.33|) and ( p.34|) are the standard commutation relations of the 
Poincare group. The equations ( p.35| ), ( |2.36| ) and ( |2.37| ) simply reflect the tensor nature 
of the operators. This set also corresponds to the Poincare group transformations. The 
derivation of these equations is explained in ref. 

We have thus shown that under the condition ( p. 17 ) the problem of constructing 
states with definite angular momentum can be formulated in terms of the kinematical 
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operators M^^, and W'^ in exactly the same way as in terms of the operators J^^ and 
5*^ which depend on the interaction. This naturally reflects the fact that the angular 
momentum of a system determines the kinematical properties of the wave function relative 
to transformations of the coordinate system. The dynamics is involved in the composition 
of the angular momentum of the system from the spins of its constituents. 

As already mentioned, the angular momentum generators in ordinary LFD contain 
the interaction. To avoid any misunderstanding, note that the interaction Hamiltonian 
does not disappear in our approach. It is moved into the angular condition ( p.31| ). How- 



ever, it is very convenient to first construct on purely kinematical grounds the general 
form of the light-front wave function for a given angular momentum. One can then find 
from dynamics the coefficients in front of all the spin structures. This way of separating 
dynamical from kinematical transformation properties is of particular interest since the 
interaction Hamiltonian is often approximate. 

We emphasize that by introducing the operator L^^{u), eq.( |2.18| ), containing the 
derivatives over uj, we enlarge the Hilbert space where the state vector is defined. It is 
now the direct sum of the Hilbert space where the "normal" field-theoretical operators act 
and of the Hilbert space, where the operator L^^{uj) acts. Hence, in this enlarged space 
the scalar product and, correspondingly, the normalization condition contains integration 
over uj with the appropriate measure dn^: 

,{\',p'\p,\)^dfi^ = 2po S^^\p-p') S^'^ . (2.38) 

The particular form of the measure dfi^j corresponds to integration over the directions of 
cD in a particular system of reference. The angular condition ( |2.17D , ( p. 311 ) just ensures 



the equivalence of the approach developed in the enlarged Hilbert space to the ordinary 



approach. In particular, the orthogonality condition (|2.38|) , which contains integration 
over UJ, has to be equal to the orthogonality condition (|2.26|) , where u is a fixed pa- 
rameter. This means, of course, that the product ^{\' ,p'\p, X)^^ does not depend on u 
at all, (though for an arbitrary vector | • • •)a; from the enlarged Hilbert space, the scalar 
product u){- ■ ■ \ ■ ■ ■)ui depends on u). However, any separate contribution of the Fock com- 
ponent to ui{X',p'\p, A)^ depends on u, while the ti;-dependence disappears in the sum over 
all Fock components. 

In the enlarged Hilbert space, the angular momentum operator is given in terms of 
derivatives on the momenta and u and do not contain the interaction. The construction 
of states with definite angular momentum becomes therefore very simple. In practice, it 
is convenient first to solve the kinematical part of the problem - the construction of states 



with definite angular momentum - and then satisfy the angular condition (|2.17|) . 

The transformation properties of a state with definite angular momentum are dis- 
cussed in chapter ^. We will see that they completely determine the structure of the 
wave function, and, in particular, the number of its spin components (six in the case of 
the deuteron). These components will be found below by solving the light- front general- 
ization of the Schrodinger equation. At first glance, this procedure is unambiguous and 
nothing remains to be determined by the angular condition. What is then the role of this 
condition? 
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Without solving the angular condition, there is an ambiguity in finding states with 
definite angular momentum. To show that, let us construct the operator: 



A = {W-uf, (2.39) 

where VF^ is the kinematical Pauli-Lubansky vector (|2.28| ). It is readily verified that 
this operator commutes with P^, M^^, W^j^y and also with the parity operator. It seems 
therefore that the state vector must be characterized not only by its mass, momentum 
and angular momentum, but also by the eigenvalue a of the operator A: 

i 0„ = a 0„ . (2.40) 

One can show that, for instance, for the total angular momentum J = 1, there are only 
two states with a = and a = 1. However, these states (pa are degenerate. Since the 
commutators of M^y and of J^j, with P^, J^^ are equal to each other (see eqs.( |2.32| )-( |2.37|) ), 
the operator 

AJ^, = M^, - J^, = L^,{lo) - JiZ' (2.41) 

commutes with and with J^^, but [AJ^,^, A] 0. The state 0' = AJ^^, (p^ is therefore 
not an eigenvector of the operator A, i.e., it can be represented in the form 0' = J2a Pa (pa- 
But it corresponds to the same mass as 0q,. We would thus always get for the deuteron and 
for any state with J = 1 two degenerate states with a = and 1, in evident contradiction 
with reality. The angular condition ( |2.17D is just distinguishing a definite superposition 
of states (pa, 

^ = Y.Cc.(i)a, (2.42) 

a 

which is such that AJnyCj) = 0. This equation eliminates the problem of the "spurious" 



states of relativistic composite systems. This procedure will be illustrated in section |T5 
by a simple example. 



2.2 Covariant light-front graph technique 

The light-front graph technique is a method for calculating the S'-matrix. In the frame- 
work of perturbation theory, the on-shell amplitude given by this graph technique co- 
incides with the one given by the Feynman graph technique. However, the methods to 
calculate them drastically differ from each other, as we shall see in this section. The 
most important difference lies in the fact that in the LFD all four-momenta are always 
on the mass shell. This three-dimensional form of the theory has enormous advantages in 
solving several problems. In applications to relativistic few-body systems, it provides a di- 
rect and close connection between the relativistic wave functions and the non-relativistic 
ones, as we shall detail in the next chapters. Moreover, the diagrams corresponding 
to vacuum fluctuations are absent. This simplifles very much the theory as compared 
to other three-dimensional approaches (like the old fashioned time ordered perturbation 
theory). Another important advantage results from the explicitly covariant formulation 
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of the light-front plane defined by uj-x = as compared to the standard formulation on 
the plane t + 2; = 0. 

We shall derive in this section the rules pertinent to the covariant light-front graph 
technique by transforming the standard expression for the S'-matrix. We explicitly show in 
appendix ^ the relations between the covariant light-front amplitudes and the amplitudes 
given by the Weinberg and Feynman rules. 



2.2.1 General derivation 



The graph technique described below was developed by Kadyshevsky |^ (see ref. for 
a review) and applied to LFD in ref. [^]. It is manifestly covariant, like the Feynman graph 
technique, and retains all the positive features of the old fashioned perturbation theory 
developed by Weinberg ||^. Following ref. p5[, we start from the standard expression for 
the S'-matrix: 



S 



Texp 



W\x)d''x 



1 + 



(2.43) 

where H^'^^{x) is the interaction Hamiltonian. The sign of the T-product (and l/n\) 
are omitted, since the time ordering is made explicit by means of the 6'-functions. The 
expression (|2.43|) is then represented in terms of the light-front time a = uj-x: 



S 



Y: I i-trHl-\x,)9 {cuix, - X,)) Wf{x,) ...9 (c^-(x„_i - x„)) K^^Xn) 



(2.44) 



The index u at i?^"* indicates that if*"* and -ff^"* may differ from each other in order to 
provide the equivalence between ( p.43|) and ( p.44|) . The region in which this can happen 



is a line on the light cone. Indeed, if (xi — X2Y > 0, the signs of u!-{xi — X2) and ti — t2 
are the same and hence = if*"'*. If {xi — x-i)^ < 0, the operators commute: 



[W^\X,),W''\X2)]=0 



(2.45) 



and their relative order has no significance. On the light cone, i.e. if [xi — X2Y = 0, 
LJ-{xi — X2) can be equal to zero while ti — t2 may be different from zero. If the integrand 
has no singularity at (xi — ^2)^ = 0, this line does not contribute to the integral over the 
volume d'^x. However, if the integrand is singular, some care is needed. To eliminate the 
influence of this region on the S'-matrix, we have introduced in ( p.44|) a new Hamiltonian 
H^*, such that expressions ( [^.43D and ( |2.44| ) be equal to each other. The form of H^'^, 
which provides the equivalence between ( 2.43|) and ( 2.44|) , depends on the singularity of 



the commutator ( |2.45| ) at (xi — X2) = 0. For the scalar fields, the singularity is weak 
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enough, and the expressions ( |2.43| ) and ( p.44| ) are the same, so that = if™*. For 



fields with spins 1/2 and 1 or with derivative couphngs, the equivalence is obtained with 
differing from if*"* by an additional contribution (counter term) leading to contact 

terms in the propagators (or so called instantaneous interaction). We shall come back to 

this point later on in this section. 

Introducing the Fourier transform of the Hamiltonian: 

H^{p) = J iij*(x)exp(-2p-x)d^x , (2.46) 

and using the integral representation for the 6 function: 

^/ / NN 1 f+°° exp (iTU!-(xi — Xo)) , 

we can transform the expression ( p.44| ) to the form: 
S = l + i?(0) 



J 2m[Ti — It) 2m[Tn-i — le) 



n>2 



n-lJ 



(2.48) 



The T variable appears here as an auxiliary variable, as defined in eq.( p.47| ); ut has the 
dimension of a momentum. 



The S'-matrix (|2.48| ) gives the state vector 0((t) = S{a) 0o for asymptotic states. 



i.e. for an infinite value of the light-front time a = oo. It determines the on-energy-shell 
amplitude. The off-energy shell amplitude is determined by S{cr) at finite cr. Similarly, 
introducing another 6'-function 6{a — u-xi) in (|2.44|) , one can find easily that the S'-matrix 
on a finite light-front plane a is represented in the form: 

S{a) = 1+ r -^^^^^:^R{ujT)dT , (2.49) 

where 

R{ujt) = V(-i)" / H^{ujT-ujTi) .f'^^ . H^{ujTi-ujT2) . . . ^ ^^""^ ii^(u;r„_i) . 
„ J 2m[Ti — le) 2m[Tn-i — le) 

(2.50) 

These formulae give the iterative solution of the equation (p.l3|) . The matrix elements of 
the operator R{ujt) for r 7^ correspond to the off-energy shell amplitudes, whereas at 
r = they give the on-energy shell amplitude. We shall precise the difference between 
off-energy shell and off- mass-shell amplitudes below in sect. p.3.1| . At a — 00, the only 



residue in ( f^.49D for r = ^ survives, and we recover eq.( |2.48| ). 



We emphasize that despite the presence of the four- vector uj in eq. (|2.48| ), the S'- 



matrix and any physical amplitude do not depend on u;, since eq.( ^.48D gives the same 
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S"- matrix, as the initial one given by eq.( |2.43| ). Similarly, the off-shell matrix R{uot) in 
eq. (|2.50|) , depends on uj and off-shell light-front amplitude does not coincide with the 
Feynman one. This is natural since R{ujt) determines the S-matrix at a given light-front 
plane uj-x = a in the interaction region. As we shall see in chapter |^, the same is also true 
for light-front wave functions. The latters depend on uj since they are always off-shell 
objects and are defined not as asymptotic states, but at any finite time (at any given 
light-front plane in our case). 



2.2.2 Spin system 

The covariant light-front graph technique arises when, as usual, one represents the expres- 
sion ( p.48| ) in normal form. Let us consider for example the simple case of an interaction 
Hamiltonian of the form H = —g(f^{x), where </9 is a scalar field. We introduce the Fourier 
transform ip of the field ip given by: 



(27r)3/2 

1 

(27r)3/2 



ip{k) exp{ik-x)d'^k 
a{k) exp{—ik-x) + a\k) exp{ik-x] 



d^k 



(2.51) 



We thus have: 

<^(A;) = [a{-k)e{-ko) + a\k)e{ko)]y/2rk6{k^ - m^) . (2.52) 
When the S'-matrix ( p.48| ) is reduced to normal form, we obtain the contractions: 

2(J^{p) = m^ip)- ■■<f{kMp)- = 0iPo)Sip' - m^)6^^\p + k) . (2.53) 

It is convenient to replace in the following 6{po) in the propagator (|2.53|) by 9{uj-p). This 
is always possible, since = m? > 0. 

We would like to emphasize at this point that the propagator ( p.53| ) contains the 



delta- function S{p 



m 



and therefore all particles are always on their mass shells. 



The reason of this property, which drastically differs from the Feynman approach, is the 
absence of the T-product operator in (|2.43|) . 

In this graph technique, the four-vectors utj in ( |2.48|) are associated with a fictitious 
particle - called spurion - and the factors 1 / [tj — ie) are interpreted as the propagator 
of the spurious responsible for taking the intermediate states off the energy shell. This 
spurion should be interpreted as a convenient tool in order to take into account off-energy 
shell effects in the covariant formulation of LFD (in the absence of off- mass shell effects), 
and not as a physical particle. It is absent, by definition, in all asymptotic, on-energy 
shell, states. We shall show below on simple examples how the spurion should be used in 
practical calculations. 
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The general invariant amplitude M„m of a transition m — n is related to the S- 
matrix by: 



u=l i=l 



(27r)32£fc. . . . (27r)32£fc, {27r)^2ek, . . . (27r)32£fc„ 



1/2 • 



(2.54) 

where, e.g., = \J mf + kf. The cross-section of the process 1 + 2— i>3 + ...-|-nis thus 
expressed as: 



(27r)4 ,..,2 ^^^^3 

da = — — —\M\^ ^ 



(Pk 



'^6^^\k, + k2 - ks - . . . - kr, 



where j is the flux density of the incident particles: 

jek.Sk^ = ^[s - (mi + ma) ^/^[s - (mi - m2)Y^^ s = (fci + fca)^ • 



(2.55) 




, COT2 

Fig. 2. The vacuum vertices. 




ill 



To find the matrix element M of order n one must proceed as follows |3^, ^ 



1. Arbitrary label by a number the vertices in the Feynman graph of order n. Orientate 
continuous lines (the lines of physical particles) in the direction from the smaller to 
the larger number. Initial particles are oriented as incoming into a graph, and final 
particles as outgoing. Connect by a directed dashed line (the spurion line) the ver- 
tices in the order of decreasing numbers. Diagrams in which there are vertices with 
all incoming or outgoing particle lines (vacuum vertices, as indicated in fig. |^) can 
be omitted. Associate with each continuous line a corresponding four-momentum, 
and with each j-th spurion line a four-momentum utj. 



^In order to stick to conventional notations, the normalization of the amplitude given by the standard 
formulae (2.54), (2.55) and some factors in the rules of the graph technique (mainly, the degrees of 2tt) 
differ from the ones previously used in ref. Besides, the numbering of vertices here and in 

^ ^ has opposite order. 



20 



2. To each internal continuous line with four-momentum k, associate the propagator 
9{uj-k)5{k'^ — m?), and to each internal dashed line with four-momentum ujTj the 
factor l/(rj — ie). 

3. Associate with each vertex the coupling constant g. All the four-momenta at the 
vertex, including the spurion momenta, satisfy the conservation law, i.e., the sum 
of incoming momenta is equal to the sum of outgoing momenta. 

4. Integrate (with d'^k/{2'KY) over those four-momenta of the internal particles which 
remain unfixed after taking into account the conservation laws, and over all Tj for 
the spurion lines from — oo to oo. 

5. Repeat the procedure described in |l|-^ for all n\ possible numberings of the vertices. 

We omit here the factorial factors that arise from the identity of the particles and depend 
on the particular theory. 

The vacuum vertices indicated in fig. ^ disappear for a trivial reason: it is impossible 
to satisfy the four-momentum conservation law for them. Indeed, the conservation law 
for the vertex of fig. |] has the form ki + k2 + k^ = uj{ti — T2). Since the four-momenta are 
on the mass shell: kl_^ = m? > 0, so that the left-hand side is always strictly positive: 
[ki + k2 + k^Y — 8m^, whereas the right-hand side is zero since = 0. However, it will 
be seen that the vacuum contributions that vanish in the light-front approach leave their 
track in a different way, making in the cases discussed below the light-front interaction 
Hi^{x) in eq.( |2.44| ) different from the usual interaction H{x) in ( [^.43| ). 



The light-front diagrams can be interpreted as time-ordered graphs. As soon as 
the vertices are labelled by numbers, any deformation of a diagram changing the relative 
position of the vertex projections on the "time direction" from left to right does not 
change the topology of the diagram and the corresponding amplitude. Therefore it is 
often convenient to deform the diagram so that the vertices with successively increasing 
numbers are disposed from left to right. This just corresponds to time ordered graphs. 

The light-front amplitudes can be also obtained by direct transformation of a given 
Feynman amplitude. This transformation is given in section p.3| of appendix El 



2.2.3 Spin 1/2 system 

The rules of the graph technique for spin 1/2 particles are similar to those given above 
except for the fact that one has to worry about contact interactions we already mentioned 
in section ( p.2.1| ). To see this from a more practical point of view, let us first consider the 
diagrams of figs, ^(a) and (b) for scalar particles and consider for a moment the space-like 
plane X-x = a with = 1, as developed by Kadyshevsky p5| . 

In accordance with the rules given above, the amplitude for the diagram of fig. ^(a) 
gets the form: 



Ma = g^ S (k + p + XnY 



m 



e[Xik + p)+n] 



dri 



(2.56) 
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p' 

Fig.3. Exchange of a particle in the s channel; (a) contribution from intermediate state 
containing one particle; (b) contribution from antiparticle intermediate state; (c) same 
as (b) but rewritten as a time ordered diagram 
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Integrating over dri by means of the (5-function, we get: 



9' 



2TiJ[X{k + p)]^ + ni? - {k + pf 



Ma 

Ti = -\\k + p) + ^ [\-{k + p)f + m? - {k + pf . 
Similarly, we obtain the following expression for the amplitude of fig. ^(b): 



(2.57) 



2T2J[Hk + p)f + m? - {k + py 



M, = 

T2 = Xik + p) + ^J[\ik + p)Y + m'^ - (k + p)^ . 
The sum of ( p.57|) and (|2.58| ) gives the usual Feynman amplitude in the s channel: 



(2.58) 



m'' — [k + pY 



In the case where A = (1, 0), we recover the contribution of the usual time ordered graphs 
in the old fashioned perturbation theory: 



Mh 



2e 



k+p 



^k+p ^k 



2e 



k+p 



^k+p + ^fc + ^1 



(2.59) 



The diagram of fig. ||(b) and the amplitude Mf, correspond to a vacuum contribution. 
This is clear if we draw it as a time ordered diagram, as it is done in fig. ^(c). 

The light-front case can be obtained by introducing first the four vector u with 
5^, replacing in the above formulae A by uj/6 and then taking the limit 6^0 



UJ 



I.e. 



— > 0. This limit is however a delicate issue. For the scalar case under consideration 
above, it can be checked that the amplitude M^ disappears. In this limit, the system 
of reference where A = (1,0) moves with a velocity close to c. The graph technique for 
a;^ = is therefore naturally related to the old fashioned time ordered perturbation theory 
in the infinite momentum frame, corresponding to the particular case u = (1,0,0,-1), 
that defines the "standard" light front t + z = 0. 

The full consideration of spin 1/2 particles shows that the above result has to be 
completed. In such a case, indeed, one should associate to each propagator the factor 
(k + m) for a fermion and (m — k) for an antifermion, where k = A;^7^ (see for instance 
section |B.3.2| in appendix 0) . 

The amplitude Ma, in eq.( |2.57| ), is thus multiplied by (pi+m) = {p + k + Xti + m) 
and Mfe in eq. (|2.58| ), is multiplied by (m — P2) = — (Ar2 — p — k — m) (for shortness we 
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omit the initial and final spinors). Replacing again \hy uo/ 5 with uj"^ 
limit 5 — >^ 0: 



n k + f) + m n 



— {k + py 



2uj{k+p) 



2u-{k + p) 



(5^, we get in the 

(2.60) 
(2.61) 



The amplitude (|2.61|) is the contact term we already mentioned in section |2.2.1| . It is 



just the track of the disappeared vacuum diagrams. The sum of the amplitudes Fa and 
gives the Feynman amplitude for spin 1/2 particles. As we shall see later on, contact 
terms are indeed essential in getting fully covariant results (i.e., independent of u). Other 
examples will be considered in the following sections. 




Fig. 4. Contact interaction contributing to the exchange of a spin 1/2 particle in the s 
channel 



To incorporate the contact term from the very beginning, one should not consider 
the diagram of fig. H(b), but add to the diagram of fig. H(a) the diagram indicated in 
fig. 1^, which is obtained from fig. |^ by deleting the spurion line and marking the internal 
line by a cross. This contribution can be also derived from a counter term added to if*"* 



in order to get H^^, as it was done for spin 1/2 in ref. ||38|. Associating the crossed line 
to a fermion carrying the momentum /, we assign to this line the factor —uj9{uj-1)/{2uj-1). 
For the diagram of fig. | this rule gives eq. ( |2.61|) back. 



We can thus formulate the rules of the covariant light-front graph technique for the 
case of spin 1/2 particles. 

1. Transform the Feynman graph of a given order to the set of the light-front graphs in 
the same way as described in point 1 of the rules for scalar particles. Without making 
change in the orientation, replace the ordinary lines corresponding to antifermions 
by double lines, so that the number of fermions minus the number of antifermions 
is conserved. Both initial fermions and antifermions are shown by lines incoming to 
a diagram, and final fermions and antifermions by outgoing lines. 

2. Consider the diagrams with internal fermion and antifermion lines labelled at their 
extremities by two successive numbers (i.e., the ends of the lines are directly con- 
nected by the spurion line, see, e.g., fig. In addition to all the diagrams we 
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already have, we create from the latter diagrams another set of diagrams, by delet- 
ing those spurion lines, which connect the ends of the fermion (and antifermion) 
lines. Put on the corresponding fermion and antifermion lines a cross (see, e.g., 
fig. Draw any diagram with incoming lines at the left and outgoing lines at the 
right. 

3. The analytical expression for the amplitude is written from the left to the right. 
The factors in this expression are written in the order where they are met when one 
goes through a diagram, starting at the right, from an outgoing fermion line, and 
continuing in the direction opposite to the orientation of the fermion lines. For an 
incoming antifermion line, one starts at the left, and continues in the direction of 
orientation of the (double) antifermion lines. 

4. To each internal continuous line with four-momentum k, associate the propagator 
{k + m)6{uj-k)6{k'^ — iv?) for a fermion, and the propagator (m — k)9{uj-k)5{k'^ — m?) 
for an antifermion, the factor —Co9{uj-k) / {2uj-k) for each crossed fermion line, the 
factor uj6{uj-k) /{2uj-k) for each crossed antifermion line, and the factor l/{Tj — ie) 
for each internal dashed line with four-momentum utj. 

5. Associate with each vertex the factor gV, where V depends on the type of coupling 
(1, ^75, 7^, etc.) in the Hamiltonian H = —gipVip---. All the four-momenta at 
each vertex, including the spurion momenta, satisfy the conservation law, i.e., the 
sum of incoming momenta is equal to the sum of outgoing momenta (incoming and 
outgoing momenta always correspond to the incoming and outgoing lines). 

6. With any outgoing fermion [antifermion] line with momentum associate the spinor 
u{p) [f (p)] and with any incoming fermion [antifermion] line associate the spinor u{p) 
[v{p)]. 

7. Integrate (with (i'^A;/(27r)^) over those four-momenta of the internal particles which 
remain unfixed after taking into account the conservation laws, and over all Tj for 
the spurion lines from — oo to oo. 

8. Repeat the procedure described in ^[|^ for all n\ possible numberings of the ver- 
tices. One should also take into account the standard sign factors appearing from 
permutation of fermion fields. 

In the case of a diagram containing boson lines connected to fermion loops, one 
should go through the outgoing boson line, and then, having reached the fermion line, 
pass it as indicated in the rules, i.e., in the direction opposite to the orientation of the 
fermion line and along the orientation of the antifermion line. Both lines in the loop, 
single and double, are oriented in the same direction, but one makes a cycle when passing 
through them. 

We emphasize that the fermion and antifermion are distinguished by the type of the 
continuous lines (single or double). Both the fermion and antifermion lines are directed 
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from the smaller to the larger number (in the time ordered graph they both propagate 
from left to right, i.e., in the direction of time increase in time ordered graph). In contrast 
to the standard versions of the Feynman rules, any incoming line, both for antifermion 
and fermion, corresponds to the initial state, and outgoing lines - to the final state. 
The fermion and antifermion lines are followed in opposite directions (opposite to the 
orientation for a fermion line and along the orientation for an antifermion line). The 
fermion and antifermion propagators differ from each other. 

In many practical applications, the contact term (or so called instantaneous inter- 
action) can be conveniently incorporated for spin 1/2 particles by replacing the spin part 
of the propagator (m ± k) (for those lines where the contact term contributes at all) by 
[m± {k — tbr)] (however, the delta-function in the propagator 5(A;^ — w?) still depends on 
the argument (/c^ — m^)). Together with the fact that the contact term originates from the 
lines which ends are labelled by two successive numbers, this rule coincides with the rule 
given in refs. [^, the contact terms modify only the propagators corresponding to 



the "lines extending over a single time interval". Two successive numbers just determine 
a single time interval. 

The above replacement to incorporate systematically the contact term can always 
be made in perturbation theory. However, it cannot be made in more complicated cases. 
In particular, the graph for the wave function, for instance, indicated in fig. ^ below, 
does not contain any external crossed lines. Hence, the contribution of the contact terms 
to the deuteron electromagnetic form factors or electrodisintegration amplitude, in the 
impulse approximation for example, cannot be incorporated by the above substitution 
in the propagators, since it would create an object with an external crossed line. The 
contact terms in the deuteron electromagnetic form factors and electrodisintegration will 
be properly taken into account in section |7.2.1| . 

If a diagram contains two fermions (but not fermion and antifermion) turning into 
a boson, like in the case of the deuteron form factor (the vertex NNd), one should pass 
through one fermion line from the right to the left, and then start again from the same 
final vertex and pass through the second line. Both lines are followed in the direction 



opposite to their orientation. Such example is given below in section |6.3.3 . 



2.2.4 Spin 1 and coupling with derivatives 

For particles with spin 1 the rules are similar to the case of particles with spin 1/2. Like in 
Feynman rules, external lines for spin 1 particles are associated with polarization vectors. 
However, in the case of spin 1 particles, there is no difference between propagators of 
particle and antiparticle. The propagator has the form 

D,,{k, m) = (-g,, + ^] e{uj-k)6{k' - m') . (2.62) 



m 



The contact term is given by: 



A,. = + - I7^^(^' - • (2-63) 
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Like in the case of spin 1/2, the contact term contributes only in the hnes extending over 
a single time interval. It can be incorporated simply by replacing the spin part of the 
corresponding propagators {Qf^u — k^ky/rn?) by [g^^ — (A; — ujT)^{k — ijjr\lrv?\. 

In the case of massless vector boson (the photon for instance) the form of the prop- 
agator, as usual, depends on the gauge. For a general gauge, the momentum dependence 
of the gauge term in the propagator induces also a corresponding contact term. We do 
not investigate here the general case and give the photon propagator for the Feynman 
gauge: 

D^,(fc, m = 0) = -g,,,Q[i^-k)8{k^) . (2.64) 

There is no contact term in this gauge. 

It is sometimes more convenient to use in this formulation of LFD the light-cone 
gauge defined by lo-A = 0. The spin-1 propagator for a massless particle in this gauge is 

D,Ak, m) = (^-g,. + ^^E^^^l^^ e{u;-k)6{e - m') (2.65) 
and the corresponding contact term has the form: 

= ■ (2.66) 



The contact term ( 2.66 ) can be again incorporated by the replacement k ^ k — lot in the 
spin part of eq.( p.65| ). 

Consider finally the coupling with a derivative, e.g., when the interaction Lagrangian 
contains a term like d^(f{x). Since all the lines are oriented, any line is either incoming 
in the vertex, or outgoing from it. Let this line be associated with the four-momentum 
k. The derivative results in the multiplication of a vertex by the extra factor i/c^ for an 
outgoing line, and by the factor —ik^ for an incoming line. The contact term contributes 
to the lines involving the derivative of a field and having the same topology, as in the case 
of spins 1/2 and 1. That means that the ends of these lines should be connected by the 
spurion line. Like in the case of spin 1 propagator, the contact term can be incorporated 
by the replacement fc^ — >■ /c^^ — a;^r, where cur is the spurion momentum, and this spurion 
line should connect the ends of the particle line, corresponding to the field with derivative 
coupling. 



2.2.5 Ultraviolet and infrared behavior 

A peculiarity of the covariant light-front amplitudes is that they have no any ultraviolet di- 
vergences for finite values of all the spurion four-momenta. All the ultraviolet divergences 
in light-front diagrams appear after integrations over tj with infinite limits ||3^. Indeed, 
the energy-momentum conservation (including the spurion four-momentum) is valid at 
any vertex. Since all the four-momenta are on the corresponding mass shells, we have at 
each vertex a real physical process as far as the kinematics is concerned. For finite initial 
particle energies and for finite incoming spurion energy, the energies of the particles in the 
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intermediate states are thus also finite. Hence, the integrations over the particle momenta 
for fixed spurion momenta are constrained by a kinematically allowed finite domain. It 
is the same reason that provides finite imaginary part of a Feynman diagram found by 

replacing the Feynman propagators 77-5 ^ r by the delta-functions —mSik"^ — m?). 

[k — m + le) 

In both cases the internal particle lines are associated with the delta-functions. 

The only source of the ultraviolet divergences in the light-front amplitudes is the 
infinite intermediate spurion energies, i.e., infinite r,-. This is the reason why divergences 
may appear at the upper limit of integration over Tj. Since Tj are scalar quantities, one 
can introduce an invariant cutoff in terms of these variables. This way of regularizing 
the divergent diagrams is another advantage of the covariant formulation of LFD. We 



illustrate this property of the light-front amplitudes in sect. p.3.3| for the example of the 
self-energy. 

For massless particles, light-front amplitudes may have infrared divergences, like in 
the case of Feynman diagrams. 

Another peculiarity of LFD is the appearance of "zero modes" . For constituents of 
zero mass, for instance, the state vector may contain components with uj-k = for non- 
zero four-momentum k. In the standard approach, this corresponds to the finite light-front 
energy k~ = k\/k^ for both fc+ = and k'j_ = 0. Zero modes can also appear in theories 
with spontaneously broken symmetry. They make the equivalence between LFD and the 
instant form of quantization in which nontrivial vacuum structures (condensates) appear 



The detailed discussion of the above problems is beyond the scope of the present 
review. 



2.3 Simple examples 
2.3.1 Exchange in t-channel 
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Consider the diagrams shown in fig. ^. It corresponds to the exchange of a scalar particle 
of mass fi between two scalar particles, in the t channel. These diagrams determine, in 
the ladder approximation, the kernel of the equation for the calculation of the light-front 
wave function. The external spurion lines indicate that the amplitude is off-energy shell. 
According to the light-front graph technique for spinless particles, the amplitude has the 
form: 



/C = J e {uih - k[)) 6 [{h - k[ + LOTi - Lorf - fi^) 

+g^ J e (cu-(fc; -k^))5 ({k[ - fci + cjTi - ut'Y - fi^ 

2 



Ti — te 
dri 

Ti — ie 



g'e{ujik,-k[)) 



/i2 _ (fci - k[y + 2Tuj\ki - k[) - ie 



g^e {u;.{k[ - k,)) 

\? - {k[ - k^f + 2T'uj-{k[ -k^)-ie' ^ ' ' 



The two items in ( |2.67|) correspond to the two diagrams of fig. ^. They cannot be non-zero 



simultaneously. On the energy shell, i.e. for both r = r' = 0, the expression for the kernel 
is identical to the Feynman amplitude: 

/C(r = r' = 0) = , ^\,.^ . . (2.68) 
li'^ — [ki — k'lY — le 



Note that the off-shell amplitude ( 2.67 ) depends on uj. This agrees with the fact that 
it is related by eqs.( p.49"D , (|2.5CI|) to the ^-matrix on a finite light-front plane in the 
interaction region. It depends therefore on this plane. In the case of a scalar amplitude, 
this dependence is given in terms of extra scalar variables (in addition to Mandelstam 
variables s and t). These variables are the scalar products of uj with the four-momenta, as 
seen from eq. ( |2.67| ). Hence, the amplitude has extra singularities as a function of this new 
variables, which can be found similarly to singularities of the Feynman amplitudes |^ . 



On the energy shell, corresponding to r = r' = and to the light-front plane 
shifted to infinity, out of the interaction region, the dependence of the amplitude on uj 
disappears. In more complicated cases, when a Feynman diagram corresponds to the sum 
a few light-front diagrams (like in the case of the box diagrams considered in appendix |B|) , 
the amplitude for a particular light-front diagram may depend on uj even on the energy 
shell. This dependence disappears in the sum of all amplitudes at a given order. In this 
case the singularities of different amplitudes cancel each other in the sum. 

The off-energy shell amplitude corresponds to a diagram with external spurion lines. 
This term - off-energy shell - is borrowed from the old fashioned perturbation theory. As 
mentioned above, the latter is obtained, if one introduces in eq. (|2.44| ), instead of uj, the 
four- vector A with A = (1, 0, 0, 0). The difference of initial and final four-momenta differs 
from zero by Ar. Due to A = the sums of initial and final spatial components of 
momenta is still equal to each other, but for r 7^ initial and final energies are not equal 
to each other (their difference is equal to r). This just characterizes the "off-energy shell" 
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amplitude, which can be a part of a bigger diagram. In the hght-front amphtude for 
arbitrary cu, all the four- vector components do not satisfy the conservation law. Their 
difference is proportional to ujt. Hence the term "off-energy shell" may be inappropriate 
for the light-front amplitude, since not only energies, but also three-momenta are not 
conserved. We emphasize again that the four-momenta are always on the mass shells, 
whereas the Feynman amplitude depends on the four-momenta which may be off-mass 
shell, but the conservation law is always valid. 

2.3.2 Compton scattering 




1 , 2 

Pi=p+k+(ax 






(a) 



(b) 



Fig. 6. (a) The direct diagram for electron Compton scattering with an electron in the 
intermediate state, (b) The contact term corresponding to the diagram (a). 

The Compton scattering amplitude from an electron (through s-channel) is indicated in 
fig. IKa), and, for the contact term, in fig.^b. According to the above rules of graph 
technique, we associate with fig. ||(a) the following expression: 

u'{p')e'*{m + pi)euip)6{ijj-pi)5{p\ — m^) 6^'^\p + k + ujt — pi)d^pi 



g\'{p')e'* T^fltM P) + 9'n'{p')e'* 
m'^ — {p + k) 



ijj 



2uj-{p + k) 



eu{p), 



(2.69) 



where e = 7/^6^ (fc) and similarly for e'*. The polarization vector of the photon is denoted 
by e^{k). The order of factors in ( p.69| ) corresponds to follow the fermion line from the 
right to the left, in the direction opposite to its orientation. 

The amplitude for the contact term contribution, indicated in fig. p|(b), is: 



-g^u'{p)e* 



UJ 



2lj{p + k) 



eu{p) 



(2.70) 



We omit here the theta-function 9{uj-{p + A;)), which is always 1. 

The expressions ( |2.69| ), ( |2.70| ) reproduce the amphtudes ( p.60| ), ( p.61j ). In the sum. 
Ma + Mb, the omega- dependent items cancel each other. 
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Fig. 7. (a) The cross diagram for electron Compton scattering with an electron in the 
intermediate state, (b) The contact term corresponding to the diagram (a). 




(b) 

Fig. 8. (a) The cross diagram for electron Compton scattering with a positron in the 
intermediate state, (b) The contact term corresponding to the diagram (a). 
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Now consider the crossed amplitude, indicated in figs. |^ and |^. The amphtude 
corresponding to fig. |^(a) has the form: 

u'{p')e{m + pi)e'*u{p)6{ijj-pi)5{p\ — m^) -6^'^\p — k' + ujt — pi)d^pi 

T — le 

+g'u'{p')e ^^_^^_ ^^^ e'*u{pWi^iP - k')). (2.71) 

The contact term contribution in fig. |^(b) writes: 

M, = -g'u'{p')e ^^_^^_ ^^^ e'*uip)9iu;ip - A;')). (2.72) 

It has the form of the second term in ( |2.71| ) and cancel it in the sum Ma + M^. 

The amplitude of fig. §(a) differs from fig. |^(a) by the time ordering of the vertices. 
If we draw the vertices in the order of increasing numbers from the left to the right, this 
amplitude appears as a time ordered diagram. It contains the creation of a e"'"e~-pair by 
the photon, i.e., includes a positron in the intermediate state, shown by the double line, 
and its subsequent annihilation. The corresponding amplitude has the form: 

u'(p')e(m — pi)e'*u(p)9{u-pi)5ip\ — rn?) S'^'^Uk' — p + out — pi)d'^pi 

T — ie 

-g'u'{p')e^-j^^^--^e'*u{p)0{ujik' - p)), (2.73) 
and for the contact term indicated in fig. |^(b): 

M, = g^u'ip')e^-^^--^e'*uip)Oi^ik' - p)). (2.74) 



The first item in eq. (|2.71| ) differs from zero in the region u-{p — k') > 0, whereas in 
eg. ( 2.73 ) this occurs in the region uj-{p — k') < 0. The sum of these two contributions 
reproduce the Feynman amplitude in the ^-channel, in the whole kinematical range. 

In order to show how to deal with antifermion states, consider now the Compton 
scattering amplitude from an antifermion. The s-channel diagram is similar to fig. ^ with 
the single lines replaced by the double ones. The corresponding amplitude has the form: 

- ^'^-(P^^^^^Z-'-W - >Mp)i^^fAp% (2.75) 



while the contact term gives: 



Mb = 9Mp)e^^^j;^^^e'*v'{p'). (2.76) 
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The order of factors in ( |2.75| ) and ( |2.76| ) corresponds to follow the antifermion line from 
the left to the right, in the direction along its orientation. Note that in eg .( p .761) the 
sign of the contact term originating from the antifermion line is opposite to the fermion 
contact term, eg .( p .701) . The corresponding cu-dependent contribution in ( |2.75|) has also 
an opposite sign and thus disappears in the sum Ma + M^. 

The amplitude for the cross diagram with an antifermion intermediate state is also 
similar to fig. and is given by: 

v{p)e'*{m — pi)ev' {p')6{ijj-pi)5{p\ — m^) — k' + ujt — pi)d^pi 

-9Mp)e \^_^^_j^,^ ev'{p')e{ujip - k')). (2.77) 

The corresponding contact term is: 

M, = 9Mp)e' \^_^^_j^,^ ev'{p')e{u;ip - k')). (2.78) 

The kinematics of the Compton scattering from fermion and antifermion, i.e., the ex- 
pression of the intermediate momentum pi through the external momenta, as given in 
eg. (|2.77|) , is here identical. Besides the spinors u and v, the difference is in the signs of 
the spin parts of the propagators and of the contacts terms. 

Finally we give the amplitude corresponding to the scattering on an antifermion 
with a fermion intermediate state (same as fig. |^, but with single lines replaced by the 
double ones and vice versa): 

/di~ 
v{p)e*{m + pi)ev' {p')6{uj-pi)5{p\ — m^) -6^^\k' — p + ujt — pi)d'^pi 

+gMp)e \^_^^,_^^ ev'{p')9iu;ik' - p)) (2.79) 

and the contact term: 

M, = -gMp)e' \^_^^, _^^ iv'{p')e{ujik' - p)). (2.80) 

The difference between ( |2.77| ), ( p.78| ) and ( p.79| ), ( |2.80| ) is in the sign of the argument 
of the corresponding theta-functions. The sum of both contributions gives the Feynman 
amplitude for Compton scattering on an antifermion. 

As mentioned above, the contact term can be incorporated everywhere by the fol- 
lowing replacement in the spin matrix of the propagator: pi —>■ pi — lot. 
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p' -co t' + coti - k 



Fig. 9. Self-energy loop. 



2.3.3 Self-energy contributions 

Another simple example is the self-energy diagram shown in fig. |^. The corresponding 
amplitude (equal to the self-energy up to a factor) has the form: 



S(p') =g^ J e{uj-k)6{k^ - m^)e {u\pi + - k)) 5 ((pi + ujt^ - kf - m^) 



dl^k dri 



(27r)^ Ti — ie ' 
(2.81) 

with pi = p' — out'. 

Let q = pi+uJTi. The integral over d^k is thus reduced to the well known calculation 
of the imaginary part of the Feynman amplitude, when all the propagators are replaced 
by the delta-functions: 



J S{k^ - m^)6 ((g - kf - m^) d^k 
Inserted in ( |2.81| ), it gives: 



TT 



2v^ 



g2 _ _ 



(2.82) 



S(pi) 



167r2 



irrfl- 



pI - 4m2 + n dTi 



pI + ri ^1 - 



(2.83) 



The logarithmic divergence is at the upper limit of the integration over ri. This is a 
particular manifestation of a general property of the light-front amplitudes, discussed 
above in sect. |2.2.5|. One can introduce the invariant cutoff in terms of ri. In this way. 



after renormalization, the standard expression for the self-energy amplitude is obtained. 
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Chapter 3 

The two-body wave function 



Since the wave functions of a composite system are the Fock components of the state 
vector, their transformation properties under four dimensional rotations are determined 
by the corresponding properties of the state vector investigated in the previous chapter. 
This allows one to construct their angular momentum explicitly. The angular condition, 
necessary to construct the angular momentum operator, is here of particular relevance. We 
pay also particular attention to the normalization of the wave function. The comparison 
with the BS wave function is made, and a simple example (Wick-Cutkosky model) is 
developed. 



3.1 General properties of the wave function 

The wave functions we consider are the Fock components of the state vector defined on 
the light-front plane uj-x = a. This means that they are coefficients in an expansion of 
the state vector (j)'^^{p) with respect to the basis of free fields: 



2 



X S^^\h + k,-p-.r) exp(.ra)2(..p)rfr- 2.^3/2^ i2n)V^V2^, 
+ (2vr)3/2 J ^/^t^,.,,,,3,3 {k„ k„k„ p, ujT)al (k,)al {k,)al {k,) 1 0) 

X 5'^'^\ki + k2 + k2 — p — ujt) exp{iTa)2{uj-p)dT 
d^ki d^k2 d^ks 

(27r)3/2^/2?^(27r)3/2y2?^(27r)3/2^/2^ + '" " 

Here A in /^■^^{p) is the projection of the total angular momentum of the system on 
the z-axis in the rest frame, where p = and cri,cr2,cr3 are the spin projections of the 
particles 1, 2 and 3 in the corresponding rest systems. 

We emphasize in ( |3.1| ) the presence of the delta-function 5'^'^\ki + ■ ■ ■ + kn — p — 
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ujT)2{uj-p)dT. Formally it can be obtained from the relation ( |2.2CI| ): 

exp(iP°-a)0^((T) = exTp{ip-a)(f)^{a + cu-a) (3.2) 

with 

Indeed, after substituting ( |3.1|) in (|3.2| ), the action of the operator exp(zP°-a) on the 
ra-body sector gives the factor exp [i{ki + ■ ■ ■ + fc„)-a] in the integrand of the l.h.s. of 
( p.2| ), whereas the factor exp [i{p + ujT)-a\ appears in the r.h.s. The delta function in (|3.1|) 
ensures the equality of these factors and, hence, the relation ( p.2|) . 

In the particular case where cu = (1, 0, 0, —1), the delta-function S^^\ki + k2—p — ujT) 
gives, after integration over r, the standard conservation laws for the (_L, +)-components 
of the momenta, but does not constrain the minus-components. 

From ( p.l|) one can see that the wave function depends on the orientation of the 
light front, from its argument ujt. This important property of any Fock component is 
very natural. As explained in the previous chapter, any off-energy shell amplitude is 
related to the S-matrix defined on the finite light-front plane in the interaction region 
and therefore depends on its orientation (see eg. ( p.49|) ). The bound state wave function 
is always an off-shell object (r 7^ due to binding energy). Therefore it also depends on 
the orientation of the light-front plane. This property is not a peculiarity of the covariant 
approach. It allows however to parametrize this dependence explicitly. We will investigate 
below this dependence for a few systems. 



3.1.1 Transformation properties 

We derive in this section the transformation properties of the relativistic wave functions 
under four dimensional rotations. To do this, we write eq.(|2.9|) in a more explicit form: 



4>i\p) - 4>',i\9P) = UAg)<Pi\p) , (3.3) 

where Ujo(g) is given for infinitesimal transformations by eq.( |2.10 ) and g is a Lorentz 
transformation and/or rotation. Here A is the projection of the angular momentum op- 
erator on the z-axis in the system at rest, p = 0. After transformation, the new state 0' 
does not correspond to a definite projection of J. We then expand (p' with respect to the 



X' 



Here Dyl{R{g,p)} is the matrix of the rotational group and R{g,p) is the rotation oper- 
ator: 

R{g,p) = L-\gp)gL{p) , (3.5) 

where L{p) is the Lorentz transformation, corresponding to the velocity v = p/po, i.e., for 
example, 

L{K){m, 0,0,0) = {mKo/^/K^,mK/^/K^). 
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The Euler angles that determine the rotation R{g,p) can be expressed in terms of the 
momentum p and the parameters of the transformation g. The exphcit expression of the 
Euler angles in terms of g and p will not be needed. 

To obtain the transformation properties of the state vector, we first substitute 
from ( ^.41 ) in ( |3.3D , and represent in the form of the expansion ( |3.1D . Since aj-|0) in ( |3.1| ) 



is the state vector of a free particle with spin j, the operator aj^ transforms in accordance 
with the law 

UAg)aUk)Uj,\g) =Y.D^lW9.k)}aUgk) . (3.6) 

Comparing the left- and right-hand sides of the resulting equation, we thus obtain for the 
two-body wave function: 

^il,,,Agk,.gk,,gp,gujT) = Y: 

A'o-; cr' 



^K%wS^^^k2,p,UT) (3.7) 

and similarly for any n-body Fock component. Here X,ai,a2 are the projections of the 
spins on the z axis in the rest frame of each of the particles. We emphasize that the 
transformation ( ^.71 ) is a purely kinematical one. As mentioned above, in the usual for- 
mulation of LED, the transformation of the system of reference changes the position of 
the light-front plane. Indeed, under the transformation x ^ x' = gx, the state vector is 
transformed from the plane t + z = in a system A to the plane t' + z' = in a. sys- 
tem A'. Introducing cu^^^ = (1,0,0, —1) we indeed obtain different planes co^^^-x = and 
= {g ^u;(°))-a; = 0. In the covariant formulation of LED, the state vector remains 
to be defined on one and the same plane u-x = uj'-x' = in both systems A and A'. 
This property ensures the kinematical transformation law (|3.7|). The dependence on the 



surface is then given by the dynamical dependence of the wave function on uj. 
3.1.2 Parametrization in the spinless case 

We will mainly concentrate on the two-body wave function. Generalization to the n-body 
case is straightforward. 

Due to the conservation law 

ki + k2 = p + uJT , (3.8) 
the light-front wave function can be shown graphically like a two-body scattering ampli- 



tude as indicated in fig. |I0[ The broken line corresponds to the fictitious spurion. We 
emphasize again that although we assign a momentum ujt to the spurion, there is no any 
fictitious particle in the physical state vector. The basis in eq. (|3.1|) contains the particle 
states only. Due to this analogy, the decomposition of the wave function in independent 
spin structures and their parametrization is analogous to the expansion of a two-body 
amplitude in terms of invariant amplitudes. We will use this analogy below. 
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Fig. 10. Graphical representation of the two-body wave function on the hght front. The 
broken hne corresponds to the spurion (see text). 



Let us consider a system of two spinless particles in a state with zero angular momen- 
tum. According to ( p.7| ), its wave function ip = ^j^=j^=j=Q is scalar and, hence, depends 
on invariant variables. We can first use the Mandelstam variables: 

s={h + k2f = {p + ujt)\ t = {p-k^)\ u = {p-k2)\ (3.9) 

with 

s + t + M = + 2m2 . (3.10) 

The wave function depends therefore on two independent scalar variables: = ipisyt). It 
will be convenient in the following to introduce another pair of variables, the relativistic 
momentum k which corresponds, in the c.m.-system where ki + k2 = 0, to the usual 
relative momentum between the two particles. Note that this choice of variable does not 
assume, however, that we restrict ourselves to this particular reference frame. We denote 
by n the unit vector in the direction of uj in this system. Note that due to the conservation 



law ( ^.Sj) , the total momentum p ^ of the system in this reference frame is not zero. In 
terms of these variables, the wave function takes a form close to the non-relativistic case. 
Making the appropriate Lorentz transformation, we get: 

k = L-\V)h = h- -^[kio - J-'"^ ] , (3.11) 



n = L-^{V)Q/\L-\V)Q\ = Vv^ L-\V)Q / u-p , (3.12) 

where 

V = p + ur, (3.13) 

and L^^iV) is the Lorentz boost. The unit vector n is reminiscent of the unit vector p/\p\ 
which appears in the infinite momentum frame. 

From these definitions, it follows that under a rotation and a Lorentz transformation 
g of the four-vectors from which k and n are formed, the vectors k and n undergo only 
rotations: 

k' = R{g, V)k, n' = R{g, V) n , (3.14) 
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where R is the rotation operator ( |3.5| ). Therefore k"^ and fi-k are invariants and can be 
expressed in terms of s and t. For the wave function with zero angular momentum we 
thus obtain 

i) = ip{k,n)=ip{k'^,n-k). (3.15) 

It is seen from ( p.l5| ) that the relativistic hght-front wave function depends not only on 
the relative momentum k but on another variable - the unit vector n. 

Finally, we introduce the third set of variables in which the wave function can be 
parametrized, in analogy to the equal-time wave function in the infinite momentum frame. 
We define the variables: 

X = uj-ki/u-p , Rl = ki — xp , (3.16) 

and represent the spatial part of i? as i? = R\\ + -Rj_, where R\\ is parallel to to and R± 
is orthogonal to u. Since R-u = Roujq — Rw-Q = by definition of i?, it follows that 
-Ro = \R\\\) and, hence, R\ = —R^ is invariant. Therefore, the two scalars on which the 
wave function should depend on can be R\ and x: 



ij = i,{Rl,x) . 



(3.17) 



Using the definitions of the variables R\ and x, we can readily relate them to k'^ and n-k: 



[n-kf 



X 




The inverse relations are 



- 2 R\+rTi? 2 
/c = -4 r - m\ 



4x(l — x) 



n-k 



R\ + m? 
xil — x) 



1/2 



(3.18) 



1 

X 

2 



(3.19) 



The variables introduced above can be easily generalized to the case of different masses and 
an arbitrary number of particles [|]. The corresponding variables g^, n are still constructed 
according to eqs.( p.llD , ( p.l2|) and the variables Ri±,Xi according to ( p.l6|) i.e., Xi = 
uj-ki/oj-p and Ri = ki — Xip. 



3.1.3 Normalization 

The state vector is normalized by eq.( p.26|) : 

^{p',X'\p,X)^ = 2po 6^'\p-p') 6''' 



(3.20) 



The Fock components are normalized so as to provide the condition ( p.2(]| ). Substituting 
the state vector ( ^TTl) in the left-hand side of eq. (|3.2(]| ), we get: 



,{p',X'\p,X)^ = 2po S^'\p-p') N^' 



A'A 



(3.21) 
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with 

N^'^ = J2 N^'^ = 6^'^ , (3.22) 



n 

tX'X 



where is the contribution to the normahzation integral from the n-body Fock com- 
ponent. We represent it in terms of all three sets of variables introduced above: 



n j3u. 



„ n n n 

(2-)^ E ^?::...„.„^?....„.„^^^HE^i)2(E^jn 



(Ti---cr„ i=i i=i i=i (27r)^2£:g^ 

„ n n n j2 

(2-)^ E ^it..„.„^it...„.„^^^nE^-)^(E-^-i)2n^ 



(3.23) 



For the state with zero total angular momentum the normalization condition has the 
form: 

E^n = l- (3-24) 



n 



In this case, the two-body contribution to the normalization integral reads, for constituents 
of equal masses: 

^2 = 7^ /^'(A:i,fc2,P,a;r)5W(A;i + A;2-p-^r)^— ^2(^-p)cir 

i)\k,n) = —— i,\R^,x)—- -, (3.25) 



(27r)3y ^ ^ ' ' Sk {2TifJ ^ ' '2x{\-x 

where = ^jj^^j2o-2(^i' ^2,^5 ^t)| j,a=o and we imply in (|3.25| ) summation over the spin 
indices of the constituents. To obtain eqs. p.25| ) we used the fact that the first integral is 
the two-body phase volume: 



/■/ ^c-l'/l^/, , ■.d^kid^k2^, - , /", ^kdflk^, ^ , 

/ (■ ■ ■)5^^\k^ + k2-p- UT)-^-^2{uj-p)dr = / (• • .)—±2{uj-p)dT , 



and then use the equality -|- 2{u-p)T = {ki + ^2)^ = 4£^, which gives: 2{uj-p)dT = 8kdk. 
The change of variables in ( |3.23| ) is made similarly. For the state with J = the integral 
( p.25| ) and any Nn do not depend on u. 

For the wave function of a system with total angular momentum J = 1/2, the 
integral ( |3.23D does not depend on u as well, since it is impossible to construct any u- 
dependent terms. At first glance, one could construct the following uj dependent structure: 
u' {p)Cju{p) / uj-p {u' and u may correspond to different spin projections). However, since 



u'{,p)li,u{p) = ■^u'{p)u{p) , (3.26) 
m 
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we get: u^' uu^ / uj-p = u^' /m = 25^^\ i.e., this structure does not depend on uo. Hence, 
any A^^ ^ is proportional to 5^^. For the sum of them we get ^ = A6^ ^, and the wave 
function for a system with J = 1/2 is normahzed by the condition A = 1. As an example, 
we give for J = 1/2 the normalization integral for a ra-body system : 

Nn = {2nf - E m:-<rA'S^'\h + --- + K-p-u;T)l[ ' 2(a;-p)dr , (3.27) 

^ o-,o-i---o-„ i=l l^^J ^^fci 

which equals to 1 in the case where the system consists of n particles only. 

We emphasize that the condition (|3.2(]| ), independent of the light-front plane, is 
a dynamical property of the state vector provided by the Poincare group. When J is 
different from 0, 1/2, the integral depends in general on u, though this dependence 
disappears in the sum (|3.22|) calculated with all Fock components. For an exact state 
vector one should still write: A^^ ^ = A6^ ^, and the normalization condition is reduced 
to A = 1, like in the spin 1/2 case. However, for spins higher than 1/2 the tu-dependence 
of A^'^ ^ does not disappear automatically if the state vector is approximated by a limited 
number of Fock components i.e. if a finite sum is retained in eq.( p.22| ). 

Consider, for example, a system with J = 1 (the deuteron for instance). Extracting 
from the Fock components the polarization vector e^^\p) of spin 1 system, we represent 
A"^'^ in the form: 

N::'' = ef\p)Ii:'^ei'\p). (3.28) 
The general structure of the tensor //^'^ is the following: 

/r = -Ang^'' + B^^p^p^ + a^UJ^ + P'i^n + ^n(p^^^ - P'UJ^) + K + ^) , 

(3.29) 

where An, Bn,Cn, Dn, En are constants, M is the mass of the composite system. For 
example, eq. ( |3.28|) has the form in the rest system: 



iV^^ = An6''' + 1^ (n^'n' - ^5^'^ ) , (3.30) 



where n'^ is a unit vector in the direction of uj. After integration of eq. ( |3.30| ) over 



according to (|2.38|) , the irreducible structure {n^' — 5^'^ /?>) gives zero, and we get the 



relation An = Nn- Using eq.( |3.24| ), we get J2An = J2 Nn = 1. As it was mentioned above 



in sect. p.l.3| , the same condition has to be valid without integration over the directions of 
UJ. This means that the exact Poincare covariant state vector should provide zero values 
for the sum of the constants C„,-D„,i?„ (J^Cn = J2Dn = J2En = 0), although each 
of them, or a partial sum of them may not be. While it is not required to do so, the 
cij-dependent structure in eq. (|3.30|) disappears after averaging over n (or u). 
The general normalization condition A = 1 thus writes, with (|3.29| ): 

A = ISy.N^'^ = -I^^l (^,. - ^) = 1' (3.31) 
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where 



(3.32) 



We emphasize that the normahzation ( |3.31| ), according to ( p.22| ) and ( |3.23| ), contains the 
sum over the n-body Fock components for all n. It cannot be generally fulfilled for a two- 
body component, if the contribution of other components is not negligible. The explicit 



expression for A2 in terms of the deuteron wave function is found in sect. 5.1.2 



3.1.4 New representation 



One can see from ( |3.7| ) that the relativistic wave function, in contrast to the non-relativistic 
one, is transformed in each index by different rotation matrices. It is therefore convenient 
to use a representation in which the wave function is transformed in each index by one 
and the same rotation operator R{g, V), rotating, according to (|3.14|) , the variables k and 
n. We define the wave function in this new representation as follows: 



JlO'lj20-2 

02) rD/^-l/"n^ u M^^A' 



X Z}^^;j,{i?(L-i(P),fc2)}<".;,,..(fci,fc2,p,a;r) , (3.33) 

where, e.g., R{L~^(V),p) is given by ( p.5|) with g = L^^iV). 

Under the transformation g, the operator R {L~^{V), ki) is factorized as follows 

R(L-\gV),gk,) = R{g,V)R{L-\V),k,)R~\g,k,) . (3.34) 

Using also the property D{RiR2} = D{Ri}D{R2}, we find that the wave function (|3.33|) 
is indeed transformed like eq.( p.7| ), in which the arguments of all the D-functions are 
replaced by: 

R{g,V)=L-\gV)gL{V) . (3.35) 



We thus have: 



^^^,^igk,,gk2,gp,guT) = {R{g,V)}D^^f^>{R{g,V)} 



X Di;f^,{R{g,V)}^^:„,{k,,k2,p,u;T) . (3.36) 

The equation ( |3.36| ) together with ( |3.14| ) shows that in this new representation and in 
the variables k, n the relativistic wave function transforms exactly as a non-relativistic 
wave function under a rotation R. This strongly simplifies the spin structure of the 
relativistic wave function, making it as close as possible to the non-relativistic one. The 
only difference is the dependence of the wave function on the extra variable n. We shall 
illustrate this dependence in the case of the Wick-Cutkosky model in this chapter, and in 
the case of the deuteron in chapter ^ 
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(-) 

In the particular case of spin 1/2, the matrix -D^.^/ coming in the transformation 
( ^.33| ), takes the form P5| : 



Dl{R{L-\V),h)} 



{kio + m){Vo + VV^) - a-V a-ki 



(3.37) 



[2{k^o + m){Vo + VV^){k^,V, - kyV + mv^)]i/2 
The transformation to the representation (|3.33| ) is similar to the Melosh transformation 



46| . We will come back to this point in sect. ^3] 



3.2 Equation for the wave function 

The equation for the wave function is obtained from the equation for the vertex part 
shown graphically in fig. |TT[ 



>- 






CO T 












Fig. 11. Equation for the two-body wave function. 

It is the analogue, for a bound state, of the Lippmann-Schwinger equation for a 
scattering state. Let us first explain its derivation for the case of spinless particles. In 



accordance with the rules given in sect. p.2.2|, we associate with the diagram of fig. |Tl] the 
following analytical expression: 



r(fci,fc2,p,cjr) = j T{k[,k'^,p,uT')e{u-k[)6{kf -m'^)e{uj-k2)S{k'^ -m^ 

x5(^) {k[ + k'^-p- ujT')d%IC{k[, k'^, ut'; ki, k2, ujt) " . "7^ 

r — le [zn)'^ 



(3.38) 



The delta-function 5^^\- ■ ■)d'^k[ is kept for convenience. The kernel /C is an irreducible 
block which is calculated directly by the graph technique once the underlying dynamics 
is known. We should then express the vertex F through the two-body wave function. 
This can be done by comparing, for example, two ways of calculating the amplitude for 
the breakup of a bound state by some perturbation: 1) by means of the graph technique 
(the result contains F); 2) by calculating the matrix element of the perturbation operator 
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between the bound state and the free states of n particles (the resuh contains We 
thus get:0 

^2, p, ut) = — , 3.39 

s — 

where s = {ki + ^2)^ = {p + ur)'^. The corresponding relation for the n-body case has the 
same form. In any practical calculation of the amplitude, we associate F with the vertex 
shown in fig. ^ and then express T in terms of by eg. ( 3. 39 ). Taking into account that 
s — = 2{uj-p)t, and keeping all the spin indices, the equation for the wave function 
form 

2 



{ki + k2f - M' 



(Tl(72 



xV<:HK, ^r; A;i, k^, UT)5^'\k[ + k', - p - 



LOT 



d k-y d ki^ 
2s h' 2e]c' 



2{uj-p)dT' 



(3.40) 



where we have defined V = —JC/Am?. This equation can be extended in a similar way for 
the many-body component of the state vector. 

For many practical applications, it may be usefuU to express this equation in terms 
of the variables k and n, i.e. express the wave function in terms of 5^^-^0-2(^15 k2,p,uJT) = 
^^^^^{k,n) . With the transformation ( p.33| ), one thus gets: 



A{k'- 



m 



(3.41) 



r^^^,ik',n)V<:Hk',k,n,M 



d'^k' 
) 



The integration variable k' in eq. (|3.41| ) is defined analogously to k in eq.( p.ll|) : 

k' = L-\V')k[, (3.42) 

where V = k[ + k'2. The equation (p. 41 ) can also be obtained from ( 3.4CI|) simply by 
transcribing it in the system of reference where ki + k2 = 0- In this system, the rota- 
tion operators in (|3.33| ) are the unit operators, the wave function ^^lO-^ turns into ^^^a^- 
However, \l/'^„, in the integrand differs from „i . The wave function \I/'{k'n) obtains 
the same form as \P{k,n) (however, expressed through the variable k') in the system 
where k[ + k!^ = 0, but not in the system where fci + A;2 = 0. This fact is marked by a 
"prime". The function is obtained from (p-^, ^, {k[, k'^,p,uT') by simply expressing 

the arguments k[, k^.p.ujr' through /c, /c', n. The relation between these variables has the 
form 



k\ = k' + 



2e, 



-n + 



{Sk' — SkY 
'2£k£k' 



{n-k')n, k[ 



10 



2 I 2 p2 _ 2 ^ 

^k' + ^fc _^ £fc^ —n-k' 



2e, 



"^^k^k' 



(3.43) 



^The coefficient in the relation ( p. 39 ) differs from previous references (see [^) due to the different 
normalization coefficient in (3.1) and the different normalization of the amplitudes given by the rules of 
the graph technique. 
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k2' = -k' + 



2e, 



-n 



(gfc' ~ ^kY 



{n-k')n, k'2 



P 



-LUT 



LUT 
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M2 



4ei 



-n, 



44 



44 - 
4efc 
-M2 



-n, 



20 



Po 



4el 



^k' + ^fc _ 

^k ' 
-M2^ 



(3.44) 
(3.45) 



4£i 



cur 



4£t 



-n, 



4£g, - 

4£i 



(3.46) 

Since the formulae ( P .431 - ^^461) are the Lorentz transformations between the systems with 
^1 + ^2 = and k[ + k2 = 0, the function \P^}^,^{k n) can be also obtained from \P^,^^,^{k n) 
by a rotation of spins by means of the D-functions: 



^/A _ /}a^)t{i^(/^-i(p')^ k'^)}^^D'^^^^{R{L-'{V'), k[)}. 



(3.47) 



For spin 1/2 particles, the matrix Z}''^^{i?(L^^(P), A;)} is given explicitly in eq.( p.37|) . 

In the simple case of a scalar particle, the equation for the wave-function in terms 
of the variables k, n has the following form: 



(4(fc2 + 



m^) - M^)ilj{k,n) 



m 



27r3 



^lj{k\n)Vik\k, n, 



d^k' 
I 

^k' 



(3.48) 



An equation of such a type was also considered in refs. ^ ^ jS^, 

In the non-relativistic limit, equation ( p. 481 ) turns into the Schrodinger equation in 

momentum space, the kernel V being the non-relativistic potential in momentum space, 

and the wave function no longer depends on n. 

We emphasize that the wave function, which is an equal-time wave function on 

the light front, turns into the ordinary wave function in the non-relativistic limit where 

c — > 00. This reflects the fact that in the non-relativistic limit two simultaneous events 

in one frame are simultaneous in all other frames. 

In the variables R±_ and x, eq.( 3.48|) can be rewritten in the form: 



x(l 



X] 



m 



,x] 



27r3 



iIj{R'^,x')V{R'^,x';R±,x, 



d^R\ dx' 



2x'(l 



X 



(3.49) 
.48|) compared 



In this form, this equation is nothing else than the Weinberg equation | 
The advantages of the equation for the wave function in the form 
with (|3.49|) are its similarity to the non-relativistic Schrodinger equation in momentum 
space, and its simplicity in the case of particles with spin. These properties make eq. (|3.48|) 
very convenient for practical calculations. 

The kernel of eq.( p.48|) depends on the vector variable n. We shall see that this 



dependence, especially the part which depends on M^, as given by eqs.( |3.45| ), ( |3.46| ) for 
instance, is associated with the retardation of the interaction. From this point of view, 
the dependence of the wave function ip{k, n) on n is a consequence of retardation. 

In absence of explicit calculations, a very useful information on the global contribu- 
tion of many-body Fock components can be obtained by considering the dependence of 
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the interaction on the total mass M of the system. The current normahzation procedure 
is then to insert in the norm operator the derivative of the potential with respect to M^. 
In the case of a J = state made of two spinless particles, the normalization is defined 
by: 



(2vr) 



ek6{k — k' 



4m2 dV{k',k,n,M'^ 



(27r) 



9M2 



■ipi^k, n) 



(3.50) 



where the second term accounts for the many-body contribution to the norm, J2n>2 ^n- 
The extension to states with J 7^ is straightforward. There are many justifications for 
the introduction of the derivative of the potential V in eq. (|3.5CI| ). A general derivation can 
be done similarly to the normalization condition of the BS function |5^. Apart from the 
fact that it appears in the description of a system of two coupled components in which 
one component is explicitly retained, it is suggested by the requirement that two-state 
vectors describing two infinitesimally close states should be orthogonal. Examination of 
eq.( p.48| ) then shows that such a property is fulfilled by introducing in the scalar product 
of the two state vectors the derivative of V with respect to as it is done in ( p.50 ). The 
extension to the normalization immediately follows from the fact that the orthogonality 
and the renormalization are defined from the same operator, namely the time component 
of a conserved current. 

The consistency of (|3.5CI|) with the general normalization condition ( |3.22| ) supposes 
that the potential V is determined to all orders in the coupling constant entering the 
interaction. In practice, this is limited to the exchange of one or two bosons, which may 
be already enough for the applications considered in the following. Let us just mention 
here that the examination of a few examples (Wick-Cutkosky model, model-perturbative 
calculation) shows that the n dependence of the integrand in ( p.50|) due to the derivative 
of the potential tends to cancel that of the first term. This is in accordance with the 



general discussion given above in section 3.1.3 



When the state is described in relativistic quantum mechanics with a fixed number 
(two) of particles with dV/dM'^ = 0, the normalization for J = is reduced to N2 = 1. 



3.3 Relation with the Bethe-Salpeter function 

To find the relation between the light-front wave function and the Bethe-Salpeter function 
we should start from the integral that restricts the variation of the arguments of the Bethe- 
Salpeter function to the light-front plane: 

/ = jd'x, d'x, S(..x,) S(..x,) t(..,.„p)exp(*.., + , (3.51) 

where ki,k2 are the on-shell momenta: k\ = k^ = m?, and $(xi,X2,p) is the Bethe- 
Salpeter function 

<l>(xi,X2,p) = (0 \T{^{x^)^{x2))\p) . (3.52) 
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Representing the 6 -functions in ( |3.51| ) by the integral form 

6{uj-Xi) = — [ exp(— zci;-a;iai)(iQ;i, 6{uj-X2) = — [ exp{—iuj-X2a2)da2 , 
introducing the Fourier transform of the Bethe-Salpeter function ^{k,p), 

$(xi,X2,p) = (27r)"^/^exp [-ip-{xi + X2)/2] ^{x,p) , x = Xi - X2 , 

<l>(/,p) = J ^{x,p) exp{il-x)d'^x , (3.53) 

where I = {li — h)/'^-, p = li + h, h and I2 are off-mass shell four-vectors, and making the 
change of variables ai + 02 = t, (02 — ai)/2 = /5 we obtain: 

/+00 r+00 
6'-^\ki + k2-p~uJT)dT / <l>{h = ki-ujT/2+uP,l2 = k2-u}T/2-ujl3)df5 . 
-00 J—oo 

(3.54) 

On the other hand, the integral (|3.51|) can be expressed in terms of the two-body 
light-front wave function. We assume that the light-front plane is the limit of a space-like 
plane, therefore the operators ^p{xi) and <f{x2) commute, and, hence, the symbol of the 
T product in (|3.52| ) can be omitted. In the considered representation, the Heisenberg 
operators (p{x) in ( p.52| ) are identical on the light front u-x = to the Schrodinger oper- 
ators (just as in the ordinary formulation of field theory the Heisenberg and Schrodinger 
operators are identical for t = 0). The Schrodinger operator ip{x) (for the spinless case 
for simplicity), which for u-x = is the free field operator, is given by ( p.51| ): 

1 /■ r ^ ^ 1 d^k 

(p{x) = ^^.^^2 I ctik) exp{—ik-x) + a\k) exp{ik-x) —^==. (3.55) 



{2n\ 



We represent the state vector \p) = (j){p) in ( |3.52| ) in the form of the expansion ( |3.1| ). 
Since the vacuum state on the light front is always "bare" , the creation operator, applied 
to the vacuum state (0| gives zero, and in the operators ^{x) the part containing the 
annihilation operators only survives. This cuts out the two-body Fock component in the 
state vector. We thus obtain: 



if'^rv'Tl ^(^1' p. ^r)5^'\ki + k2-p-uT) dr 
2[uj-ki)[uj-k2) J -00 



(3.56) 



Comparing (|3.54|) and ( p.56| ), we find: 



ilj{k,,k2,p,uT)= / ^k + (3u;,p)d(3 , (3.57) 



7t{uJ-p) 

where the argument p in (|3.57| ) is expressed through the on-shell momenta fci, /c2 as p = 
ki + k2 — ujT, in contrast to off-mass shell relation p = li + I2. For the Bethe-Salpeter 
function $(/i,/2) parametrized in terms of the momenta /i,/2 the formula ( p.57|) reads: 

ij{kuk2,P,UJT) = (^•^l)(^"^2) /"*'°°$(/^ = ki-UJT/2 + Uj(3j2= k2-UJT/2-Uj(3)d(3 (3.58) 

Tr[LJ-p) J-00 
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In ordinary LFD, eqs.( |3.57| ) and ( |3.58| ) correspond to the integration over dk_. This 
equation makes the hnk between the Bethe-Salpeter function $ and the wave function i(j 
defined on the hght front specified by u. It should be noticed however that eq.( p.58[ ) is 
not necessarily an exact solution of eq. (|3.48| ), since, as a rule, different approximations are 
made for the Bethe-Salpeter kernel and for the light-front one. In the ladder approxima- 
tion, for example, the Bethe-Salpeter amplitude contains the box diagram, including the 
time-ordered diagram with two exchanged particles in the intermediate state, as indicated 
in fig. ^ in appendix |B.3| . This contribution is absent in the light-front ladder kernel. 

The quasipotential type equations for the light-front wave function derived by re- 
stricting arguments of the Bethe-Salpeter amplitude to the light-front plane z + t = and 
corresponding electromagnetic form factors were studied in refs. |]55| , |56| . 



3.4 Application to the Wick-Cutkosky model 



3.4.1 Solution in the covariant formulation of LFD 

As a simple, but nevertheless instructive, example, we shall derive in this section the light- 
front wave function of a system consisting of two scalar particles with mass m interacting 
through the exchange of a massless scalar particle calculated in the ladder approximation. 
This is the so-called Wick-Cutkosky model. The diagrams that determine the kernel are 
shown in fig. ^ The kernel is given by eq.( p.67|) with /i = 0. Going over from the 
kernel JC to the potential V = — /C/(4m^), introducing the constant a = g"^ / {IQirm?) , and 
expressing ( |2.67|) by means of the relations (|3.43| - |3.46|) in terms of k, k', n, we obtain ||57| : 

V = -Aira/K^, (3.59) 



where 



ik' -kf - {n-k'){n-k) 



n-k ' n-k 

^k' ^k 



(3.60) 



For k, k' <C m, eq.( |3.59| ) turns into the Coulomb potential in momentum space 



Vik',k) 



{k' -kf 



(3.61) 



For a -C 1, \eh\ = \M — 2m\ = ma'^/A <^ m, the wave function is concentrated in the 
non-relativistic region of momenta. The non-relativistic wave function of the ground state 
in the Coulomb potential has the form: 



'TxniK 



5/2 



{k^ + K 



2\2 



(3.62) 
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where k = ym\eb\ = ma/ 2 (it is normalized, however, according to ( p.25| ) with Sk ~ m). 
The integral over d^k' in ( |3.48| ) is concentrated in the region k' n. Therefore, aX k^ n 
the momentum k ' in V(k, A;, n, M^) can be ignored, and from ( p.48|) we find: 



^ ^ _ mViO,k n,M^) r . 

(27r)3(A;2 + /t2) 7 "^^ ^ 

Note, that in ( p.63| ) / ip{k')d?k' = {27T)^ilj{r = 0). The wave function ( |3.62| ) is normalized 
to 1. Substituting in the r.h.s. of eq.( p.63|) the expressions (|3.59| , p.60| ) for V and (|3.62|) 



for ip, we obtain 

Hk, n) = TfX ■ (3-64) 



£k j 

This relativistic wave function of the ground state with zero total angular momentum is a 
good approximation of a more exact one in the range k > k. Corrections of order a log(a) 
should be considered in the range k < k (see [|8|). Although the kernel ( p.59| ), ( |3.60| ) 
contains the modulus \n-k'/ek' — n-k/ek\, one can show that the solution has no "cusp" 
at n-k = 0. This cusp in ( |3.64|) appears due to our approximations. 

One can check on this simple example that it is the retardation of the interaction that 
is the dynamical reason for the dependence of the wave function on the variable n. The 
non-relativistic Coulomb expression for the kernel ( p.61| ) does not contain retardation 



and does not depend on n while the relativistic kernel ( |3.59| ) contains retardation and 
depends on n. This leads to the dependence of the wave function on the argument n. 
It may seem at first sight that the dependence of the wave function on n is not due to 
the contribution of the many-body sectors. Indeed, the coupling of the two-body sector 
with the other sectors must contain a coupling constant, whereas the parameter that 
determines the dependence of the wave function on n is the nucleon mass and does not 
contain the coupling constant. This argument should be taken with some care however. It 
can be seen from eq. (|3.48| ) that the operator fc^ + — M^/4 acting on the wave function 
is equivalent to the potential acting on this same wave function. It is thus proportional 
to the coupling constant. Such a result underlies the disappearance of the n dependence 
of the integrand of eq. (|3.5CI|) previously mentioned. 

We emphasize that retardation leads to both the fi-dependence and the presence of 
the carriers of the interaction in the intermediate state, which contribute to the many 
body sectors. However, these two effects, being important in full measure in a truly 
relativistic system, can manifest themselves in a different way in weakly bound systems. 
Neglecting the many-body sectors does not necessarily entails to neglect the n-dependence 
of the wave function at /c ~ m. The above mentioned Q;log(a) correction to eq.( |3.64| ) for 



instance, which originates from the last n-dependent term in ( p.60|) , is n independent 



and has no counterpart from the many-body sector. As we shall see also in chapter ^, it 
is necessary to take into account the n-dependence of the wave function even when one 
restricts to the two-nucleon sector. 
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The wave function of the 2p state can be found analogously. In the representation 
(PI), it has the form [0: 



1 



(3.65) 



X A;Fi,(^/A;) + Fi,(n) 



(2gfc - Mf 



[n-k) 



2m 



e{-n-k) - 9 (n-k) 



The wave function corresponding to the angular momentum / = 1 contains the spherical 
function Yix{n). This is an illustration of the fact that the vector n participates in 
the construction of the total angular momentum on the same ground as the relative 
momentum k. The dynamical difference between the solution with k\\n and ± n is 
obviously related to the property that some of the components of the angular momentum 
J, before using the angular condition, depend on the interaction. 



3.4.2 Solution in the Bethe-Salpeter approach 

The exact expression for the Bethe-Salpeter function in the Wick-Cutkosky model is found 
in the form of the integral representation [^, ^ and, for zero angular momentum, reads: 

i r+i g{z,M)dz 
= "Ti^i-i {m^ - My A - - zp.l - .6)3 • (^-^^^ 



Substituting ( p.66[ ) in (|3.57|) , and making the change of variable /5 — > P/uj-p, we find: 

_ _Mi_ii) iB r . (3.67) 

2W>/' J-oo J-1 p2 + K2)(l -2J„) +,3(20 — z) — ieY 

where = 1 — 2x, = w? — 

For fixed z ^ zq, the integrand, as a function of /5, has one pole of third order, and 
therefore the integral is zero. For z = zq, the integrand does not depend on (5 and the 
integrand over /3 diverges. To find the value of the integrand ( ^.67] ), we first calculate the 
integral over dz from Zq — e io Zq + e ai t ^ retaining only the term proportional to e, 
and then over d(3. As a result, one obtains [|60i : 



o(l-2x,M) 

i) = . 3.68 

25y^x(l -x) (A; 2 + ^2)2 



The spectral function g(z,M) is determined by a differential equation [^, Q and has 
no singularity at 2; = 0. The exact wave function (|3.68|) is therefore analytic at a; = 1/2. 
The approximate explicit solution found in |59| for g{x, M) has the form: 

g{z, M) = 'l^Tx^K'l^{\ -\A) ■ (3-69) 
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Inserting ( p.69| ) in ( |3.66| ) and integrating over one can recover the approximate solution 
of the Bethe-Salpeter equation given by eg . (17.21) . The discontinuity of the spectral function 
g{z, M) at z = is again the result of our approximation, since the solution ( |3.69| ) 
corresponds to an asymptotically small binding energy. Substituting ( |3.69|) in ( p. 681) , we 
reproduce the expression ( |3.64D for the wave function. The discontinuity of ( |3.69D at z = 
results in the "cusp" of (|3.64|) at n-k = 0. 



3.5 Angular momentum and angular condition 



The wave function (|3.65|) is consistent with the general structure of the light-front wave 
function of a system with total angular momentum equal to 1 (for spinless constituents): 



'ilj\k,n) = f,Yix{k/k) + f2Yix{n) 



(3.70) 



with /i2 = fi,2{k'^,n-k). The states with higher angular momentum are constructed 
similarly, using Clebsch-Gordan coefficients. In this case, it is clear that the angular 
momentum operator in the representation (|3.33| ) has the form: 



J = —i[k X d/dk] — i[n x d/dn] . 



(3.71) 



The same expression for J was found in refs. It can be shown to be quite general 

for spinless particles. Acting on the wave function ( p.64| ) for zero angular momentum and 
on any function depending on the scalars k ^ and fl-k, the operator ( p.71|) gives zero. The 
generalization for the case of constituents with non-zero spins is evident: 



J 



-i[k X d / dk] — i[n X d/dn] + si + s*2 



(3.72) 



where Si_2 are the spin operators of the constituents. 

As we explained in sect. |3.1.4| , the relative momentum k and the spin operators 
Si,S2 have identical transformation properties. In the standard approach, the angular 
momentum operator obtains the same form after performing the Melosh transformation 
p6| , |29[| , once the n-dependent term is neglected in (|3.72|) . In addition, the factor n-k in 
fi,2{k'^,n-k) turns into fc^. At first glance, the dependence of the wave function on kz 



can be interpreted as a violation of rotational invariance (see ref. I^TJ). Indeed, this wave 
function is not eigenfunction of the operator L^, where L is the usual angular momentum 
operator: 

L = -t[kxd/dk] . (3.73) 

The rotations changing the position of the light-front plane t + z = are dynamical, 
the corresponding angular momentum operator contains the interaction, and does not 
therefore coincide with eq. (|3.73|) . 

In the covariant formulation, the action of the dynamical angular momentum op- 
erator has to coincide with the action of the operator ( 3.71 ). This is provided by the so 
called angular condition ( |2.17| ) derived in section ^171. 
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As explained in section the solutions of eq.( |3.48| ) with non-zero angular momen- 
tum are always degenerate if one does not apply the angular condition. This condition 
eliminates this unphysical degeneracy. This can be seen easily in the example of the wave 
function (|3.7CI| ). 

For this aim, it is convenient to represent ( |3.70D in the form: 



'ijj(k,n) = fik + /an . 
By analogy with eq.( |2.39| ), one can construct the operator: 

A = {n-jf. 



(3.74) 



(3.75) 



Due to eq. ( |3.71D , this operator can be represented in the form: A = {n-L)"^. Any n- 
dependence of the kernel does not violate conservation of the projection of the operator 
L on the direction n. Hence, the operator A is conserved. As a scalar, it commutes with 
the angular momentum operator J: [J, A] = 0. Therefore, any solution of eq.( |3.48| ) with 
definite angular momentum is characterized also by the eigenvalues of the operator A: 



Alpa = a Ipa ■ 

For J = 1, there are two eigenstates with a = and a = 1 respectively: 

^po = n{n-k)hQ{k'^,n-k) , 
^pl = k — n{n-k) hi{k'^,n-k) , 



(3.76) 



(3.77) 
(3.78) 



where /io,i are arbitrary scalar functions. In the case where /io,i do not depend on n-k, the 
states ipOA are eigenstates of the usual angular momentum operator squared: L'^ipQ^i = 
1(1 -|- l)ipo,i- Rewritten in terms of the spherical functions, they still have no definite 
projection of the usual angular momentum on the 2;-axis. Thus, ipQ is proportional to 
J2m o-mYimik I k) (with ttm = ^iml^)) and, hence, is a superposition of states with different 
projections m. In the non-relativistic limit = hi, and the degeneracy of the states 
( p.771 ), (|3.78|) is nothing but the usual degeneracy relative to projections of the angular 
momentum. 

In the practical case of the approximate two-body equation ( p.48| ), the states ( |3.77|) 
and ( |3.78| ) are not degenerate but correspond to two different energies. Both of them 



are nonphysical. It is difficult to project explicitly the general angular condition ( |2.18 
on the two-body sector. However, the model solution ( p.65| ), found from the known 
Bethe-Salpeter function, does not contain any ambiguity, and hence, satisfies the angular 
condition. It has the form of eq. (|3.74|) , i.e., is a combination of two states with a = 
and a = 1 (neglecting at this level of approximation the energy splitting). This example 
illustrates the result, proved in section [211] , that the angular condition allows to construct 
the physical solution from unphysical degenerate solutions. Another example is given in 
3^ . The angular condition does not eliminate the n-dependence of the wave function. 



The angular condition in the ladder approximation was found in 
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Chapter 4 

The nucleon-nucleon potential 



The understanding of the A^A^ interaction should ultimately rely on QCD. However, the 
non-perturbative character of this theory in the low energy range makes this goal quite 
far away. For many purposes, using effective degrees of freedom such as mesons in the 
t-channel or nucleons and baryon resonances in the s-channel in order to describe this 
interaction will still be relevant, provided it is used in a given domain of energy and mo- 
mentum. In the present context, this approach is natural because i) the meson exchange 
picture provides a reasonably good description of the NN interaction and ii) it allows one 
to consider specific relativistic corrections. The effects due to a deeper understanding of 
hadronic physics will have to be considered, involving in the simplest case off-shell effects, 
the effective character of some meson exchanges or of hadronic form factors. 



4.1 Mesonic degrees of freedom in nuclei 

The relevance of mesonic degrees of freedom in describing the NN interaction stems 
from the study of peripheral partial waves phase shifts at low energy. These ones are 
dominated by the exchange of the pion, due to both its low mass and the low energy 
domain. In related area, dealing for instance with electromagnetic properties of light 
nuclei, the relevance of the pion degree of freedom relies on meson exchange current 
contributions to the deuteron electrodisintegration near threshold where it allows to 
achieve agreement with experiment at momentum transfers around ~ 0.5 (Gev/c)^ 
|63| . As this contribution is not required by a pure phenomenological description of the 



A^A^ interaction, it is a clear evidence for the role of the pion in the description of nuclear 
forces. In both cases (A^A^ interaction and deuteron electrodisintegration cross-section), 
the role of the pion is related to the peculiar nature of this particle within QCD. The pion 
is the Goldstone boson associated to the spontaneous breaking of chiral symmetry which 
provides its low mass. The 7r-exchange is now a common component of all the models of 
the A^A^ interaction. 

Beyond single pion exchange, one expects contributions due to 2,3,... pion ex- 
changes. They have shorter and shorter ranges, 1 fm or less, and it becomes more and 
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more difficult to identify their contribution to the A^A^ interaction as the number of ex- 
changed pions increases. At this point, models differ from each other. 



Some models |p3, following the philosophy of the Reid Soft Core (RSC) potential 



|65| , simply parametrize this part in terms of Yukawa potentials, but without precise 



relation to meson exchanges. As a support for such an approach, one may argue that 
nucleons have a radius of the order of 0.8 fm (as far as this radius can be identified to 
the proton charge radius) and consequently, below a distance of 1.5 fm, nucleons began 
to overlap, making the description of A^A^ interaction in this range in terms of exchange 
of mesons quite doubtful. 

A second approach relies on the use of dispersion relations to calculate the ir- 
reducible 27r-exchange contribution to the A^A^ scattering amplitudes. The interaction 
potential is obtained by removing the 27r-exchange contribution resulting from the it- 
eration of one TT-exchange which is automatically generated by the dynamical equation 
used to calculate the two-body wave function. This allows one to extend the potential 
to distances as low as 0.8 fm. Notice that the full 27r-exchange scattering amplitude 
has covariance properties that are lost by the separate vr-exchange and 27r-exchange con- 
tribution to the A^A^ potential. The short range part is largely phenomenological, with 
the consequence that the Lorentz structure of the A^A^ amplitude may not be correctly 
determined, quite similarly to the previous approach. An important point to be noticed 
is that this approach relies on physical vrA^ scattering amplitudes with the consequence 
that no form factor is introduced at the nNN {nNA, . . .) vertices. Off-shell effects, which 
have a shorter range, are incorporated in the phenomenological part. On the other hand, 
the irN scattering amplitude satisfies properties expected from chiral symmetry, although 
a pseudoscalar ttNN coupling is used. The large contribution from the excitations of 
NN pairs (Z-diagrams) to the vrA^ scattering amplitude expected in this case is cancelled 
here by a contribution involving mainly the A resonance intermediate state (calculated 
via dispersion relations). Finally, although it does not incorporate explicit p-exchange 
contribution, this approach can account for it, through the 2tt interaction in a P-state in 
the t-channel. 

A third approach is based on a field-theoretical description of the meson-nucleon 
interaction, including vertex form factors 0, In the simplest approximation, the 

A^A^ interaction results from a sum of contributions involving single meson exchanges: 
TT for the long range part, "a" giving attraction at intermediate distances (1 fm) and 
accounting for the exchange of 27t in a S-state, p accounting for the exchange of 27r 
in a P-state, and u for the short-range repulsive part. Other well known mesons such 
as 7],ai,6 may also be considered, but are quantitatively less important. While some 
phenomenology is involved in the vertex form factors or through adjustments of coupling 
constants, this approach has the great interest to provide a parametrization of the A^A^ 
interaction in terms of explicit exchanges of mesons, which are expected to be important, 
if not dominant, in some channels. The fit to A^A^ scattering is quite satisfactory and 
there is no real indication at present that such an approach breaks down, even at short 
distances where it could fail. Because it keeps the full Lorentz structure attached to the 
different exchanged mesons, this approach is quite appropriate to the present purpose. 
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The above approach can be improved by considering two-meson exchanges. This is 
especially important for a better account of the " cr" exchange contribution, whose effective 
character is well known. As in the dispersion relations approach, the iterated one meson 
(pion in practice) exchange has to be removed from the full 27r-exchange contribution to 
the NN amplitude. The above contribution also involves A resonances (and other higher 
resonances) in intermediate states in a quite similar way as Van der Waals forces in atomic 
physics are generated. The non-local character attached to such contributions is likely to 
be accounted for only in an approximate way by single meson exchanges (as "a" and p 
mesons) . 



4.2 The non-relativistic NN potential 

The non-relativistic NN potential is to be used with the Schrodinger equation. In mo- 
mentum space, which is more appropriate for an extension to relativistic calculations, it 
reads: 

-m + / {k\v\k')ij{k') = Em , (4.1) 

— * — * ^ — # — * 

where k represents one half of the relative momentum of the two nucleons, k = ^{ki — ^2), 
and similarly for k'. In r-space, where many models were developed, the corresponding 
equation would be: 

ij{r) + / d^r'{r\V\f')^{f') = Ei){r) . (4.2) 

m J 

Typically, in a meson exchange theory, {k\V\k') takes the following form in the simplest 
case: 

Vik, k') ^ {k\V\k') oc f , (4.3) 

p? + {k- k'Y 

while in r-space, one would get 

V{f,r') = {f^V\f') oc 9^S^^\r-f')^^^^^^ . (4.4) 

where f represents the relative distance between the two nucleons, r = ri — r2. The 
same interaction that is non-local in momentum space, eq. ([4.3|) , is local in r-space. This 
explains why the last one is often preferred in practical calculations. 



4.2.1 Spin structure of the non-relativistic A^A^ potential 

Quite generally, the NN non-relativistic Galilean invariant potential contains 5 terms 
with different spin structure (assuming that strong interactions conserve isospin). It may 
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be written in momentum space as: 

V{k,k') = y^(fc,^') + \/5s(fc,fcVr^2 + 1^5oAfc')^('?i + ^2)-[^><F] 

+VT{k, k') [{k - k'far^2 - 3(A? - k')-ai {k - k')-a2} 

+Vso2{k, k')ar[k x k'] ar[k x k'] , (4.5) 

where Vc{k,k'), Vss{k,k'), Vso{k, k'), Vxik^k') and ^502(^,^0 stand respectively for 
central, spin-spin, spin-orbit, tensor and quadratic spin-orbit forces respectively. They 
are predominantly scalar functions of [k — k 'Y (with further slight dependence on /c ^ and 
k '^) and contain isospin independent and dependent terms: 

V(k, k') = V\k, k') + V^{k, k')frT2 . (4.6) 

Other structures, which have an off-energy shell character, are discarded on the basis that 
their effect can be re-absorbed in the shorter range part of the interaction model, mostly 
phenomenological. Such terms have been considered in ref. ||69|| . 

In r-space, and assuming a dependence of V{k,k') on {k — k'Y only, the above 
interaction takes the following form: 

V{r,r') = 5^'\r- r') \^c{r) + Vss{r)Bra2 + Vso{rf-^^-l 

I -2 ^r(^2 I + Vso2{r) + ■ • • | f> , (4.7) 

where the points indicate the existence of other terms involving V502, but usually ne- 
glected. The potentials Vc{r) and Vss{r) can have a linear fc^ dependent term [70| , 



which, by using an appropriate transformation can be dealt with in coordinate space 
without too much difficulty. In momentum space, such terms can be kept to any order 
without further difficulty. These terms have quite different origins. Some simply come 
from the mathematical structure of Dirac spinors. Other ones come from the dependence 
on the s- variable arising from two- meson exchange {y/s is the energy of the system in 
its cm. system). This s-dependence can be turned into a k^ dependence if one assumes 
s ~ 4(A;^ + m^). This procedure amounts to neglect off-shell effects that can be incorpo- 
rated phenomenologically through a fit to experimental data of the short range part of 
the NN interaction. 

Typically, the effects under consideration are of order k'^ /m? and look like kinemat- 
ical relativistic corrections. Knowing their origin is however not sufficient to fix them in 
the NN potential. The equation, which the NN potential has to be associated with, has 
to be precised in order to have an unambiguous prescription. Thus, for the Schrodinger 
equation (or equivalent ones such as the Lippmann-Schwinger equation for the T-matrix 
for instance), the appropriate definition of V{k^ k') is: 

/rn u'^i {k)Oiu'''^ {k') u"'' {-k)02u''2 {-k' 



Am^ifi^ + {k- k'y) 
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where the spinors are normahzed as uu = 2m, hence the factor 4m^ in the denominator in 
( [4.8|) , and O represents the vertex describing the interaction of nucleons with the mesons 
of interest. The factors ^m/ek and ^m/ek', that are neither 1, as in the non-relativistic 
hmit nor m/sk and m/ey as one may expect from the normahzation factors relative to the 
two nucleons in the initial and final states, are required to satisfy a unitarity condition . 
Furthermore, their omission in earlier calculations of the NN potential was responsible 
for a large energy dependence of the 27r-exchange contribution (obtained by removing the 
iterated one pion exchange from the full 27r-exchange contribution to the NN scattering 
amplitude). The presence of the factors \Jm/ek and \Jm/ek' in (|4.8|) is important in 
comparing relativistic calculations to non-relativistic ones. They already include some 
sizeable relativistic kinematical corrections. 

Anticipating a comparison with the light-front equation for the A^A^ system, as given 
by eq. (|3.48| ) for instance, we can rewrite the Schrodinger equation ( [4.1|) with eq. ([4.8|) for 
the potential and = feg/m, as follows: 

(4.9) 

Quite generally, the on-shell A^A^ amplitude, from which the A^A^ potential is derived, can 
be expressed as a sum of independent invariants. Their number, five, is in relation with the 
number of independent terms in the non-relativistic limit given by eq.( [4.5| ). Their choice 
is not unique however. They may be built from the different tensors 1(5), ^^{PS), 'J^{V), 
7^75(^4) and a^y{T). Other choices, in closer relation to the 27r-exchange contribution 
can also be made ||6^. They differ by off-shell effects, which may be accounted for in a 
phenomenological way by the fit of shorter range contributions to experimental data, as 
already mentioned. 



4.2.2 Physical inputs 

The A^A^ potential due to the exchange of single mesons which we consider here can be cal- 
culated from the following Lagrangian densities describing the meson-nucleon couplings.0 
(z) Pseudoscalar mesons (tt, rj) (PS coupling): 

= i g ^^^^ (j)(P') ; (4.10) 

(m) Scalar mesons (a, 6): 

= g tp^(l)(--') ; (4.11) 

(m) Vector mesons (p, u): 

= ^P [gY<Pi^ + (//4m)a^'^(9^</)(^) - S.^^f))] ^ . (4.12) 

^In the present paper we use the definition of coupUng constants corresponding to g^j^j^/An k, 14, in 
contrast to refs. |7l|, |7^, where the definition with g^p^N ~ 14 has been used. 
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For simplicity, the isospin degrees of freedom have been omitted. 

Other types of couphng may be considered, such as the pseudo-vector one for the ir 
(or 7]) meson for instance: 

(w) Pseudoscalar mesons (vr, rj) (PV couphng): 



2m 



(4.13) 



For a transition involving on-mass-shell nucleons described by positive energy spinors, 
couplings {i) and (iv) are equivalent. They strongly differ off-mass shell as it is well 
known. The first one can give rise to larger effects due to the so-called Z-diagram it leads 
to, while the second one does not. The choice of the coupling is deeply connected with 
the way in which chiral symmetry is chosen to be realized. 

As far as nucleons described by positive energy spinors are retained in the descrip- 
tion of the NN system, the potential is unique. As will be seen later in section |4.3.1| , 
different potentials can be obtained on the light front depending on the choice of the par- 
ticular coupling used at the meson-nucleon vertex (like for instance the PS or PV nNN 
couplings). Using the standard approximation made in the meson propagator consisting 
in taking ko — k'Q = 0, neglecting therefore retardation effects, one gets in the NN c.m.s.: 



V{k,k') 




9-kNN 



u{k)'~^^u{k') u{—k)'~^^u{—k') 



9(tNN 



4m2(/i2 + (A; - k'Y) 
u{k)u{k')]^ [u{-k)u{-k')]^ 



ri-T2 



4m2(/i2 + (fc_fc')2) 



Am'^{lil + {k-k'Y 



xu{—k) 



-u{k) 

I 

/ 



/ 



{9 + f)i,-,^{h + K), 



u{k') 



2m 




(4.14) 



where ki, k2 {k[, k'2) are the initial (final) on-shell nucleon momenta. Only a few significant 
exchanges have been displayed in (|4.14| ). Other ones, t], 6 or p can easily be obtained 
by removing or inserting the appropriate isospin dependence. Hadronic form factors may 
also be introduced by giving to the meson-nucleon coupling constants a dependence on 
the meson four-momentum. 



4.2.3 Choice of the parametrization 

There are various single meson exchange models in the literature. They are not all equiv- 
alent, with the consequence that they should be used within a defined scheme. One of the 
most ambitious model is the field theory motivated Bonn model (Bonn-E) . Its energy 
dependence, which has a relationship with retardation effects included here in the kernel 
of the interaction, has prompted the authors to derive energy independent models (now 
denoted QA and RA) that can be used more easily. The emphasis on reproducing the 
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deuteron D state properties of the original model has biased this derivation. The model is 
unable to account for the observable mixing angle ei beyond 100 MeV. These drawbacks 
have been corrected in later versions (QB, QC, RB) [^] by relaxing the constraints on 
the deuteron D state properties that are known to be model dependent ||73l . 



Other models differ by the values of the coupling constants, including the associ- 
ated form factors. This is expected for short-range contributions which have an effective 
character anyway, but this also occurs for the nNN coupling constant, which may be 
sensitive to isospin breaking effects |6^. The last one is important for a fine description 
of the NN interaction but unimportant in the present review. The difference with Bonn 
QB and QC models , which motivation was shortly explained above, involves changes 
in both the nNN coupling and those determining short range contributions. 

The difference between momentum and coordinate space models has conceptually an 
other origin. Working in one or the other should indeed be equivalent. In practice however, 
this is true as far as the linearly fc^-dependent terms, which are tractable in r-space models, 
represent an accurate description of the full momentum dependence of the A^A^ interaction 
model. Recent developments, dealing with the well known vr-exchange, show that off-shell 
effects involving fourth order terms in the nucleon momentum are relevant for describing 
the NN interaction in coordinate space [B^. Their neglect explains a large part of the 



differences in the bare predictions obtained from the Bonn QB and Paris |7y] models, 
which reproduce the same NN scattering data otherwise. In particular, the deuteron 
D-state probability of coordinate space models is expected to be systematically between 
0.5 % and 1 % higher than in momentum space models. 

Other differences may be mentioned. The genuine term, (fi-[k x k'] (T2'[k x k'], in 
eq. (K:.5|) is approximately replaced by the term ayl a^l + cttI (Tyl in coordinate space 
models. This can introduce some bias in fitting the potential parameters to A^A^ scat- 
tering data. The full 2'7r-exchange contribution considered in the Paris model does not 
seem to bring significant improvement upon a single meson exchange model, probably 
because these ones contain enough parameters (coupling constants together with vertex 
form factors) to accurately account for these A^A^ data. This is indirectly confirmed by 
the Argonne model which also does well with these data although it does not incor- 
porate k"^ dependent terms as in the Bonn R, Paris or Nijmegen models and, on the other 
hand, fit the deuteron quadrupole moment to the observed one without any regards to 
the contribution of meson exchange currents. 

Altogether, it seems that A^A^ scattering data put strong constraints on the A^A^ 
interaction models and within a given scheme do not allow for much freedom. Differences 
are closely related to the scheme under consideration and, most often, involve terms of 
order (fc/m)^, that are of relativistic origin. For some part, the present models may 
be unitary equivalent, implying the existence of corrections pertinent to the model under 
consideration for any calculation of the physical observables. The choice of some particular 
model depends on a prejudice ||69|, models that keep track of these theoretically motivated 
{k/mY terms should deserve a special attention however. 



59 



4.3 Meson exchange interaction on the hght front 



As already mentioned, the approach to the non-relativistic NN interaction is semi-phe- 
nomenological. This means that the meson-nucleon field theory is used to obtain a form 
of the A^A^ interaction kernel (OBEP, as a rule), and the parameters of this kernel are 
found from fitting the experimental data in the framework of the Schrodinger equation. 

A similar approach can be based on the relativistic quantum mechanics in the front 
form. Besides the Hamiltonian, one should introduce the interaction in all the generators 
of the Poincare group, which change the orientation of the light-front plane 0. 

This approach has been developed in a number of studies (see, in particular, [^, 
0, 0, |7|, |7B|, |7|, |7S|, 0). We do not describe them here in details (see [g^, ^ for a 
review) . One should mention that in this approach the full program of fitting the phase 
shifts and the relativistic kernel has not yet been realized. 

The kernel in these references is considered as purely phenomenological, i.e. it is 
not derived from light-front field theory, even in the framework of OBEP. It is assumed 
that it does not depend on the orientation of the light-front plane. However, if one 
uses the field theory to construct the off-energy shell light-front kernel, this one always 
gets a dependence on the light front through the four-vector uj or its spatial part n (see 
examples in sect. |4.3.1| ). It seems reasonable therefore to use this kernel as an input in the 
light-front version of relativistic quantum mechanics. This cannot be done by a simple 
replacement of the ri-independent kernel in the above versions of the relativistic quantum 
mechanics by the n-dependent one. This would destroy the commutation relations, found 
for the n-independent mass operator. In order to restore the covariance, one should find 
a solution of the Poincare group algebra with this new n-dependent mass operator. Due 
to covariance, the exact on-shell amplitudes in this approach will not depend on u (with 
the accuracy of numerical calculations). 

It is quite probable that in both phenomenological light-front approaches (with the 
u - independent and dependent kernels), in spite of the completely different behavior of 
the off-shell amplitudes, the phase shifts and the binding energies can be equally fitted, at 
least in the two-body problem. This implies that there exists an unitary transformation 
between corresponding sets of the Poincare generators. However, this does not imply 
their practical equivalence, since this unitary operator has to be dynamical and highly 
complicated, since it has to change the off-shell behaviour of all the partial amplitudes. 
It hardly can be useful in practice. 

It seems therefore reasonable to develop a phenomenological scheme inspired from 
field theory and to start with the OBEP kernel in the light-front relativistic approach. This 
way seems rather promising, since it reconciles the main features of the field-theoretical 
kernel on the light front with relativistic covariance. It has not yet been developed in its 
full form. 
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4.3.1 OBEP on the light front 

The contribution due to single meson exchange, which is shown in fig. ^, can be obtained 
from the meson-nucleon interaction described by the Lagrangians given by eqs. (CT) - 
( [4.13 ). The analytical expression for the pion exchange in the pseudo-scalar coupling, for 
instance, is a straightforward generalization of eq. (|2.67|) . It has the form: 
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(4.15) 



where F is a phenomenological form factor, assumed to depend on the same argument as 
the denominator of the kernel. Here and below it is introduced "by hand" . 

The 7] exchange contribution can be obtained from ( |4.15 ) by removing the isospin 
factor and replacing glj^j^ by g'^j^^. The contributions due to a- and ^-exchanges are 
obtained from the previous ones by removing the 75 matrices and changing the overall 
sign. As for the vector meson exchange, it contains a peculiarity because of the contact 
terms due to the vector meson propagator and the coupling with derivative. We therefore 
explain this point in more details. The starting expression for the vector meson exchange 
contribution, obtained by applying the rules of the graph technique to the diagram of 
fig. ^, has the form: 
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Ti — te 

where k is the meson momentum. The factors (i) and (— z) at {k — ujti) are associated 
with outgoing and incoming lines, according to sect. ^.2.4| . The factor ujti, everywhere in 
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the difference {k — ujti) in the expression ( |4.16| ), incorporates the contact terms for both 
the vector meson exchange and for the couphng with derivative. 

After integrating over r and introducing the form factor F, we get for the potential 
V= -/C/4m2: 
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(4.17) 



where K"^ is given by eq. ( |3.60D . 

The contribution of the isovector exchange (p-meson) is obtained by multiplying 
( [4.1?1 ) by the factor ti-T2- The terms in ( |4.17 ) proportional to r, r' and rr' were omitted 
in ref . g| . 



For reasons already explained, the pseudovector ttNN coupling is sometimes pre- 
ferred to the pseudoscalar one. To illustrate the treatment of terms corresponding to 
derivative couplings of the meson field, we just give the result here: 
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(4.1^ 



The comparison with the expression obtained with the pseudo-scalar tcNN coupling evi- 
dences differences which depend on the four- vector a;^ and vanish in the limit of on-energy 
shell nucleons (r = r' = 0). 

It is also interesting to compare the interaction on the light front, eq.( |4.15| ), with 
the Bonn-Q model (|4.14] ), which contains some relativity. Differences involve again terms 
depending on the four-vector u^, which vanish on-energy shell. Further differences, due 
to the time component of the four-momentum of the meson in its propagator (retardation 
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effects), vanisli on energy sliell. Tlie difference in tlie normalization factors Jm/sk and 



m/e^, discussed in sect. |4.2.1| , is related to tlie particular equation for the wave function. 



4.3.2 Beyond OBEP on the light front 

The interaction due to single meson exchanges can produce contributions depending on 
the arbitrary four vector uj. While this dependence vanishes for the Born amplitude 
(which is calculated for on-shell particles), there is no guarantee that this property holds 
to all orders. This implies an asymptotic scattering behavior, and therefore a scattering 
amplitude which depends on cj. A quite similar problem arises in other approaches based 
on a description referring to the plane determined by o" = Ax with the particular choice 
A = (1,0,0,0), but here the difficulty is more striking as results may depend on the 
arbitrary orientation of the three- vector uj{n) while physical ones should not. For example, 
the relativistic wave function originating from the vr-exchange, calculated in a perturbative 
way in the first 1/m order, reads |72]: 
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(4.19) 

The quantity is positive in the case of a bound state such as the deuteron, and it is 
negative for a scattering state. 

The defect of the wave function (^4.19|) is that it has a n-dependent residue at the 
pole at A; = in, that can contribute to an observable amplitude. A similar defect in the 
scattering state wave function appears in the dependence of asymptotical states on n, i.e.. 



in the on-energy-shell amplitude, in contradiction with the results of sect. p72 . 

Knowing that the sum of the contributions of all irreducible diagrams to all orders 
should reproduce results obtained from Feynman diagrams, which are independent of u, 
one can guess that higher order terms in the interaction are essential for removing the 
above unpleasant feature. In this section, we explore their contributions and the role they 
may play. As will be shown below, the defect of the wave function ( |4.19| ) is corrected by 



incorporating the contact term in the A^A^ kernel [Bl . 
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Fig. 12. Irreducible diagram with contact term. 
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To give some insight about the relevance of these extra terms in the interaction, 
we will concentrate on a selected one, where the contact term discussed for the spin 1/2 
particle (eg. ( |2.60D ) is operating (see fig. |T2D. This will be done by retaining the relevant 
lowest 1/m order terms, and for contributions in relation with the single vr-exchange. At 
this order, only a- and cj-exchanges contribute, so in fig. O one wavy line corresponds 



to TT-exchange and the other wavy line corresponds to cr and u. The corresponding 
contribution to the kernel has the form: 
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where V'^^^ represents the dominant contribution to the interaction at zeroth order in 
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The insertion of dV^'"^^ in an equation such as ( p. 411 ) provides an extra, n- 
dependent, contribution to the wave function. In a perturbative calculation it is given 
by: 
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where ip^^\k") is a zeroth-order non-relativistic wave function. 

One should also take into account all the iterations of the above contributions with 
the non-relativistic potential V'^'^'^, since these iterations remain in the first 1/m order. 
They are represented as 
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where the integrand incorporates the function dip^'"^ itself. Therefore eq.( [4.22| ) is an 
inhomogeneous equation for . 
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Applying eg. ( p.41| ) at the lowest order considered here allows to replace the integral in 

2 ^ 

the last bracket by —jj^ipQ{k'). It is then easily seen that this term adds to the first one 
in the bracket to cancel the front factor — - at the r.h.s. of ([4.24|) . One thus gets: 
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The pole in eq.( [4.19| ) at k = in {k = ko in continuous spectrum case) has now disappeared. 
This result is important since it guarantees that the observables, such as asymptotic 
normalization or phase shifts, do not depend on the arbitrary direction n, though the 
wave function |4.25| depends on it. 

In r-space eq. ( [4.25|) takes a simpler form: 
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which evidences a radial part in close correspondence with that appearing in the pair 
term contribution to meson exchange currents or also in relativistic components in other 
approaches. 

To calculate the contribution dil)^'"^'^ one has to substitute for 5il)'^^'>'^ + dil)^'"'^'^ 
in the integrand in eq.( [4.22| ) the wave function ( ^.25| ). The result of this substitution 
coincides with —dil)'^'''^'^- Hence, the equation (|4.22| ) turns into 
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This is the homogeneous equation relative to 5il)2'"^'^ + dip^'"'^'^ . However it does not 
determine any eigenvalue, since the energy E m = —mE is already fixed by the 
equation for the complete wave function. Therefore its solution is zero: 



<5V^r'""(^,n)+5V^J''^+-(fc,n) 



0. 



(4.28) 



The significance of this result may be obtained from a comparison with other relativistic 
approaches |T2[. In these approaches, a term like ( |4.26| ) is associated with a pair of 
nucleons with positive and negative energies. Their interaction due to a and u exchanges 
produces a correction proportional to — (instead of + V"^ as in eq.( [4.22|) ). The 
second term in eq.( 4.20| ) has therefore the effect to remove the contribution that should 
be absent in any case. 

In this section, we have considered some of the second order n dependent contribu- 
tions to the irreducible interaction V . This has been done to lowest 1/m order. These 
contributions which are the first ones of an infinite series seem to play an important 
role to make the theory consistent, accounting in particular for the n independence of 
the asymptotic wave function, which are related to physical properties. They also show 
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that solving the equation of motion without using the appropriate interaction may lead 
to wrong results. This consideration makes the relationship with other relativistic ap- 
proaches closer. The present study is, to a large extent, exploratory. It suggests that the 
interaction should fulfill some constraints whose determination may be helpful for further 
work. 
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Chapter 5 

Applications to the nucleon-nucleon 
system 



As a direct application of the formalism developed above, we consider in this chapter the 
two nucleon system. Since its non-relativistic phenomenology has been developed over 
more than twenty years, it allows a sensible discussion of the new features brought about 
by relativity. 



5.1 The deuteron wave function 



5.1.1 Spin structure 

The construction of the deuteron wave function is a nice illustration of how to construct 
a state with non-zero angular momentum in the covariant formulation of LFD. For this 
nucleus, = 1"*" (pseudovector particle), the decomposition of the wave function ^^20-1 
in independent spin structures has the general form 



with 



<P^^^,^{h,k2,p,ujT) = y^e'^{p)u'''{h)<P,U,u^'{h) , (5.1) 



{ki - k2)f, , 1 , {ki- k2)i,uJ 

/m"^ m uj-p Zmuj-p 

+V'57r^ 75W7(^i + ^2)1.(^1 - k2)pU^ + ^ r!^^^ , (5.2) 

and m is the nucleon mass. The light-front deuteron wave function is determined by 
six invariant functions v^i-e, depending on two scalar variables, e.g. s = {ki + k2Y and 
t = {p ~ kiY . The extra spin structures in front of v^a.g are constructed by means of the 
four-vector uj. Other possible structures (e.g., e^M2a)7^t/cWi, etc.) are expressed through 
the six structures given by ( ^.2|) . The wave function ( p.l|) is kinematically transformed 
under rotations and Lorentz transformations g of its arguments according to eg. (|3.7|) . 
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The transformation to the new representation ( p.33| ) has the form: 



X DiJi,{i?(L-i(P),fc2)}<.i(A:i,fc2,P,^r), (5.3) 
where, e.g., R{L~^ (V) , p) is given by (p.5[), in which g = L^^iV), and V is given by ( p.l3| ). 



/ 1 \ 

The matrix D^l, in (|5^ ) is defined by eq. ( p.371 ). 



The transformation law of the wave function ( |5.3| ) is given by eq.( p.36D . In this new 



representation given by (|5.3| ), the rotation operator is the same for all the spin indices 
and coincides with the operator rotating the variables k, n in eqs.( p.l4D . 



In the variables k and n, the problem of constructing the wave function is the same as 
in the non-relativistic case: we have to construct the general form of a pseudovector from 
the variables k,n, the Pauli matrices and two-component spinors. Since ai,a2 = ±1/2, 
the wave function is a 2 x 2-matrix. In this respect, the only difference from the non- 
relativistic case is the presence of an additional vector n that increases the number of 
independent structures. The decomposition of the wave function ( |5.3| ) takes therefore the 
form 

•^.2^1 (k, n) = ^wl^^\k, n)aywl , (5.4) 

with 

Mn) = /i^^ + /2^(^^-a) + /3^(3n(n-c?)-cr) 
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2k 

^ -^^/f ^ ^ •^^"^^^^ X n] X (?] , (5.5) 

where w is the two-component nucleon spinor normalized to w'^w = 1. The relation 
between and ip is the same as the relation between the spherical function Y-^{n) and 

n. The scalar functions /i_6 depend on the scalars k"^ and n-k oi k = \k\ and z = cos{nk). 
For an isospin zero state, the Pauli principle results into: 

/l,2,3,5(^, Z) = /l,2,3,5(^, -Z), fifi{,k, z) = -fi^&{k, -z) . (5.6) 

The reduction of the number of variables can be summarized as follows. The initial wave 
function ( |3.7| ) depends on 4 four-vectors, i.e. on 16 variables. Four mass-shell constraints 
reduce this number to 12. The four- dimensional conservation law ( p.8|) reduces it to 8. 
Lorentz-invariance eliminates the dependence on the velocity v of the system of reference, 
i.e., reduces the number of variables to 5. The wave function ( |5.5|) just depends on the 5 
independent components of two vectors A;, n, (n^ = 1). Because the rotational invariance 
eliminates their dependence on the three Euler angles, the scalar functions fi depend only 
on the two variables [k and z or R\ and x). 



68 



In the region k <^ the functions /s-e which are of relativistic origin be- 
come neghgible, /i 2 do not depend anymore on z and turn into the S- and D-waves: 
fi ^ Us,f2 ~ —ud- From the general decomposition (|5.5|) , one recovers the usual non- 
relativistic wave function: 



1 1 



3k{k-a) 
fc2 



(5.7) 



The representation ( ^.2|) is however convenient in the calculation of physical observables, 
like electromagnetic form factor for instance. The decomposition (|5.5| ) is more adapted 
for comparison with the non-relativistic limit and will be chosen in this chapter to present 
numerical results. The relations between the functions ipi-Q and /i_6 can be found by 
comparing eqs.(|5.1|),(|5.2|) with eqs.(|5.4|),(|5.5|) in the reference system where ki + k2 = 0. 
Indeed, both representations coincide in this system, since all the D-matrices in ( ^.3|) turn 
into unit matrices. These relations are given in appendix 0. 

The number of components of the wave function is a priori not the same in different 
relativistic approaches. For example, the Bethe-Salpeter amplitude of the deuteron is 
determined by eight components |^ 



(see sect. ^.2|) , whereas the Gross wave function [|T3 
contains four components. It does not follow, however, that these components do not 
lead to physical observable consequences. In a given framework, they are unambiguously 
related to measured cross sections. The convenience of one approach or another de- 
pends on the problem under consideration. The relation between the Bethe-Salpeter and 



light-front components of the deuteron functions is given in section 5.2. We emphasize 



that within the LFD the number of components is the same both in the non-covariant 
approach and in the covariant formulation. In the non-covariant approach, these compo- 
nents manifest themselves as six independent matrix elements of the wave function ^'^0-2 
with m = 0, ±1, ai^2 = ±1/2 instead of two matrix elements in the non-relativistic case. 
The representation of the wave function (|5.5|) in terms of the well defined spin structures 
and corresponding components /i_6 is however much more transparent than in terms of 
those matrix elements. 

Note that when a deuteron state at rest, and described from two positive energy 
spinors for its constituents, is seen from a different inertial frame, it automatically acquires 
four extra components depending on the velocity of this frame. This also holds for the ^5*0 
scattering state (one extra component in this case). The equality of the total number of six 
components thus obtained with the number of components in the light-front approach is 
not a coincidence, since it is related in both cases to the appearence of an extra vector. The 
extra components in the light-front approach just represent the dynamical counterpart 
arising from the boost of the system at rest to some velocity v, and the vector n is the 
counterpart of the direction v when f — > c. 
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5.1.2 Two-body contribution to the deuteron normalization 

For the contribution of the two-body Fock component, the normahzation integral ( |3.23|) 
is rewritten as: 

N^'' = / y ^fa. a ^f. , . + h-p- UJT)^^2{uJ-p)dT. (5.8) 

^ (27r)3 J ^2<T2ji<xi jaaajiai \ 1^ 2 y ] ^^^^ ^^^^ \ y) \ ) 

Substituting the wave function ( |0| ) into the integral (p.8D, we get: 
where 



and is given by eq.( 

Substituting in ( p.31| ) the integral given by ( p.9| ) with the wave function (|5.2|), one 
can obtain the normalization condition for the deuteron wave function in terms of the 
components ipi-e- Expressing the functions (pi^Q through fi_Q by the formulae ( |(J.3| ) - 
( p.8|) given in appendix we can represent it in terms of the /i_6 functions. It is however 
simpler to do this directly in terms of /i_6- For this goal, consider the integral ( |5.8|) in 
the reference system where fci + ^2 = 0. The deuteron wave function is given by ( |5.5|) . 
We then get for ^5 the basis where A', X = x,y, z: 

= 7^ / Tri^phk. n)Uk. n)}— . (5.10) 
The general structure of the integral ( p.lO| ) is given by ( p.30| ) and A2 is found from ( p.31| ): 

A2 = ■ (5-11) 

Substituting in (|5.10|) the wave function ( p.5| ), we get: 
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where 



/ = fl + fl + fl + fl{^ + z') + fl{l-z') + fl{l-z') 

-hhH - 3^2) + 4/2/4^ + 4/3/4;^ (5.13) 

and z = cos{qn). In the case where the kernel of the interaction, V, is independent of 
M, one has A2 = 1. Keeping the dominating contribution in the region k <^ m, where 
the functions /a-e are negligible and /i = us, f2 = —ud, we recover the non-relativistic 
normalization condition for the deuteron wave function: 



^ J{ui + ui)edk = i 
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5.1.3 Equation for the wave function 

Keeping all the spin indices, the equation for the wave function '^'^-^^^ form already 

indicated in eq. (|3.40| ). In terms of the variables k and n, this equation can be rewritten 
as: 



4{k'^ + 
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27r3 
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(5.14) 



According to our discussion of section |3^ , the function is obtained from 



^■^,^„,^{k[, k'^,p,LJT') simply by expressing the arguments k[, k'2,p,ujT' through k,k',n, as 
given by eqs.( p.43| - p.46| ). Expressing the function tp through ^^^g.^ with the help of the 



D-matrices, in eq. (|5.3|) , the wave function ijj'{k\n) can then be obtained from ijj{k',n) 
by a rotation of spins according to: 



^' = /^(^)t{i?(L-i(p'), k'2)}^D^'^HRiL"\P'), k[)}, 
where the matrix D^^^{R{L^^(V), k)} is defined by (|3.37| ). 



(5.15) 



5.2 Connection with the Bethe-Salpeter amphtude 

The Bethe-Salpeter function of the deuteron $(p, /ci, /C2) is a 4 x 4 matrix depending on the 
two off-shell nucleon momenta ki and /c2 and the on-shell deuteron momentum p = ki + k2 
with p"^ = M^. Both nucleons in the vertex d ^ np are final, whereas, as a convention, the 
index a in the matrix corresponds to the final nucleon and the index j3 corresponds 
to the initial nucleon. It is therefore convenient to separate from at the right the 
charge conjugation matrix Uc- 

In order to construct the Bethe-Salpeter function ip^^ , we will start with the "initial" 
expression a^ + b{k-^)/m and then multiply it by ki and k2 from right and left respectively 
and from both sides. Here ^ = e^{p) is the deuteron polarization vector with p-C, = 0, and 
k = {ki — ^2)72. We thus obtain: 
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(5.17) 

The functions ai,bi depend on the scalars k-p and k"^. One can see from this equation 
that the deuteron Bethe-Salpeter function, ip^^ , contains eight spin structures. This fact 
agrees with another classification |1^, ^ Q using the partial waves "^^^^Lj, where the 
index p corresponds to the so called p-spin, which operates on the positive- and negative- 
energy states exactly in the same manner as usual spin operates on the usual spin states. 
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In these terms, the Bethe-Salpeter function of the deuteron is determined by the following 
eight coupled waves |T7[| : 



J C+ n+ i pe 3 po 



1 po 
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3 pe 



PI 



(5.18) 



The indices " + " and " — " mean "up" (++) and "down" ( ) p-states of both nucleons, 

the indices e (even) and o (odd) mean the triplet and singlet p-spin states. From the Pauli 
principle it follows that the function $ satisfies the symmetry relation: ^ai3{ki, k2,p) = 
^0a{k2,ki,p). Hence, %Ij^^{1,2) = —Uclip^^ {2,l)YUc- To provide the symmetry relation 
in a simplest way, we transcribe ( p.l7| ) in the form: 



9i + 92^ + 93 + 9a — -k2 + 2m) 

— — - -.z 



m 



^ m + k^o , (^2 

+5-5 \ 



m 



96^ + 97 



m 

M 

m 



{ki + m) 



+ 



m 



^(fci + ^2). (5.19) 



gi{—k-p,p'^) for z = 1 — 4, 6, 7 and 



m 

The symmetry implies: gi^k-p^p^ 

95,sik-P,P^) = -95,8i-k-p,p'^)- 

The relation between the Bethe-Salpeter and the light-front wave function is given by 
eg. ( 13.581 ). We show below that p.58| ) together with the deuteron Bethe-Salpeter function 
( ^.19| ), coincides with the deuteron light-front wave function ( |5.2| ). 

We substitute ( |5.19|) in (|3.58| ). The momenta ki,k2 in eq. (|5.19|) are off the mass 
shell. However, after the following substitution in ( |3.58| ): 

ki ^ li = ki ~ ujt/2 + uj(3, k2 —>■ I2 = k2 — ujt/2 — ujI3 

we have to put ki, k2 on mass shell. Due to that, we may use the Dirac equation u{k2){k2~ 
m) = {ki + m)Ucu{ki) = 0. After applying the Dirac equation, the function ip^'^ in the 
integrand of (|3.58|) writes: 



1 



^""'{k, - ujt/2 + (3uj, k2 - ujt/2 - (3uj,p) ^ gi-[{ki) + {uji)(3] + 92^ 

m 

Uiuj - uji)T] + (?4^[(H) + {toi)[3]2[3uj 



+g^-[2{uji)(3 



m 



+g,-[(iu;-uji)(3-{uji)T] 



m 



9& 



rV4) - 9s-[{ki) + {uJ-mujT . (5.20) 



m 



The function g-/ does not contribute. One can see that the structures at (71,2,4,6,8 in (|5.20|) 
coincide with the corresponding structures at v^i-4,6 in ( ^.2|) . The expression {^lj — uC,) 
at gs and g^ in ( p. 20] ) can be transformed to the form ■j5e^,^p^^fj,piykpUJ^ at ip^ in (p. 2]). For 
this aim we decompose ^2(^0) — ojC, )ki in the complete set of 4 x 4 matrices: 



M = k2{iuj - uji)h = A + 575 + C"^7m + + 



(5.21) 



The coefficients are given by: A = |Tr(M), B = |Tr(75M), 2ii^^i, = |Tr(i(T^i.M) and 
= Dfj^ = 0. Sandwiching eq.( [5.2lD between the spinors u{k2) and Ucu{ki), and using 
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again the Dirac equation, we find: 

+ \{s - M^)(x - \){uji) + \{s - M')m'^]U^u^ , (5.22) 

where, as usual, s = {ki + k2)'^ = 4(/c^ + m^). This structure contributes to (p^ and to other 
structures. By this way, we reproduce exactly the six structures in the light-front deuteron 
wave function (|5.20|) . Comparing their coefficients, one can easily find the expressions for 
(fi in terms of the integral of Qi over (3. The arguments of the functions Qi = gi{k'^,p-k) 
are replaced in the integrand by: 

p-k -l(x-^)is-M')+m'p', 

where we introduce the dimensionless variable /?' = /3(a;-p)/m^, and, hence, depend on 
/c^, n-k = —2ek{x — 1/2). This provides the explicit link between Bethe-Salpeter and 
light-front wave functions. 



5.3 Two nucleon wave function in the = 0^ scat- 
tering state 

This relativistic scattering state wave function has been calculated in ref. []7T|. Before 
coming to the relativistic case we remind some definitions for the non-relativistic contin- 
uous spectrum wave function in coordinate and momentum spaces. 



5.3.1 Non-relativistic wave function 

We start with the known S-wave continuous spectrum wave function in coordinate space, 
ipp{r), having the usual asymptotical behavior: 

*V(r)l„„ « <(r) . 51!*!:) + ;£!: = ,..5!!*L±^) . (5.23) 

where / = e*^ sm6/p is the scattering amplitude. The corresponding energy is Ep = p'^/m. 

The Fourier transformB ipp{k) of 'ipp{r) is expressed through the off-shell scattering 
amplitude f{k,p): 

Mk) = f Mr) exp{-tk-r)d'r = '^5{k - p) + , (5.24) 

J p rZ p Z6 



^We use the same notation both for the wave function ^p(r) in coordinate space and for its Fourier 
transform 'ipp{k). The difference is indicated by the arguments. 
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where the S-wave off-shell amplitude f{k,p) satisfies the Lippmann-Schwinger equation. 
The variable k is the argument of the Fourier transform, whereas p is related to the 
eigenvalue Ep. 

Because of the asymptotical behavior (|5.23| ), valid for any short-range interaction, 



the momentum space wave function ( [5 .241 ) has a singular term aX k = p and a quickly 
decreasing off-energy shell dynamical part f{k,p). 

It is convenient to extract the global phase factor from the wave function: 

Mr) = e''^^Vp(^) (5.25) 

and from the amplitude: 

f{k,p) = e''^P^g{k,p) . (5.26) 
The function ifp{r) is then real and has the asymptotical form: 

^,M|_«^,-M = 2!!*!:±^1. (5.27) 

The function g{k,p) is real with g{p,p) = sin 6 /p. 
The Fourier transform of Lp^ thus gives: 

= / expi-tk-r}d'r = %5{k - p) + ^'^l . , (5.28) 
f J ^ p'^ K — p — le 

where / = f{p,p) is the on-shell scattering amplitude given above. Subtracting ( |5.28| ) 
from ( |5.24|) , one finds: 

9ik,p) = ^ + {e-p') HiMr) - Viir)f^^rdr . (5.29) 

p Jo K 

By the relations (|5.24|) and ( p.26| ), g{k,p) determines the continuous spectrum wave func- 
tion in momentum space. The difference, (Pp{r) — ipp{r), decreases faster for r greater 
than the interaction range (exponentially for the Yukawa-type potentials) and the integral 
( ^.29| ) can be easily computed. 

5.3.2 Spin structure of the relativistic wave function in LFD 

This wave function is the relativistic light-front generalization of the non-relativistic np 
wave function in the ^5*0 state. The partial scattering amplitude (scattering through 
J'f = 0+ state) can be factorized as follows: 

^;^V",^=e^^X,..xL,e5o), (5.30) 
where xC^o) is the final ^S'o-state spin function on energy shell normalized to 1 : 

Xa.aXSo) = ^^^Vr'{pi)l,Ujf'{p2) , (5.31) 
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In the non-relativistic hmit in the center-of-mass system of the final nucleons, xi^So) 
turns into its non-relativistic counterpart: 

xCSo) ^ ^wlaywl = q°^.^^<(lX(2) , (5.32) 

where w„ represents the two-component nucleon spinor (see appendix and Ci°^i^_^ is 
the Clebsh-Gordan coefficient. 

The amplitude A does not depend on the spin projections of the final nucleons ai 
and a2. We have separated the phase factor e*"^, where 6 is the ^5*0 phase shift. We expand 
the off-energy shell amplitude A for forming the 0+-state in the most general invariant 
amplitudes |l7T| : 

2muj Am? 



-4 = Ui'^^UcU2ai + ui 



75[/cM2a2 , (5.33) 



uo-p s 

where s = (fci + ^2)^, ui^2 = ^'^^'^(^1,2) and 01^2 are scalar functions. In the representation 
( p.3| ) (without the matrix D^*), the amplitude ( |5.33| ) can be decomposed as follows |]7T| : 

1 X / ia-[k X n] \ ^ 
= ^^'^1 1 ^1 + k J ' 

where the scalar functions gi^2 depend on k and z = cos nk. This procedure is completely 
analogous to the one used above for the deuteron wave function. The relativistic amplitude 
( ^.34| ) contains not only the ^5*0 state with spin zero (the first item), but also a state with 
spin 1 (the second item). Therefore, it can be misleading to use the notation ^Sq. We 
refer here to the J'^ = 0+ state. From the Pauli principle it follows: 

9i,2ik,z) = fi'i,2(fc, -z) . (5.35) 



The relations between the invariant amplitudes in the representations ( |5.33| ) and ( |5.34| ) 
are the following: 

ai = — gi, a2 = z-g2 ■ (5.36) 



5.3.3 Equation for the wave function 

The inhomogeneous equation for the relativistic light-front amplitude has a form similar 
to that determining the bound state wave function ( p.40|) : 



^{Pl,P2]ki,k2,UJT) = IC{pi,P2;ki,k2,UJT) + j^-:^ J J='{pi,P2]k[,k2,UJT') 

dr' cPk' d?k' 

xlC{k[, k'2, lot'; ki, ^2, ujT)S''^\k[ + k'^- ut' -ki~k2 + ujt)— tttt^tt^ • 

(r - ze) 2ek[ 2ekL, 

(5.37) 
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The delta-function, together with the integration on one of the momenta, is here kept 
for symmetry. Substituting here the amphtude (|5.30|) and using the variables k, n in the 
system of reference where "P = we find: 

1 , ia-[k X n] , f AV d^k' . 
^{gi + ^T -92) = -- e''xV-—^ — — . 5.38 

with V = — /C/4m^. The products ^ind AV contain summation over spin indices. 



5.4 Numerical results 



In order to have a first estimate of the relativistic corrections to the deuteron wave func- 
tion, one can calculate the components /i — /g perturbatively, starting from the non- 
relativistic solution us and ud- 



-m 



27r3(4(A;2 + m2)-M2) ^t;: 



ii^'NR^y)aia'V^'J{k',k,n,M){ay)^'^, 



cPk^ 



(5.39) 

where ■j/'^^ is given by eg. ( ^.15|) with the replacement of ip hj ipNR defined by (|5.7] ). The 
approximation (|5.39| ) is based on the fact that due to the small deuteron binding energy, 
the deuteron wave function is concentrated at non-relativistic values of k' and the integral 
in the r.h.s. of (p.l4|) is dominated by the same domain. From eq. (|5.39|) one can find /i_6 
(see ref. |^ for more details). 

The comparison of eq.( p.39| ) with the equation to be used with the Bonn-Q models 
(as written in eq.( [4.9| ) for instance) indicates that they are identical in the limit where 
n dependent terms and boost effects under the integral are omitted. In this limit, the 
light-front wave function should be equal to yjek/m ip^°"'"-~'^ (k) . This one could therefore 
provide an approximate zeroth order in solving eq.( |5.39| ). 

The same approximate solution can be found for the ^Sq scattering state ||7l|]. In 
this case, one replaces A in the integrand of eq. ( p.38[ ) by: 



A 



— u[-f5UcU2gi{k',p) 

TT 



(5.40) 



with gi{k',p) = {ek'/m)g{k',p) and g{k',p) given in ( |5.29| ). 

We use the one-boson-exchange interaction with the parameters of the Bonn-QA 
model 0] . This one was used in earlier calculations |4^, ^ and in view of the exploratory 
character of the present work, we have not used a more realistic model (section [4.3.2| ). 
The non-relativistic deuteron wave function was also taken from ref. 0. For the ^5*0 
scattering state, we use the Paris potential [|7D] and the wave function was taken from 
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Results for the deuteron components |^8| are displayed in figs.|T^-[T5|. All the wave 
functions are normalized by ^/A^, as given by ( |5.12| ). In figs]!^ and the functions 
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Fig. 13. The function /i, in (MeV/c) in absolute value, as a function of /c at a fixed 

value of z = cosnk = 0.6. /i is positive at the origin and changes sign twice. The dashed 
line is the non-relativistic S wave function calculated with the Bonn potential, and the 
solid hne is the relativitic result calculated in perturbation theory (see text). 
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□ 



Fig. 14. The same as in fig. |T^, but for tlie function /2. 



Tlie wave functions gi^2 [^1 for the ^5*0 scattering state, for Ep 
MeV/c), in (MeV/c)~^, are shown in fig. |16. 



1.5 MeV (p = 37.52 



In view of the results given in refs. |7T|, one can make the following remarks 



1. In the region k > 0.4-0.5 GeV/c, the components /s of the deuteron wave function, 
and g2 of the ^5*0 scattering state are comparable to the non-relativistic components. 
They are dominant above k = 0.5 GeV/c. These components are calculated here 
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Fig. 15. The relativistic components /i — /e, in (MeV/c) as a function of k, Bit z — 0.6. 
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Fig. 16. The relativistic wave functions of the ^Sq scattering state, in (MeV/c)~^, 
for Ep = 1.5 MeV, at z = 0.6 (sohd hne). The dashed hne is the non- relativistic ^Sq 
scattering state wave function. 
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with the pseudo-scalar representation at the irNN vertex. We shall come back to 
this point in the next chapter. 



2. The functions /i, — /2 and gi are close to the input non-relativistic wave functions 
in the region < k < 0.5 GeV/c. They are however strongly modified for larger 
momenta. This indicates that relativistic corrections are important in this region. 
Some of the differences at low k, for /2 for instance, may disappear once the low en- 
ergy A^A^ scattering data or the asymptotic normalization, Ad and As, are required 
to be produced. 

3. Calculating the phase shift by means of the relativistic function, gi{p,p) = sin5/p, 
we get S = 48° in comparison with the value 61° found from the non-relativistic input 
Paris wave function. The deviation of 13° indicates that the relativistic effects on the 
low-energy phase shift are sizeable and advocate for a new fit of the parameters in the 
relativistic kernel. This deviation shows the actual uncertainties of the relativistic 
kernels. 

4. The relative contribution of the relativistic components /a-e to the normalization 
integral is of the order of 1%. 

5. At /c of the order of 1 GeV/ c all the components of the wave function strongly depend 
on z. This is in contrast with the components /i,/2 and gi in the non-relativistic, 
small k, region. 

6. The n dependence of the wave function is a property of the wave function related to 
the off-energy shell amplitude. At k = p the scattering amplitude should not depend 
on n, otherwise it would violate rotational invariance. In practice this dependence 



could appear due to approximations. It was found in that, at /c = p ~ 37.5 
MeV/c, g2 is negligible compared to gi and the z-dependence of gi is of the order of 
10"'^. These facts indicate that in this respect these calculations are self-consistent. 

The existence of the relativistic extra components /s and g2 which dominate already 
at moderate momenta may seem surprising given the success of the non-relativistic phe- 
nomenology in the two nucleon systems. This fact is indeed related to the dominance of 
meson exchange currents in the deuteron electrodisintegration amplitude. It shows that 
the calculation of the light-front wave function omitting its n-dependence cannot be con- 
sidered as realistic. We shall come back to this point in section |7.3.1j , where the physical 
significance of these components will be investigated in leading order in a 1/m expansion. 
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Chapter 6 

The electromagnetic amphtude 



The physical amphtude and electromagnetic form factors should not depend on a partic- 
ular choice of orientation of the light-front plane. This is however only the case in a given 
order of perturbation theory or in an exact calculation. In any approximate calculation, 
the light-front electromagnetic amplitude depends on the orientation of the light front, 
i.e., on the four-vector uj. This dependence is non-physical for any on-shell amplitude. 
We will show in this chapter, for systems with spins 0, 1/2, 1 and transitions between 
them, how to extract the physical form factors from the cu-dependent electromagnetic 
vertex. The covariant formulation of LFD will prove to be very powerfuU in resolving the 
ambiguities arising from the light-front orientation. 

6.1 Factorization of the electromagnetic amplitude 

The form factors are experimentally extracted from the scattering amplitude of a charged 
particle (usually the electron) from a bound system. One of the contribution to the 
amplitude is shown graphically in fig. 

For simplicity, we consider a two-body system to start with, but all our results are 
obviously applicable to many-body systems. In the Feynman approach, the amplitude of 
this process is factorized in the product of the electromagnetic vertex of the projectile, 
the propagator of the virtual photon and the electromagnetic vertex of the bound 
system. This latter is then expressed through the physical form factors. The complete 
set of light-front diagrams should, of course, be factorized and reproduces this structure. 
However, the separate amplitude represented in fig. ^ does not in general factorize in 
this formalism. Indeed, the expression for this amplitude, corresponding to the photon 
propagator and the adjoining spurion lines has the form: 

M= [{■■■)d iuip' - p)) 6 Up' -p + ut;- ujT.f) "^^^ . "^""^ . . (6.1) 

The dots represent the parts of the amplitude which contains the matrix element of the 
current Jp. They will be discussed in the next section. 
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p p 

Fig. 17. Particular amplitude corresponding to the scattering of an electron from a bound 
system. 



For definiteness, we suppose that uj-{p' — p) > 0. Integrating eg. (|6.1|) over dr^ by 
means of the (5-function, we get, with q = p' — p: 



M 



dT2 



[2{uj-q)T2 - g2 _ ie] (rg - ie) 
This amplitude is not a product of the term with the form factors. 




(6.2) 



P ki 

Fig. 18. Electromagnetic vertex of a bound system. 

However, as seen from ( |6.2|) , the factorization can be obtained under the special 
condition 

u-{p' — p) = uj-q = , (6.3) 

implied as the limit uj-q +0. Under this condition, the amplitude ( |6.2| ) obtains the 
form M ~ f{(l)/(f'- The integration over dT2 in (|6.2| ) is absorbed by the form factors. 
The corresponding electromagnetic vertex can then be represented by fig. ITBI with the 
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off-mass-shell photon momentum ^ 0. In the standard formulation of LFD with 
oj = (1, 0, 0, —1), the condition (|6.3|) turns into the usual condition = 0. This condition 
can only be achieved for < 0. 



6.2 Extracting the physical form factors 



The problem of extracting the physical form factors from the elementary electromagnetic 
amplitude in the LFD is rather general and appears for a system with any spin. As 
it was discussed in chapter 0, the sum of all the on-energy-shell amplitudes in a given 
order of perturbation theory does not depend on the four- vector uj, whereas any given 
amplitude can be tu-dependent. This is so also for the on-energy shell electromagnetic 
amplitudes. The amplitude corresponding to a separate diagram can indeed depend on 
u. The diagram of fig. O for instance is not the only one contributing to the form factor. 



Other contributions are shown in fig. IS. 





Fig. 19. Amplitude which cannot be expressed through the bound state wave function. 



The matrix element of the exact current, Jp = {p'\Jp{0)\p), can be represented as 
the sum of all possible contributions: 



(6.4) 



where the superscripts is a short notation for indicating the various levels of approxima- 
tion. In the case of the deuteron form factors for instance, the terms in ( |6.4| ) correspond to 
the impulse approximation, meson exchange current contribution, etc. Any contribution 
in (|6.4| ) consists of the sum of two parts: a u independent one denoted by F^"'\ and a u 
dependent one denoted by B'^^cu): 



J(n) ^ pin) ^ ^(r. 



(6.5) 



The explicit form of eq. (|6.5| ) for systems with spin 0, 1/2 and 1 is given below in sect. 
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Since the matrix element of the exact current does not depend on uj, the cj-dependent 
parts cancel each other: 

B^^\co) + Bf\u) + . . . + 5W(a;) + . . . = , (6.6) 

and, hence, the matrix element is determined by the sum of the cu-independent parts only: 

= + + . . . + Fj") + • • • . (6.7) 

It follows that the problem of restoring the independence of any physical amplitude 
on 10 is reduced to the separation, and omission, of the "pure" u-dependent terms. 

We emphasize that this definition of the form factors, after separating out u- 
dependent contributions, still leads to approximate form factors. However, our procedure 
guarantees that the form factors do not receive any spurious, non-physical contributions. 
Since these spurious, cu-dependent, contributions can be as large as the physical ones, 
in particular at high momentum transfer, our procedure is the most adequate when one 
wants to draw any conclusion on a particular model, or to compare various relativistic 
approaches. Moreover, it leads to an unambiguous determination of the form factors, 
independently of the orientation of the light front or of the choice of the component of 
the current which is used to calculate the amplitude. This is not the case in the usual 
formulation on the light-front plane t -|- z = 0, as it is already known for the deuteron 
form factors |^5|, ^ |. 



In the case of spin and 1 systems, the criterion of irreducibility for the cj-dependent 
contribution is rather evident (though with a subtle point discussed in sect. |6.3.3|) . It can 
be read directly from the decomposition of the current: the terms proportional to uj do not 
contribute to the physical form factors. In the case of spin 1/2 particle, the unambiguous 
separation of the purely w-dependent part requires some care and is discussed in section 

lo:^ . 

In a more general language, the extraction of the physical form factors can be 
explained as follows. The physical electromagnetic amphtude Jp = {p'\Jp{0)\p) does 
not depend on the hypersurface where the state vector \p) is defined and, hence, does 
not contain any term proportional to ujp, provided the current operator Jp(0) is the full 
operator. This current operator has to contain the interaction. Indeed, the commutation 
relation between Jp(0) and the four- dimensional angular momentum J^p is similar to 
( p.33| ), since it is universal for any four- vector: 

-[Jp{0), Ja(3] = gppJM - 9paJp{Q) . (6.8) 

% 

Since Jq,^ in (|6.8|) contains the interaction, the current JpifS) has to contain the interaction 
too. The consistency of the transformation properties of the current and of the state vector 
ensures the independence of the matrix element of the current Jp(0) on the quantization 
plane orientation. 

In practical calculations however, this consistency is violated twice: i) the diagram 
of fig. ITB| corresponds to the free current and ii) the state vector is approximated by its 
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two-body component. Due to these reasons, the vertex Jp obtains spurious contributions 
proportional to ujp. They are however completely explicited in our covariant formulation, 
while they are hidden in the usual formulation of the LFD with t + 2; = 0. 



6.3 Light-front electromagnetic vertex 

The calculation of any amplitude with composite systems in the initial and final states 
is done according to the rules of the graph technique discussed in section |2.2| , with the 
coupling constants at the vertices replaced by the vertex functions V and introduced 



in section for the final and initial systems respectively. We detail in this section how 



to perform the calculations for spin 0, 1/2 and 1 bound states. 
6.3.1 Spin system 

Let us consider the structure of the vertex of fig. IT3 for a system of spin 0. It has the 



general form: 

Jp = {p + P'),F{Q^) + a;^i^±^i?,(Q2) , (6.9) 



in which the factor {p + p'Y / {2uj-p) is separated for convenience. The tilde on top of 
the electromagnetic amplitude Jp here and below indicates that this amplitude is an 
approximate one and contains cj-dependent non-physical contributions, as compared to 
the exact physical amphtude with no tilde. As usual, we have = —q^ = —{p' — p)"^. 
The invariant functions F and Bi could depend in principle on uo-p and uj-p' . However, 
the four-vector uj is defined up to an arbitrary number, and, hence, the theory is invariant 
relatively to the replacement uj — > atu, where a is a number. The form factors F and Bi 
can therefore depend only on the ratio uj-p' joj-p. But since uj-q = 0, we have uj-p'/u-p = 1, 
and the functions F and Bi (and similar functions for the case of spin 1/2 and 1) depend 
on only. 

The general physical electromagnetic amplitude of a spinless system is given by: 

Jp ^ (pVp(O)b) = ip + p')pFiQ') . (6.10) 



The main difference of the amplitude ( |6.9| ) with respect to ( |6.10| ) is the presence of an 



additional contribution, proportional to Up. Note that the expression (|6.9| ) satisfies the 
conservation of the current Jpq^ = since u-q = 0. This conservation is automatically 
fulfilled in this particular case due to the simplicity of the structure (|6.9| ), but it cannot 
be imposed a priori. We will see below that for a system with spin 1, some cu-dependent 
terms violate indeed current conservation. 

In the simple case where the constituents are also spinless (Wick-Cutkosky model), 
the amplitude indicated in fig. ^ reads: 

Jp = j [{p + i^T2-ki)p+{p' + UTi-ki)p]e{u\p-ki))5{{p + UT2-kif -rn^^ 
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X T1T2- — — 7-—^ , 6.11 

(ri - te) {t2 - It) {2Txy 

where the vertex functions ri,r2 depend on the momenta in the vertices 1 and 2. In- 
tegrating by means of the (5-functions over dri dT2 dko, and expressing F through ip by 
means of eq. (|3.39|) , we get: 



where ri = [m^ — {p — ^1)^] / (2u;-(p — ki)) and T2 = [m^ — {p' — ki^] / {2uj-{p' — ki)) 

From ( |6.12|) one clearly sees that Jp contains terms proportional to Up and, hence, 
has the structure of eq.(|6.9[). However, to avoid any misunderstanding, we emphasize that 
the contribution ujti + ujt2 is not the only source of u;- dependence. Even in the absence 
of this term and in the case where the wave function ip does not depend on n, the term 
proportional to ujp survives because of momentum conservation at the electromagnetic 
vertex. 

In the spinless case, the prescription to extract the physical contribution to the 
form factor, F^Q"^), from the a;-dependent amplitude Jp can be obtained immediately by 
multiplying both sides of eq.( p.9|) by Up. We thus get: 

FiQ') = ^ ■ (6.13) 



With ( |6l2D , ( |67[3D we obtain: 

m'') = T^Ji^{R\.A^m^-^^)\^)f^^y (6.14) 
\2i^Y J 2x(l — x) 

We have represented here, and in the following, the four-momentum transfer g by g = 
(goj ^1 %) with A-cJ = and q\\ is parallel to Q. Since uj-q = 0, we have = —q^ = A^. 

In the usual light-front formulation, with u = (1, 0, 0, —1), eq.( |6.13| ) corresponds to 
expressing the form factor through the J+ component. This is well known, and eq. (|6.14|) 
has been found in ref . ||8^ . However, this procedure cannot be extended to the calculation 
of physical form factors of systems with total spin 1/2 and 1, as we shall see in the next 
sections. 

One can similarly calculate for comparison the non-physical form factor Bi. 

, 2^ ^ f {Rl + - M\i - xf - xR^-K + gV2) 

' (271)3 7 (4M2 + g2) 

xV^(^l, x)^{{R^ - xA)\ x)^^^^^. (6.15) 

Here and below we use M for the total mass of the system, m being used for the mass of 
its constituents. 
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For a system with small binding energy |eft| -C m, the wave function is mainly 
concentrated at a; ~ 1/2, ~ '"^Icfel- For such a system, the function Bi at <^ 
has smallness ~ le^l/m in comparison with F or ~ Q'^/M'^, if Q'^/M'^ > le^l/m. However, 
at Q"^ > 711? the function Bi becomes comparable with F . We will illustrate this in section 
|7.1| in the Wick-Cutkosky model. 



6.3.2 Spin 1/2 system 

The electromagnetic physical amplitude for a spin 1/2 fermion is determined by the two 
standard form factors: 

Jp = Fi{Q'^)u'ypU + ^-^^^^-^u'apuq^u = u'T pU . (6.16) 

Here u' = u°''{p'), u = u°'(p) are the spin 1/2 fermion bispinors and M is the fermion 
mass. The charge (Ge) and magnetic (Gm) form factors are expressed through Fi and 
F2 as follows: 

Ge = Fi-^F2, Gm = F, + F2. (6.17) 

The covariant light-front spin 1/2 electromagnetic vertex was firstly considered in 
ref. 1^ and investigated in detail in ref. ||88| (see also ref. |l89|). 



Like in the case of spin 0, the dependence of the corresponding amplitude on the 
extra four-vector uj increases the number of independent terms in the decomposition of 
the amplitude. Besides the two structures entering in eg. ( |6.16| ), one can construct the 
following three structures: 

1 , 1 , 1 , , , 

— u uju Pp , — u u Up , 7 tttM uu Up , (6.18) 

u-p u-p y^-p) 

where P = p + p' , lu = ^^^7^. 

Other structures, like non-gauge-invariant ones: uuqp, uuuqp, u'apyU u'^ and also 
iu''~^^u Cppu-yP^q"^^ (not independent from the previous three), contradict T-invariance. 

At first glance, it seems that all the structures in (|6.18|) exhibit no tu-independent 
parts, since, like in the previous section, they are proportional to u/u-p in the first and 
second degree. This is indeed true for the last two structures which are explicitly propor- 
tional to Up/ u-p. 

This is however not the case for the first structure u'uuPp/u-p which is directly 
proportional to Pp, as one of the physical amplitudes (after Gordon decomposition). It 
could in principle contain a part which has to contribute to the physical form factors and 
therefore cannot be simply omitted together with the full structure. This becomes evident 
in the particular situation where p' = p. In this case, due to the identity u{p)'yp_u{p) = 
uupp/M, we get uuuPp/u-p = uuPp/M, i.e., a non-zero cu-independent contribution only. 
It indicates that this part exists also when p' 7^ p. 

In order to construct the appropriate structure, we demand that it is orthogonal to 
the physical contribution given in ( |6.16D . Since there are only two physical contributions. 
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any other contribution proportional to Pp and orthogonal to them is completely indepen- 
dent and therefore non-physical. Thus, eliminating this structure from the vertex, we 
shall not loose any physical contribution to the form factors. 

For the two structures u'U^u and u'V^u^ we define here the orthogonality by the 
following condition: 

Tr[{p' + M)UP{p + M)Vp] = (6.19) 

with Vp = JoV^Jq. This definition is motivated by the procedure, detailed in the next 
subsection, to calculate the form factors. 

The only alternative to the structure u'Cju/u-p is the term u'{ijj/uj-p — a)u with 
arbitrary a, since all the possible scalar structures sandwiched with the spinors are pro- 
portional either to u'Cju or to u'u. The value of a is then determined by the orthogonality 
conditions. This gives: 

a = — — = — , wliere ri = . (6. 20) 

P2 (1 + ?^)M' ' 4M2 ^ ' 

Indeed, at this value of a,0 the following two orthogonality conditions are simulta- 
neously satisfied: 

Tr 



Tr 





AM 


oj-p 


p2 


Cj 


AM 


^u-p 


p2 



{p' + M)iaP''q,{p + M) — Pp = . (6.21 




The reason why a unique value of a satisfies two conditions is a consequence of the fact 
that the expression {uj/uj-p — 1/((1 + vi)M)), being substituted in ( |6.21| ), is proportional 
to e, where e is the following four-vector: 

^ - 2P^) . (6.22) 



2uj-p V P2 
orthogonal to all available four momenta: 

e-p = e-p' = e-q = . (6.23) 

For convenience, we shall normalize e by the condition = — 1. Note also that the term 
u'{u/uj-p — 1/((1 -|- ri)M))u is proportional to u'eu. 

The two last structures in (|6.18| ) are not orthogonal to 7p and iap^q^ . However, they 
are proportional to ujp, and, hence, do not contribute to the cu-independent structures 
associated with Pi,P2 and vice versa. In this respect, the situation does not differ from 
the spin and 1 cases. We emphasize that in this way the structures proportional to 
ijjpjijj-p and Pp are separated by different formal criterions. 

Another way to separate the physical and non-physical structures is the following. 
One can simply subtract the two physical structures with undefinite coefficients from 
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u'Cju Ppl (uj-p) and then multiplying at the left by u' and at the right by u and taking 
sum over polarizations get a 4 x 4 - matrix Mp = {p' + M)Up{p + M). This matrix can 
be decomposed in the full set of 4 x 4 -matrices, similarly to ( |5.21| ): 



Mp = Ap + Bp-f5 + Cpp^^ + Dppi^^-i^ + Ep^,ia^\ 

Now the matrix Mp can be analyzed in terms of the tensors Ap^ Bpp, etc. We require them 
to be irreducible or proportional to uj. The tensor Ap turns out to be proportional to Pp, 
and therefore is set to zero. The tensor Bp equals to zero. The trace of Cpp is set to zero 
in order to construct an irreducible tensor. The antisymmetric part of it is proportional 
to UJ. These two conditions fix the coefficients, and we obtain the final structure Up which 
exactly coincides with [u/uj-p — 1/((1 + ri)M)]Pp. With these coefficients, the tensor 
turned out to be proportional to u, and the tensor Eppi, is irreducible and proportional 
to UJ for the antisymmetric part. In this analysis, all the three nonphysical structures in 
( |6.18| ) are treated on an equal footing, since the tensors obtained after decomposition of 
the last two are proportional to Up and therefore are non-physical. 

It should also be mentioned that in the Breit system, where P = 0, the term u'eu 
obtains a particularly transparent form. In this system, the spatial part of e is the unit 
vector directed along uj: e = uj/ujo = n, whereas cq = 0. Therefore, the term u'eu becomes 
proportional to n: 

u'eu = —u' {p')n-^u{p) , (6.24) 

with n orthogonal to both p and p' = —p. 

In an arbitrary system of reference one can find in u'eu an cu-independent part. 
Omitting this full structure, we inevitably omit also this tu-independent part. Its form 
depends on the criterions we used above (orthogonality and irreducibility) . These criteri- 
ons, strictly speaking, are however a matter of convention. They are imposed as natural 
and reasonable, but are not mathematically derived from some more general principles. 

The above consideration shows that we achieve (with the structure u'euPp) optimal 
separation of nonphysical and physical contributions. 

Collecting the various possible contributions, we obtain the following expression for 
the spin 1/2 light-front electromagnetic amplitudJ 



Jp = u TpU , 



2 



iF2 f UJ 1 \ M M 

T^^P^l'' + B^i — - 7T— ^ P, + B,—ujp + B,—- 
2M ycu-p {l + r])MJ uj-p (cu-p) 



(6.25) 



The electromagnetic vertex ( |6.25| ) is gauge invariant since Jpq'' = (with u-q = 0). As 



mentioned above, the possible non-gauge-invariant terms are forbidden by T-invariance. 



similar expression for Tp was given in | |85[ |. In this reference however, as weU as in [g9|, the structure 
LuPp/uj-p was used instead of the structure proportioanl to Bi in ( 6.25| ). In this case, the form factor 
GsiO) does not coincide with the normalization integral. 
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The physical form factors 



We can now express the physical form factors Fi and F2 through the full vertex function 
fp. We multiply Jp by [u"' {p')Yu'^{p)]*, [u''' {p')iaf"'qi,/ {2M)u''{p)]* , etc. and sum over 
polarizations. We thus obtain the following quantities: 



ci = Tr[OpY] , C2 = Tr[0piaP''q,]/{2M) , C3 = Tr[Op{u/u-p - 1/(1 + r])M)]P'' , 
C4 = Tr[Op]ujPM/uj-p , C5 = Tr[OpUj]ujfM^/{uj-py , (6.26) 

where 

Op = iP' + M)fp{p + M)/{AM^) . (6.27) 

By this way, we get a linear system of five equations for Fi, F2, Bi_^^ with the inhomo- 
geneous part determined by Ci_5. The way of calculating ci_5 through traces with the 
matrix Op just corresponds to the operation used in the definition of the orthogonality 
condition ( |6.19|) . The values of Ci„5 can be calculated in specific models for the structure 



of the composite system, as shown in the next chapter. 

Solving the above mentioned system of equations relative to Fi and F2, we find: 

Fi = 4(^1^)2 [(c3 + 4c4-2ci)(l+r7)+2(ci + C2)-2c5(l + r/)2] , (6.28) 

7 [(c3 + 4c4-2ci)(l+r/)+2(ci + C2)-2(c5 + C4)(l+r/)2] .(6.29) 



4r/(l + 77) 



In spite of 1] in the denominator in eq.( |6.29|) and below, there is no singularity at = 0. 
From eqs.( |6.17|) , we can now easily obtain Ge and Gm- 



Ge = C4/2 , (6.30) 
Gm = 4^(/^ ^) [(C3 + 2c4 - 2ci)(l + r]) + 2{c, + C2) - 2c,{l + r])^] , (6.31) 

with ci_5 given in ( |6.26|) . These formulae are quite general in the light-front approach. 
They are applicable to any model of the spin 1/2 fermion structure. 

Comparison with other approaches 

In the usual formulation of LFD on the plane t + z = 0, the form factors of spin 1/2 



systems are found from the plus- component of the current (see, e.g., [|^, 0, ^ p2l), 
i.e., with our notation, from the contraction of Jp in eq.( |6.25D , with ujp. This contraction 
eliminates the contributions ci/i?2,3, but does not eliminate the term with Bi. The form 
factors F[ and F2 deduced in this way are thus given by: 

J- = -'[^^^ + §---^^^ + 25^(^-(r^)]- 

= u'[F[lp+'^^,.q'']uoj^ . (6.32) 
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where 

F[ = F, + , F2' = F2 + -—Bi . 6.33 

l + rj l + rj 

Equation (|6.32| ) shows that the structure of J-uj (or 7+) indeed coincides with the standard 
representation of the nucleon electromagnetic vertex, eq.( |6.16|) . However, F[ and in 



( |6.32| ) are not the physical form factors, but their superposition (|6.33|) with the non- 
physical contribution Bi. 

The expression for Bi can be found from the above mentioned system of equations, 
leading to: 

Bi = J (C3 + 4c4 - 2ci)(l + 77) + 2(ci + C2) - 4c5(l + ^if] . (6.34) 

Together with eqs.( |6.33| ), one has for F[ and Fg: 

F[ = C5/2, = (C5 - C4)/(2r/) . (6.35) 
From eqs. (|6.35|) and (|6.17|) one can find the corresponding G'^ and G'^f. 



G'e = Ge, (6.36) 
G'm = Gm + 2Bi = [C5(l + r^) - c^]/{2ri) . (6.37) 

The difference between the two approaches - the first one separating out all the uj- 
dependent terms to define the physical form factors, and the second based on the compo- 
nent J-|. of the electromagnetic current and incorporating non-physical item with Bi - is 
in the magnetic form factor. This difference is of relativistic origin and has to disappear 
for form factors of extremely non-relativistic systems (i.e. in the limit where the average 
internal momenta go to zero) at finite . It also disappears in an exact calculation 
leading to all i?j = 0. We shall estimate this difference in a simple quark model for the 
nucleon structure in the next chapter. 



OfF-energy shell effects in electromagnetic form factors 

Though all the four-momenta are on the mass shell, this does not mean that there is no 
off-shell effects in the form factors at all. If, for example, the nucleon form factors are 
used in the calculation of the deuteron form factors, they are part of a larger diagram and 
contain external spurion lines. They are off-energy shell. Such a form factor is shown in 

fig. 

In the standard formulation of LFD this corresponds to the fact that the minus- 
components of the nucleon momenta and the photon momentum are not related by the 
conservation law. Hence, the off-energy shell form factor can depend on the minus- 
components of the momenta, in addition to its Q^-dependence. In the covariant formula- 
tion of LFD the graph for the off-energy shell form factor contains external spurion lines. 
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p p 

Fig. 20. Off-energy shell form factor. 

The off-shell effects are parametrized in terms of the corresponding spurion four-momenta. 
Besides Q^, the off-shell form factor F depends therefore on two variables p3| : 



F = F(Q^s,s'), where s = {p - urf , s' = {p' - uor'f . (6.38) 

Since u-q = 0, no other independent scalar products can be constructed. 

The spin structure of the off-energy shell electromagnetic vertex differs from (|6.25|) 
by three extra items: 



„ iFo ( Cj I \ M 



where all the form factors depend on Q^, s, s' and are symmetric relative to the permuta- 
tion s <^=^ s'. The last three spin structures have been already mentioned after eq.( |6.18 



Their opposite T-parity is corrected in (|6.39| ) by the factor [s' — s), changing the sign due 



to permutation of initial and final momenta. We emphasize that, in contrast to ( |6.25| ), 
the a;-dependent items in (|6.39|) appear not because of approximation, but due to the off- 
energy shell continuation. Of course, in an approximate calculation of the off-energy shell 
form factors these item appear due to both reasons. The form factors i?i_3 in an exact 
electromagnetic amplitude turn into zero on the energy shell r = r' = 0— >s = s' = M^. 

The off-shell electromagnetic vertex of a spin-zero particle differs from eq.( |6.9|) by 
the extra term proportional to (s' — s)qp/M'^. 

The extra form factors and the dynamical dependence of all the form factors on the 
off-shell variables s, s' are usually unknown and are neglected in almost all calculations. 

6.3.3 Spin 1 system 

The electromagnetic form factors of a system with spin 1 (the deuteron, for instance) are 
defined by the following decomposition of the electromagnetic vertex: 

(AVpIA) = (p'){Pp[-^i(gV'^ + -^2(g^)|^] +^i(g^)«g'^-^?,V)}e^(p) 

= et\p')Ti:^el{p). (6.40) 



93 



Here e^(p) is the deuteron polarization vector, p and p' are the initial and final deuteron 
momenta, A and A' are the corresponding helicities, P = p + p' and q = p' — p. Charge 
(Fc), magnetic {Fm) and quadrupole (Fq) form factors are expressed through jFi, JF2 and 
Qi as follows: 



Fc 

Fm 
Fn 



-^i-y[-^i + 6^i-^2(l + r/)] 



(6.41) 

(6.42) 
(6.43) 



where rj = Q'^/AM'^. 

In the absence of polarization, the cross section of electron-deuteron scattering has 
the form: 

1 



da 

In 



da 

dn 



A{q')+iaxi' -e B{q' 



(6.44) 



where 9 is the electron scattering angle in the laboratory system, and {da/dVL)^ is the 
Mott cross section: 

a2cos2(e/2) 



' da\ 

dn I 



4^2 sin^(^/2) 



(6.45) 



The functions A{q^), B{q^) are expressed through the form factors as follows: 



A{q') 



L{l + r^)F'^{q^) . 



(6.46) 
(6.47) 



Unpolarized experiments are able to separate A and B, but do not allow one to 
extract all three form factors. To do this, one should measure a polarization observ- 
able. Due to its spin one, the deuteron can have, besides vector polarization, quadrupole 
polarization. The corresponding observable is T20 given by: 



T2o[A{q')+tan^-eB{q^] 



1 

V2 



9 



1 



r]FcFQ + -r]'F' + -r] [1 + 2{1 + r]) tan' -6 F^ 



M 



(6.48) 



The deuteron electromagnetic properties are reviewed in ref. p4 



The decomposition ( 6.401) is valid for the full electromagnetic current operator. Con- 
sidering the one-body current J^, as indicated in fig. we get the following analytical 
expression: 



m 



(2vr) 



Tr[<j)'^'{k'^ + m)Vp{k2 + m)<i)''{ki - m)] 



e{uo-ki)e{uj-k2)e{uo-k'^) d^k^ 



uo-ki/uo-py 



(6.49) 
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where Tp is the electromagnetic vertex of the nucleon given by eq.( |6.16| ) (with M replaced 
by m). The matrix is given by eq. (|5.2|) , while the matrix (f)"^ is obtained from (p^, by 
the replacement ^2 ^2- This expression is derived using the graph technique presented 
in chapter It is similar to the calculation indicated in sect. |6.3.1| for spinless particles. 
The difference with ( |6.11D , ( |6.12D is due to the non-zero spins. We explain below, how 



eq. (|6.49| ) is derived. 



According to the remark in the end of sect. 2.2.3| , we start in the graph ^ from the 



outgoing deuteron line (i.e., from the vertex A^A^ — > d'), follow the line in the direction 
opposite to the orientation of any nucleon line and reach its end, i.e., the vertex d NN. 
Then we start again from the vertex A^A^ — >■ d, and follow another nucleon line, also in 
the direction opposite to its orientation. We thus obtain: 



x{ki + m)f3a. (6.50) 



The first row in ( |6.50| ) corresponds to the upper nucleon line in the diagram, the second 



one to the lower line. We attach the deuteron wave function to the upper nucleon line. The 
first row contains the factor [7o0''^f^c7o]^ which originates from the complex conjugated 
final wave function of the deuteron: 

[e';{p')u{k',)<j>'^VMh)r = ef{p')u{h)[^o<P'^U,lo]^u{k',). 

We keep in ( p.5(J| ) the matrix indices a and (3 explicitly. Since both nucleon lines are 
equivalent, the order of indices (3, a is the same. This means that one of the factors in 
( |b.5(J| ) (we take the second one) is the transposed matrix. With the wave function {^72) 
we get: 

and, hence, obtain the factor: 

that gives the trace in eq. (|6.49|) . 

Like in the case of a spinless system, the light-front vertex ( |6.49| ) depends on the 
four-vector uj. Its expansion in terms of form factors contains therefore extra terms. This 
expansion has the general form [^: 

{XW\) = ^efip')j;:^etip), where = + B^/ (co) , (6.51) 

where Tj^" has the structure given by eq.( p.40| ) and is determined by the physical form 
factors, J-'i, T2 and Q\. We recall that tilde is used to distinguish approximate amplitudes 
and matrix elements, containing tu-dependent non-physical contributions, from the lo- 
independent ones. The tensor B^^ contains the uj dependent terms: 



M2 



Big^" + B2^ + B^M'- + B, 



M2 ' {uj-pY ' 2uj-p 
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(6.52) 



q^'uj" + q^'u'' 
2uj-p 

Here P = p + p' and Bi, ...,Bs are invariant functions. The tensor B'^'^ as well as Tj^'^ is 
symmetrical with respect to the simultaneous permutations p ^ p' {q ^ —q) and fi ^ u. 
The structures at Bj and Bs in eq.( |6.52| ) do not satisfy the conservation of current, that 



is qPJ^ ^ ^ 0. 

The total number of independent spin structures can be calculated without their 
explicit construction. Owing to the cu-dependence, the helicity indices in the vertex are 
not constrained by the equality A = Ai + A2 (A is the photon helicity, Ai, A2 are the hadron 
helicities). Therefore we have 3x3x3 = 27 initial matrix elements. P-invariance reduces 
this number down to 14. T-invariance (or C-invariance in annihilation channel) eliminates 
in addition 4 structures and finally 10 independent structures remain. This calculation 
automatically takes into account the conservation of the electromagnetic current which 
would lead to a relation between B^ and Bg. However, cj-dependent terms in eq. (|6.52|) do 
not obligatory satisfy this condition. The functions Bj and Bg are therefore independent 
and we obtain 11 structures (or invariant form factors): three of them are in T^^ and the 
remaining eight ones are in B^^. 

From the point of view of the usual light-front approach with t + z = 0, the depen- 
dence of the vertex ( |6.51D on the orientation of the light front means lack of relativistic 



covariance. Hence, the formula ( |6.40D , based on this covariance, is not any longer valid. 
Therefore, the number of independent matrix elements in eq.( |6.51|) is not reduced to 3, 



but is larger. The increase of the number of matrix elements just corresponds to the 
increase of the number of independent spin structures in eq. (|6.52|) . In the formulation of 



the light-front with t + z = 0, the contribution of these structures is not separated out 
from the physical form factors. In the covariant formulation, the dependence on the light 
front is parametrized explicitly by eqs.( |6.5lD , (|6.52|) . 



However, like in the spin 1/2 case, and in the absence of explicit evaluation of 
the a;-dependent terms, one can see some ambiguity in the procedure: what has to be 
considered as a non-physical contribution? Consider for example the term proportional 
to B^: PpC*^ {p')uj ^ujye^{p) / {oj-pY . Multiplying it by e', e* and summarizing over A', A, we 
get: 

P V e^7»')e*^7»')^^^e" r»)e*V») - P ( - ^^'"'^ - "^^^'^ + 

' h (^-p)' ' ~ V i^-py (^■p)M' {oo-p)M^ + j • 

(6.53) 

It contains, apart from tu-dependent contributions, the term p'^^/M'^. This one cannot 
be removed from (|6.53|) and be included in the physical part, since it is required by the 



transversality of the tensor (|6.53| ). The orthogonality condition, similar to the spin-1/2 
case (obtained by the generalization of eq.( |6.19| )), results in the charge form factor Fc 
which does not coincide at g = with the normalization integral, in contrast to the spin- 
1/2 case. We therefore do not impose this condition in the spin-1 case, and define the 
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non-physical structures as proportional to the first and higher degrees of cu, including the 
terms where uj is contracted with the polarization vector. After that, the physical form 
factors can be unambiguously separated from the spurious tu-dependent contributions. 
This allows one to obtain the contribution to the physical form factors in terms of the 
electromagnetic amplitude in an explicit form by a generalization of the formula ( [6.131 ). 
This approach is confirmed by comparison with the Bethe-Salpeter solution in the Wick- 
Cutkosky model (see sect. |7.1.2| ). 



Contributions to physical form factors 

The explicit expressions for the form factors are obtained by solving eq. (|6.51|) relative to 



jFi, J^2 and Qi. The procedure is described in ref. |85]. The result is given below: 



-^1 

2M2 



P 2u-p 



9fj.u 



+ P 



2 ^ii^v 



(6.54) 



^ 2{uj-p)q^ 



[9, 




(6.55) 



+ PuUJ 



2{uj-p)q^ 
{uj-p)q^ 



2iu-p) 



{uj-pY 



{uo-p)q^ 



(6.56) 



These expressions determine the electromagnetic form factors. In spite of the fact that uj 
enters the r.h.s. of eqs. (|6.54[ - |6.56|) , these expressions do not depend on uj. 

The form factor JF^, eq.( |6.54[ ) which coincides at = with -F(7(0), does not coin- 
cide with the normalization integral ( |3.31 ). This is due to the fact that the normalization 
is obtained on the basis of the decomposition ( p. 291 ), where the non-physical structure 
proportional to i?, contains g^" besides uj^^uj^ . This is required by the irreducibility of 
this structure and is compatible with the angular condition. In the decomposition of the 
electromagnetic vertex, we have not separated out this term from the terms proportional 
to uj^uj'^ , as we explained above. We believe that the concordance of -Fc(O) with the nor- 
malization integral is beyond the accuracy of the current calculations which omits all the 
Fock components except for the two-body one. In practical calculations of the deuteron 
elastic form factors, the difference is very small. 



Comparison with other approaches 



The deuteron electromagnetic form factors were already calculated in LFD in refs. 

m 



95, 



Wl\ , |100| , |101| |. More generally, spin-1 form factors were considered in refs. 
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B6| , |102| , |103| , |104|| . Using the exact decomposition ( |6.40| ), one can express the matrix 
elements < A'|J_i_|A >= J\i\ in terms of the form factors: 



'A'A 

Ji-i 
Jio 



;i-2r/)J-i-2r/2^2 + 2r/6;i 



(6.57) 
(6.58) 

(6.59) 
(6.60) 



The four matrix elements in ( |6.57[ - |6l60| ) can be expressed through three form factors and, 
hence, are not independent from each other. They satisfy the relation [03]: 



(1 + 2r7) Jn + Ji-i - 2v^Jio - Joo = , 



(6.61) 



which can be checked by direct substitution. This condition is called "angular condition" 
(do not confuse with the angular condition for the state vector discussed in sect. p. 1.3] ). 

In principle, with the exact matrix elements Jyx one can calculate any triplet of 
matrix elements and find the form factors. It is not so however for the approximate ones 
J\i\. The following triplets of the matrix elements were used in the literature: Jn, Ji„i, Jqo 
and Jii, Ji_i, Jio in ref. (the so called solutions A and B respectively); Ji_i, Jio, ^oo 
in ref. [[102|| . For example, taking the triplet Jn, Ji_i, Jqo and solving the system of 
equations (|6.57|) , (|6.58|) and ( |6.60D (with J\i\ J\'\) relative to J^i, one finds: 

-GK-A 



-r-GK-A 
-'2 

rGK-A 



-Jll — Ji 
1 

'2r] 



1 - 2ri)Jn + Ji-i - Ji 



00 



(6.62) 



In refs. [^, |100|| the form factors were expressed through different combinations of four 



matrix elements. 

The matrix elements in the r.h.s. of eqs.( |6.6^ ) calculated in LFD are however not 
given by the decomposition ( |6.40|) , but should be calculated according to eqs.( |6.5lD , ( |6.52|) . 
The matrix elements of ( |6.51|) have the form: 



Ji,o 
Jo,o 



-72^(^1 - + 6^1/2) + ^^2^6 . 
-(1 - 1ri)Tx - 2TfT2 + 2T^gx - 2r]Be 



B5 + Br 



(6.63) 



The matrix elements Jn, Ji_i have the same form as Jn, Ji_i, whereas Jio, ^oo differ 
from Jio, Jqo by the items containing the nonphysical form factors B^, Bq, Bj. Other 
nonphysical form factors -Bi_4 and Bg do not contribute to these matrix elements. 

The matrix elements Jyx do not satisfy the condition ( |6.61|) . Substituting Jyx in 
eq.( |6.61| ) instead of Jyx, we get: 



A = (1 + 2r7) Jn + Ji_i - 2J2r]J,o - J 



^00 



iB, + By 



(6.64) 
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No. 


Version 


Form factor T\ 


Form factor 


Form factor ^1 


1 


GK-B [0 




T2 


^1 - ^6 


2 


GK-A g 


«l 


Tx 


^2 


^i-^6-2^A 


3 


BH 


102 








^1 ^6 l + 2r?'^ 


4 


D 




^2 + ^A 




5 


CCKP |i 


)5| 


^^+2(1 + .)^ 


-^2 




6 


FFS 


|10( 




-^1 


2(1 + 


The same as for CCKP 



Table 6.1: The set of form factors T\^T2-,Q\ calculated in six versions discussed in the 
section |6.3.3| . 



This function was calculated for the deuteron in ref. |^ and for p-meson in refs. |103| . 
Note that A in (|6.64 ) should not be confused with the perpendicular component of the 
momentum transfer q. 

The substitution in ( |6.62|) of the matrix elements Ja'a from (|6.63|) gives: 



■GK-A 



GK-A 



tGK-a _ 

J- 2 — -r 2 



To 



1a, 

2ri 



(6.65) 



This solution gives T\^T2 without any nonphysical contribution (though still approxi- 
mate), but Q^^~^ contains the nonphysical form factor and, in addition, the function 
A = — (-B5 + B'j\ eg. ( |6.64| ), responsible for the violation of the condition ( |6.61|) . 

The different solutions were compared numerically in ref. ||103|| (for the case of p- 

m 



The results are summarized in table |0 



meson form factor) and analytically in ref. 
Calculations by analytical formulae from this table exactly coincide with the numerical 
calculations ||103| . 

one can see that all the methods give superpositions of the 



From the table 6.1 



physical form factors with the nonphysical ones. Most of them contain the nonphysical 
contribution A and differ from each other by its coefficient. This is so not only for the 
methods considered above, but also for any other method based on the J+ component. 
If A = 0, the form factors do not depend on the prescription, as it was pointed out in 
103|| . But Q\ still contains B^. 

current com- 



(0) 



Note that the terms B^^ contribute to the matrix elements of the J[ 
ponent if at least one of the projections A, A' equals to zero. In ref. |^ it was proposed to 
avoid using the zero projections of A, A', but consider the matrix elements with A, A' = ±1 



(0) 



of the current components J. 
for ijj = (1, 0, 0, —1), one gets: 



and J?) 



{ J^^-*, Jy^^- Taking into account that a;+ = 



ef {p')BTe'M = 0, el'' {p')BTe'M = for A, A' = ±1 
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We see that this method indeed gives the solution for the form factors, in which the 
nonphysical contributions vanish entirely. 

In both methods the complete separation of the physical form factors from nonphys- 
ical ones is realized at the expense of incorporating, besides J+, other components of the 
current. A similar situation takes place with the form factors of a spin 1/2 system P8| . 
These other components can receive contributions not contained in J_|_ (for example from 
the instantaneous part of the fermion propagators, which do not contribute to J_|_ due to 
7^ = 0). The instantaneous contributions (contact terms) have been already taken into 
account in the deuteron electrodisintegration |7^. As we shall see in the following sec- 
tion, these contributions have a very important physical interpretation in terms of meson 
exchange currents in the non-relativistic approach. They can therefore be interpreted also 
as a two-body current. Considering the J+-components only, these contributions would 
be lost. This also shows that the consideration of other components of the current, beside 
J+, is inevitable if one wants to keep all physical contributions, and only these, to the 
form factors. 



6.3.4 The transitions 1+ - 0+ and 1 - 0" 

The pseudo-vector - scalar transition 1+ — 0+ corresponds to the deuteron electrodisin- 
tegration amplitude near threshold (a few MeV in the center-of-mass system of the np 
pair). In this case, the final state is dominated by the singlet S-wave and the amplitude 
does not depend on the direction of the relative np-momentum. The amplitude of the 
vector - pseudoscalar transition 1~ — 0~, corresponds, for example, to the p — tt tran- 
sition, and has the same structure. For definiteness we will discuss here the deuteron 
electrodisintegration. 

In the impulse approximation, LFD was applied to the deuteron electrodisintegration 
amplitude in ref. |l(J5] and, in its explicitly covariant form, in ref. |l(J6|. The corresponding 



amplitude can be written as: 

M2^3 = -^u{K)Yu{K)el{p)Fl^xCSo) , (6.66) 

where F^p describes the vertex c? + 7* — > ^5*0. The function xC^o) is the spin function 
of the final nucleons in the ^So-state we discussed in section |5.3| . The vertex can be 
represented as follows: 

F,p = ^e,,.,g'^P^A , (6.67) 

where P = p + p' , p' = pi + p2 and A is a dimensionless scalar function: A = A(Q^, u). 
The tensor ( |6.671 ) is the only one which has the appropriate properties expected from 
covariance, transversality and parity conservation to describe the vertex d'j* — >■ np{^So)- 
Let us emphasize that the vertex ( |6.67] ) automatically satisfies the transversality condition 
Ff^pqf = 0. In this particular case, the Lorentz invariance leads to gauge invariance for 
the amplitude. This is a consequence of the particular ^Sq final state. The cross-section 
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of inclusive electron-nucleus scattering is usually represented in terms of the structure 



functions Wi and W2 (see, e.g., ||107 



da ( da\ , ,0 , , , , 

(W2 + 2tan2(-)iyi) . (6.68) 




dVt.dE' Xdnj^' '2 



The structure functions W2 and Wi can be expressed through the unique scalar amplitude 
A [HI: 

W^2 = Tr|^|Ap, (6.69) 

1^1 = (1 + uyQ^)W2 . (6.70) 

Here p* is the c.m.-momentum of the final nucleons. All the dynamical information 
about the process is therefore contained in the transition form factor A. We stress that 
this result follows only from the assumption that the final ^Sq state dominates and is 
valid independently of the presence or absence of final state interaction, meson exchange 
currents and isobar configurations. Their contribution only influences the amplitude A. 

The particular amplitude for the vertex dj* — > np(^S'o) in any approximate calcu- 
lation, denoted by F^p, does not coincide, however, with F^p in ( |6.67| ). The tensor F^p 

we can 



depends in that case on u. Following the arguments presented in section ^ 
decompose on the general invariant amplitudes. It has the form: 

Fpp = ^^epf^u^q^p^A + Cp^u-yq"^'^ Bi + Cppy^p'^uj'^ B2 (6.71) 



where = epa/3-y ^'^ q'^ P^- The decomposition ( |6.71|) contains the symmetrical structures 



in front of the functions -83,4,5. Corresponding antisymmetric terms (like V^gp — Vpg^) are 
not independent and can be expressed through the first three contributions. We em- 
phasize that Fpi, depends on u even for the components of the wave function which do 
not depend on uj. In the latter case, u enters through the rules of the graph technique. 
Similarly to the case of the deuteron form factors, the decomposition ( |6.71|) enables us to 



separate the w-independent parts from nonphysical w-dependent ones so that one can un- 
ambiguously extract the physical form factors from the initial tensor F^p. From eq. ( |6.71| ), 
we immediately find: 



2 

m 



This simple formula is the analogue of formulae ( |6.54D - ( |6.56| ) for the deuteron form 
factors. 



Comparison with other approaches 

The transition form factor A given by ( |6.72 ) cannot be reduced to the p = + component 
of the amplitude F^p. For example, in the reference system where go = 0, and g'is parallel 
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to the X-axis, eq.( |6.72| ) is reduced to: 



^ = 7^(^.+ -^+.)- (6-73) 

This formula includes F+y with p = y. 

The contribution of the plus- component is given, as usual, by the contraction F^pU^. 
If F^p would coincide with F^p in eq.( |6.67|) , i.e., would not contain any non-physical 



contributions, we indeed could find the form factor from the plus-component by the 
equation: 

A = -^^^e"^''^u;^q,p,{F,pu;n ■ (6.74) 
However, applying this equation to F^p, we get: 

27T?^ 

A' = -^^^—^e''^''^u;aq.p-y{F^pUjP) = A - B, . (6.75) 

Like in the cases considered above, we get a superposition of the physical form factor A 
and nonphysical one B^. Rewritten in terms of the components, the contraction (|6.75|) 
obtains the form: 

A' ^ =I^F,. . (6.76) 

The difference with eq. (|6.72|) is reduced to the replacement of Fy^ in ( |6.76| ) by the an- 
tisymmetric combination {Fy+ — Fj^y)/2 in ( |6.73| ). This seems very natural since the 
physical amplitude (|6.67| ) is antisymmetric relative to permutation of the indices p, /i, 
and, hence, any symmetrical contribution is a spurious one and has to be excluded (in- 
dependently of the fact that it depends on oj or not). This is another reason to separate 
it. We emphasize that the difference between the amplitudes A and A' is of relativistic 
origin and disappears in the non-relativistic limit. 



6.3.5 The transitions 0" - 1+ and 0+ - 1" 

The pseudo-scalar - pseudo-vector transition 0^ — 1^ corresponds, for example, io ix — Ai 
and to K — Ki transitions. The scalar - vector transition amplitude 0^ — 1^ has the same 
form. 

The general transition amplitude reads: 

{PS{p')\Jp\PV{p, A)) = m{p-q)g^p - ppq^ + F^iqpq^ - g%^)]ei(p) . (6.77) 

It is determined by two form factors Fi and F2. 

In an approximate calculation, this amplitude should incorporate the a;-dependent 
contributions. These are given by: 

{PS{p')\Jp\PV{p,X)) = G^A(P) ' (6-78) 
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where 



(6.79) 



The terms A', C", S2 do not contribute to eq.( |6.4| ), since p^ep{p) = 0. They appear however 
before contraction with e^{p). The term Cp'^qp appears, since an approximate amphtude 
is not gauge invariant. It should be also separated. Its appearance has nothing to do 
with the cu-dependence or LFD. In the Feynman approach its separation is equivalent to 
the construction, by hand, of a gauge invariant amplitude by the following subtraction: 
^ Gpp - {Gp:pqP')qp/q^. 
One can construct other gauge non- invariant terms, proportional to p'^Pp, q^Pp and 
qp.p'p, but they are expressed through the structures included in (|6.79|) . The term VpVp 
with Vp = epa/3'y^'^q^p'^ is also expressed through them. 







Gpp 
(g2)2 









G 



pp 



qPi^P 
(uj-p) 



I / 2 2 2^ 

+ imps -m^v-Q ^ 



pv 



2{uj-py 



(6.80) 



^ps 



r, U^UJP , 2 2 
(ini^T-^ + \"^pv -1T^ps + Q 



fqi" yOof^pP + ujPp'>' 
+ q - 



[uj-p) 



(uj-p) 



''2{uj-py 

+ qf'qf - q^g^'P 



2 qPu^' 



(6.81) 



Here rups and rupy are the masses of pseudo-scalar and pseudo-vector mesons. 
Comparison with other approaches 



In ref. [[I08|| the form factors are found from the component p = +. Contracting ( |6.77|) 
with Up, we get: 

{PS{p')\Jp\PV{p, X))u;P = [Fi((p-g)^^ - {u;-p)qn - F^q^ uP]el{p) . (6.82) 



Contracting (|6.78| ) with ujp, we reproduce the equation (|6.82|) with new form factors F[, 
given by: 

F[ = Fi, F2 = F2 — Bi . 

Hence, the component provides Fi without any spurious contributions, while the form 
factor F2 is determined by: 



p' - p R - PP 
-^2 - -^2 - i5i - -^7 r 



2 ^ u 

m^„, pp 



pv 



uj-p 



(6.83) 



It includes the spurious contribution Bi 
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Chapter 7 

Electromagnetic observables 



The formalism and methods presented above are apphed to different physical systems and 
processes: pion, nucleon and deuteron electromagnetic form factors, p — vr and ed — enp 
transitions. With standard assumptions and approximations, we easily reproduce results 
given in the literature. Some of these results, however, like the asymptotic behavior of 
the pion form factor and meson-exchange current contributions in ed — enp amplitude, 
obtain a new and clear physical interpretation in our approach, since they directly reflect 
the presence of extra relativistic components in the corresponding wave functions. 

7.1 Electromagnetic form factors in the Wick- 
Cutkosky model 



The Wick-Cutkosky model allows to get the explicit form of both the Bethe-Salpeter 
functions and the light-front wave functions for states with angular momentum J = 
and J = 1. This gives us the opportunity to compare the electromagnetic form factors 
calculated by two methods: by the standard one, with the use of the Feynman rules and 
through the Bethe-Salpeter amplitude, and by means of the light-front graph technique 
through the light-front wave function. This comparison is especially non-trivial in the 
J = 1 case, where the physical form factors should be calculated after separating out the 
non-physical (u;-dependent) amplitude from the physical one. 

7.1.1 Spin 

For a spin-0 system consisting of spinless constituents, the form factor is expressed through 
the light-front wave function by eq. (|6.14 ). Its expression through the Bethe-Salpeter 



function corresponding to the Feynman graph of fig. |2l| reads: 



{p + p'),F{t) = I J (p + p'-2k),<l>i^p-k,pM^p' - k,p'){m' - k')-0^,. (7.1) 
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P k p' 

Fig. 21. Feynman electromagnetic vertex of a two-body system. 



|t|/m2 


Fit) 


BS 


LF 


NR 





l.OOE + 00 


9.67E - 01 


l.OOE + 00 


0.01 


5.14E-01 


A.98E - 01 


5.17E-01 


1 


6.02E - 04 


5.75E - 04 


6.23E - 05 


25 


l.OlE-06 


9.76E - 07 


1.05E - 06 


100 


6A0E - 08 


Q.20E - 08 


6.55E - 08 


225 


1.28^-08 


1.24E - 08 


1.29E - 08 



Table 7.1: Electromagnetic form factor of a system with spin J = in the Wick-Cutkosky 
model for a = 0.08. Columns BS, LF and NR represent the results of the calcula- 
tion using the Bethe-Salpeter amplitude ( [7.2|) , the light-front wave function ( p.64| ) and 
non-relativistic wave function (|3.62| ) respectively. 



The Bethe-Salpeter function is given by eq.( |3.66| ). Substituting here the spectral 



function ( |3.69| ), one can find the explicit form of the Bethe-Salpeter amplitude: 



'IC 



m 



2 



m 



rn? — {-p — — iO 



(7.2) 



where c = 2^v7rm/t^ with k = Jmleb 



ma/2. 

The analytical (in the asymptotical region 
factors calculated through the light-front wave function, eg. ( |6.14| ), and through the Bethe 
Salpeter amplitude, eq.( f7Al) , was carried out in ref. |85| 
these two approaches are compared in table |7?1| . 

Both approaches give the same asymptotical behavior of the form factors at \t\ ^ 



and numerical comparison of the form 
. ( |6.14| ) , and through the Bethe- 
The form factors calculated by 



2. 



Fit) 



16a^m^ 
t2 




(7.3) 
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where a = j (IGvrm^), g is the couphng constant in the Wick-Cutkosky modeL 

In contrast to ref. the LF form factor in table ^]T] has not been renormahzed at 
t = 0. In the light-front approach, the form factor at low t, as well as the normalization, 
is smaller than the Bethe-Salpeter one by 4Q;/37r. For the normalization, the discrepancy 
has its origin in the contribution of the derivative of potential to the norm operator, as 
indicated in eq.( |3.5CI|) . For the form factor, this corresponds to omitted contributions 
indicated in fig. ^ (the so-called recoil contributions). If corrections to the form factor 
and to the norm (which coincide with each other at = 0) are included, then the 
renormalization of the light-front result, as done in ref. is not required. 



At high t, the results obtained in LFD and in the Bethe-Salpeter approach are very 
similar. Non-relativistic calculation is also close to them, when (a/27r) log(|t|/m^) <^ 1. 
This is a peculiarity of the Wick-Cutkosky model. In the asymptotical region where 
(a/27r) log(|t|/m^) > 1 LFD and Bethe-Salpeter form factors are still close to each other, 
while they differ from non-relativistic one. 

The asymptotics of the non-physical form factor Bi, which appears in the decom- 
position ( |6.9| ) and is given by eq.( |6.15| ), has the form: 



(7.4) 



In the intermediate region, Bi exceeds F : Bi{t) / F{t) ^ 2 while it is negligible in the 
non-relativistic region. According to |^ , Bi becomes comparable with the physical form 
factor F{t) at \t\ ~ rn? . Therefore its separation is important. In the spin-0 case it 
is separated automatically, if one proceeds from the p = + component of the current. 
However, this is not so for the spin-1 case. 



7.1.2 Spin 1 

For a system with angular momentum J = 1 made out from two spinless constituents, the 
light-front matrix elements J^'^ of the current are obtained from the spin-0 case equation 
( |6.12| ) with the replacement of the spin-0 wave functions if)' , ip by the spin-1 wave functions 



Tp'^^ipy. In four-dimensional notations, Tp^j, has the form: 

<p\hM,P,i^r)=e^^\p)r, r = ^1^^^^^^ + ^2— ■ (7.5) 
^ 2 ixJ-p 



The explicit form of the components ipi^2 in ( [i".5|) can be found from the comparison of 
( [7.5| ) with ( |3.65| ) in the reference frame where ki + k2 = 0. Substituting J^'^ in eqs.( |6.54|) - 
( |6.5(j| ), in order to separate the nonphysical contributions, we find three form factors 

•^1) -^2) ^1- 

On the other hand, in the Bethe-Salpeter approach, the electromagnetic vertex is 
directly related to the physical form factors by the decomposition ( |6.40 ) and can be 



obtained from (|7.1|) by the replacement of the spin-0 Bethe-Salpeter amplitude by the 
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Fig. 22. Amplitudes corresponding to admixtures of the three-body sector in the state 
vector. These amphtudes correspond to retardation effects. 
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Fig. 23. The form factor Qi in the Wick-Cutkosky modeL The sohd hne represents the cal- 
culation with the Bethe-Salpeter function (|7.6|) and by eq. (|6.56|) with the light-front wave 
function ( ^.65| ). Dashed line illustrates the form factor incorporating the non-physical 
contribution by eq. (|6.65| ) (solution A from ref. 



H) 



The form factor Qi calculated in the two approaches in ref. pH] is shown by the solid 



line in fig. The two calculations give the curves indistinguishable from each other. In 
order to show the importance of a proper separation of the non-physical contribution we 

(dashed 
It 



show also in fig. ^ the form factor corresponding to the solution A from ref. 
line). It is given by eq.(|6.65|) and contains the non-physical contributions A and 
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considerably differs from the correct form factor. According to ref. |]86|, this difference in 
the model is mainly due to A, whereas the Bq contribution is negligible. 

In the asymptotical region, these Bethe-Salpeter and light-front form factors are 
given by the following analytical expressions fSl 




l + ^log(-^) , (7.7) 
i + ^ (7.8) 

(7.9) 

We emphasize that, generally speaking, using two-body Bethe-Salpeter function 
does not correspond to two-body Fock component, but incorporates implicitly many-body 
sectors of the state vector. 



7.2 Applications to the quark model 

We shall investigate in this section a few simple examples on the application to the quark 
model of the formalism developed in this review. We shall discuss in some details the 
pion wave function and form factors, and give some results for other observables. 

7.2.1 The pion form factor 

Pion wave function 

The spin structure of the pion wave function as a system consisting of quark and antiquark 
in the J'" = 0~ state is identical to the np wave function in J'" = state, eq. ( p.33| ). The 
negative parity is generated automatically by opposite internal parities of the quark and 
antiquark and therefore does not change the spin structure. The pion wave function has 
thus the form: 



.1 . ^ 

A1— + A2 

m uj-p 



-f5v{h) = u{k2)0v{h) , (7.10) 



where O = {Ai/m + A2U}/uj-p)'y^, and m is the quark mass. For simplicity, we do not take 
into account isospin. 

The representation of this wave function in terms of the variables k and n is almost 
identical to ( |5.34D : 

1 . / ia-\n X k] \ /^, , , x 

^ = 7f ^2 I 91 + -92 1 w, , (7.11) 

with the following relations between the invariant functions: 

Ai = -^^(9i + j92), A2 = j92. (7.12) 
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The normalization condition is a particular case of ( p.23| ): 
where 



(7.13) 



A1A2 



D 



1 d^k 



(fR I dx 



(27r)3 (1 - x)2£fc, (27r)3 {2tt)^ 2x{l - x) ' 

Substituting in eg. ( [7. 13] ) the wave function represented by ( [7. 101 ) and ( |7.11| ), we get: 



(7.14) 



No 



m^x(l — x) 



D = j[gl + {l~z')gl]D , (7.15) 



where z = cos{kn). Obviously, in the case where the interaction potential between quarks 
is independent of M, as often assumed in phenomenological approaches, N2 should be 
equal to one. 



Pion form factor 

The diagrammatical representation of the pion electromagnetic current in the impulse 
approximation almost coincides with fig. |1^. For the interaction of a virtual photon with 
a quark, the horizontal line has to be replaced by the double line of antiquark. This 
current has the form: 



j;-? = / Tr[-0'{k'2 + m)j,{k2 + m)0{m - ^OIt^^TT^ 



(7.16) 



where O depends on the initial momenta, and O' depends on the final momenta. The 
quark current jp is taken as: 



(7.17) 



In this equation /i and /2 are the quark form factors. The minus sign in Tr[—0' ■ ■ ■] 
is from the fermion loop. We remind that both lines in the loop, single and double, 
corresponding to the quark and antiquark, are in the same direction. But since, according 
to the rules of the graph technique, they are followed in the same and opposite orientations 
respectively, one makes a loop when going through them. 

The pion form factor for the 7g interaction is given by: 



Tr[—0'{k'2 + ^)('-^'j)(^2 + m)0{m — ki)]] 



2uj-p(l — X 

The pion form factor for the 'jq interaction is found similarly: 

1 



■D 



F^%Q^) = / Tr[-0{m - k2){uj-j){m - k'^)0'{m + ki)] 



2uj-p{l — x] 



D 



(7.18) 



(7.19) 
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The momenta k2, k'2 in ( |7.16| ) are the quark momenta, whereas they are the antiquark 



momenta in (|7.19|) . 

After a trivial generahzation for the case of different masses of the quark and anti- 
quark, these formulas can be applied to the i^-meson form factors. The final form factor 
is expressed through the momentum transfer A (with = A^) and the integration 
variables R±,x with the use of the invariants indicated in appendix 0. 

Let us for a moment put A2 = 0. We thus have: 



x(l — x)m? 



(7.20) 



With Ai expressed through gi by ( [7.12| ) the form factor ( [7.20| ) exactly coincides 



with the expressions given in ref. ||109|] . We emphasize that the calculation of (|7.20|) , and 
more generally of any form factor or transition amplitude is, due to the covariance of our 
approach, a simple routine for analytical computer calculations. 

When A2 is different from zero, but for pointlike quarks, i.e. with /i = 1, /2 = in 



the current (|7.17|), the calculation of (|7. 18|) gives 



F^5(Q2 



x(l 



AiA'^ + 2(AiA'2 + A'l^s) + 4x(l - x)A2A'^ 



(7.21) 



where A[ 2 = ^1,2 (-R± ~ xA,x). We shall use this expression in the following to evaluate 
the asymptotical behaviour of the pion form factor. 

Note that for the interaction of the virtual photon with the antiquark, in eq. ( [7.19|) , 
we get: F'^^ = —F^''. Hence, the full form factor of a system consisting of identical qq is 
zero: F^'^{t) + F'^^ = 0. We get opposite sign of F^^ automatically, without introducing 
negative charge of antiquark. This negative charge is automatically taken into account in 
this formalism. It can be obtained from the antiquark contribution, multiplying it by a 
sign factor depending on the process under consideration. 

For the p — it transition, the quark and antiquark contributions have the same sign. 
The transition amplitude can be obtained in this case from the vr elastic amplitude 
Jp by replacing in eqs.( |7.1(^ ) the initial pion wave function O by the wave function of 



the vector meson 
function (15.11). 



The form of this wave function coincides with the deuteron wave 



Asymptotical behaviour of the pion form factor 

It is well known that QCD provides a 1/Q^ asymptotical behaviour for the pion form factor 
4|. As it was shown above, the pion wave function contains in general two components. 
The first one leads to the usual non-relativistic wave function, while the second one is 
of purely relativistic origin and w-dependent. We shall show in this section that the 
asymptotical l/Q"^ behaviour of the pion form factor is entirely determined by this last 
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a;-dependent component of the pion wave function ||110|| . This piece is essential in order 



to explain the difference with the form factor of the J = state in the Wick-Cutkosky 
model, which asymptotical behaviour is l/Q^. 

The asymptotical behaviour of the form factor can be related to the asymptotical 
behaviour of the wave function. We shall calculate the latter in a qq model with a one- 
gluon exchange kernel. The equation for the wave function corresponding to the diagram 
indicated in fig. |TT] has the form: 

^(A;i,fc2,p,a;r) = -— ^ j ^/') + ur' - k[ - k',)2{u-p)dr' 

xu{k2hf,{k'2 + m)e{uj-k'^)5{k'^ - m^)d%0\m - k[)e{uj-k[)5{k'^ - m^)d^k[Yv{ki) 
xK:{k[,k'2,ujT';ki,k2,ujT). (7.22) 

Here M is the pion mass. We use the gluon propagator in the Feynman gauge and /C is 
the scalar part of the gluon propagator (after separation of —g^u)- It coincides with the 
scalar t-channel exchange amplitude, given by eq. (|2.67|) (where one should put fi = 0). 



The matrix O' is the "inner" part of the pion wave function as defined in eq. (|7.1CI|) . Prime 
means that it depends on the momenta in the loop. For simplicity, we omit here any color 
degrees of freedom which are irrelevant for our consideration. 

To find the asymptotical behavior of the wave function, we shall use the iterative 
procedure already explained in chapter ^ to find the deuteron wave function. We suppose 
that O' is concentrated in a finite region of the quark momenta and the integral ( |7.22| ) 
is dominated by this domain. These momenta are negligible relative to the asymptotical 
ones we are interested in. This means that we can evaluate all the factors except for O' 
and the delta- functions at zero relative quark momentum. This corresponds to -R'j^ = 
and x[ = x'2 = 1/2. In that case, the quark four- momenta are equal to each other: 

K = K = {p + UJTo)/2, 

where tq = (4m^ — M"^) / {2uj-p) . We thus get the following replacement in (|7.22 

^ r^^^) {p + ur' - k[ - k'2)2{uj-p)dr' 



{2^f 

xd{uj-k'2)6{k'i - m^)d% O' d{uj-k[)5{k'^ - m^)d% Oq, (7.23) 

where Oq is the integrated value of the wave function. 

We substitute in ( |7.22| ) the one-component wave function obtained from (|7.10|) with 



A2 = 0. Together with the replacement (|7.23|) , this gives in ( [7.22|) O' = Ao'yr^/m, where 



Aq is a constant. The second component A2 of the wave function is then automatically 
generated by eq.( [7.22|) . 



We thus find for the pion wave function: 



2 A 

^ = ^u{k2){p + uTo-Amh,v{k,)K: . (7.24) 

S — 
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Substituting here p from the relation p + oor = ki + k2 and using the Dirac equation, 
we finally obtain: 



2A 



(7.25) 



S-M2 

The pion wave function ( |7.25| ) has the form of eq.( |7.1CI| ) with the following two components: 



s — 



(7.26) 



_j_ ^2 

with s = — Pi The kernel /C given in ( |2.67| ) can be rewritten in terms of the variables 

Ri ,x. For X < x' it reads: 



}C{R'^,x';R^,x,M'^) 



+ m' + -(R^-^R'^ 

X \ X / X \ X 



x'(l — x') 



(7.27) 



with /i = 0. For x > x' K{R'j_, x'\ R±, x, M^) is obtained from (|7.27|) by the replacement: 
X ^ x', R± R'j^. As indicated above, one should then put R'^ = 0,x' = 1/2. 

We will see below that this wave function enters in the form factor at x = 1/2. Then 
s ^ AR]_. As follows from ( |7^ ), the kernel /C at R'^^ = 0, x' = 1/2, x = 1/2 obtains the 
simple form: /C = / R\. The components of the wave function which dominate in the 
asymptotical region are then given by: 



^, = ^^, A2 
i?1 



R^. ■ 



(7.28) 



In this region, the component A2 decreases more slowly than Ai. 

The form factor is now easily calculated through the components Ai,A2 according 
to eq. ( [7.21| ) . The integral for the form factor is dominated by the region where one of the 
two wave functions (initial or final ones) takes its non-relativistic value, whereas the other 
wave function corresponds to the asymptotical one. This is equivalent to the standard 



calculation discussed in the literature in the usual formulation of LFD Hence, 
similarly to eq. (|7.22|) , we can replace A[ , A'2 in ( [7.21j ) by 



A[ = Aoj,6^'\R^)6{x - 1/2), A'^ = 0. 



and take the asymptotical values (|7.28 ) for Ai and A2. We thus find 



Ar 



(2Ai + A2 



A, 



(7.29) 
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The substitution of A2 from ([7:28| ) in (^^29]) with R^ = A/2 and = gives: 



From ( [7 .291) and (|7.30|) we see that the asymptotical behaviour oc of the pion 



form factor is indeed determined by the extra component A2 of the pion wave function, 
related to the cj-dependent spin structure. This originates from the slower decrease of 
A2 oc l/R± in comparison to Ai oc l/R± (see eq.( |7.28| )). The enhancement of A2 is 
determined by the factor t{uj-p) = {s — M^)/2 ?a 2i?5_ iii ( [7.25| ). This extra factor 
originates from the part ujt in the conservation law p + ut = ki + k2 used to derive 
eq. (|7.25|) and therefore accompanies namely the a;-dependent spin structure. 



7.2.2 The nucleon form factors 

As we have already seen in the case of the two-body systems, it is convenient to parametri- 
ze the covariant light-front wave function in terms of the following three sets of variables: 
(i) The four-momenta ki of each constituent. In these variables, the wave function 
of any three-quark system has the following form: 

2p = ilj{ki,k2,k3,p,LJT]ai,a2,cr3,a) , (7.31) 

where a is the nucleon spin projection and o"i.2,3 are the quark spin indices, with the 
following conservation law |^: 

ki + k2 + k^=p + ujT . (7.32) 

(m) The three- vector variables qi which are identical, in the c.m.s., to the constituent 
momenta. They are constructed similarly to the two-body case, (|3.11|) : 

gl = L-\V% , (7.33) 

where V = p + ujt = ki + k2 + k^ and L^^{V) is the Lorentz boost operator defined in 
( p. lip . The sum of gives zero, as it should be for the relative momenta: qi+q2 + qz = 0. 

{iii) The four- vectors: Ri = ki — XiP, where Xi = uo-ki/uo-p. They can be represented 
as Ri = {Rio, Ri±, Ri\\), where Ri_L-u = 0, Ri\\ is parallel to uj. We thus have the following 
relations: 

i?l_L -I- i?2± + -RsX = , 
Xi + X2 + X'i = I . 

We do not pretend here to a calculation of the form factors in a realistic model for 
the nucleon. In order to illustrate our approach and estimate the effect of eliminating the 
unphysical Bi contribution from the physical form factors, we consider a simple model in 
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which two neutral "quarks" labeled 2 and 3 are coupled to zero total angular momentum, 
so that the spin of the nucleon is determined by the charged "quark" labeled 1. We put 
the masses of all the quarks equal to each other. The corresponding wave function reads: 

i){ki,k2,k3,p,UT;a,ai,a2,a3) = [M''2(A;2)75f/cM''3(A;3)][M''i(A;i)M'"(p)]^o(^i, ^2, ^3,^,^^^) , 

(7.34) 

where Uc = 7270 is the charge conjugation matrix. The spinor structure of ( [7.34| ) coincides 
with one of the invariants given in ||111|| . In order to show the method more clearly, we 
do not symmetrize the wave function ( [7.34| ) and do not introduce isospin. So, we have 
only one (charged) nucleon - "proton" . In eq. ( [7.34| ) ipQ is a scalar function. 




Fig. 24. Nucleon electromagnetic amplitude in the three-quark model. 

We consider the charged quark 1 as a structureless particle coupled to the electro- 
magnetic current: jp = u{k[)'y pu{ki) . The electromagnetic nucleon vertex ( |6.25| ) is shown 



in fig. |2J and has the form: 



Jp = u'TpU = j[u{p'){k[ + m)'^{k+m)u{p)]Tr[{k2 + m^ . (7.35) 

Here ipQ and ip^ depend on the initial and final momenta respectively. Details of the 
calculation can be found in ref. PB|. 

At zero momentum transfer, the difference (j'^^(O) — Gm(0) obtains the simple form: 



6-^,(0) -G'm(0) = 25i(0) 



Here {k ^) is the average quark momentum. The difference between G'j^ and Gm is of 
course of relativistic origin. 

We indicate in fig. ^ the ratio R = Gm/G'^.j of the physical form factors Gm 
obtained with our formalism and the form factor which would be obtained with the 
usual procedure (using the J+ component of the current), as discussed in section 3.3. 
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Both form factors are calculated with a simple harmonic oscillator radial wave function 
|88| . We recall that the electric form factor Ge is the same in both approaches. 

The difference between Gm and G"^^ is of relativistic origin. It is already sizable at 
= for the magnetic moment (enhancement of about 15%). The two form factors 
largely differ at high momentum transfer, where corrections of the order of 50% are to 
be expected. The physical form factor Gm-, i-e. the form factor free from any ambiguity 
from the position of the light front, is always larger than in this model. According to 
fig. ^ the magnetic radius is larger for Gm than for G^. 



7.3 Application to the nucleon-nucleon systems 

7.3.1 Light-front dynamics and meson-exchange currents 

Non-relativistic phenomenology 

The qualitative, and to a large extent also quantitative, understanding of electromagnetic 
observables in the few-body systems in terms of meson-exchange currents is now well 
established [^. This is particularly transparent for isovector magnetic transitions like 
the deuteron electrodisintegration cross-section near threshold. These two-body currents 
appear in that case already at order 1/m, while they are of higher order for charge and 
isoscalar magnetic transitions. Last but not least, these currents are also required by the 
low-energy soft-pion theorems. 

In order to extend these non-relativistic analyses to a wide range of momenta, let 
say for > M^, a relativistic framework is necessary. In order to be fully consistent, the 
relativistic analyses have to give, as a limiting case, the well-known results in terms of 
meson-exchange currents. While genuine two-body currents like the mesonic contribution 
indicated in fig. ^ has to be considered explicitly in both LFD and non-relativistic 
formalisms, some of them may appear differently depending on the underlying framework. 
Physical amplitudes should however be the same. 

This is the case for the so-called pair term which, in the non-relativistic approach, is 
indicated on fig. (we assume here that the nNN vertex is of the pseudo-scalar type, we 
shall come back to this point at the end of this section). This contribution cannot survive 
on the light front because of the absence of vacuum fluctuation diagrams, together with 
the particular constraint u-q = 0. A carefuU analyzis of the deuteron electrodisintegration 
amplitude near threshold in LFD is therefore of particular interest. 



Relativistic impulse approximation 

The vertex d + j* —>■ np{^So) is represented graphically in fig. 2S. It corresponds to the 
following amplitude: 



Tr[75(fc2 + m)r)f (/C2 + m)(j)^{ki - m)] (Pki 



(1 — uj-ki/ uj-pY'' 



(7.36) 
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Fig. 26. Two-body meson exchange current: the mesonic current. 
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Fig. 28. The electromagnetic vertex 7*0? np{^SQ 



where c = i2^^^^'^^l^^^/m, 0^ is the deuteron wave function (^]2|), v^/ is the final state 
wave function and Tp is the nucleon electromagnetic vertex ( |6.16| ). The superscript V 
refers to the isovector part of the form factors which is the only one which contributes 
here. In the plane wave approximation, we simply have: 



2mTr[75(/c/ + m)rp {k2 + m)(/)^(A;i - m)] . (7.37) 

The expression of the electromagnetic current used here involves the Pauli and Dirac 
form factors Fi and F2. They differ off-energy-shell from the Sachs form factors used in 
ref. [ p,06|| . The present one for its Fi part is consistent with the minimal substitution in 
the Dirac equation. Sachs form factors may be used, but extra terms have then to be 
introduced in the current to get the exact result. 

The amplitude A, determining the cross section ( |6.68|) is found by the formula ( |6.72|) . 



In the plane wave approximation and for zero c.m.-energy in the final state, the amplitude 
A calculated with the deuteron wave function (h.2\) has the form: 



5/293/2 

A = i2Gl^2 - Glcp,) (7.38) 



with k = A/2. We introduce in ( |7.38| ) the charge and magnetic nucleon form factors ac- 



cording to ( 16.171 ). The deuteron components v^i, 3,4,6 do not contribute to ( |7.38| ). Express- 



ing ip2, through /i,/2,/5 by eqs.( |C.4|) , ( |C.7| ) (and neglecting the small contributions 
from fs and f^) we get: 

A = 2m'/' iGl[usiA/2) + -^uniA/2)] - Gy^f,(A/2)\ , (7.39) 
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where A = \fQ^. 

Neglecting /s in eg . (17.391) , we recover the usual expression in the plane wave ap- 
proximation. 



Meson exchange currents and 1/m expansion 



In order to make the connection with the non-relativistic phenomenology, it is instructive 
to calculate the /s component in leading order in a 1/m expansion. Following the iterative 
procedure used in refs. and recalled in section |^ , the /s component has the 

following form in momentum space and in 1/m order |7^: 



/sffc X n\ilk 



X 



8 V37r2m2 J (fc 2 + ^2) ^2 + _ ^ /)2 

[{k~k') X n]j . 



-^us{k')5ij - ^(S^ - 5ij)uD{k' 



(7.40) 



In the limit where the momentum k' is restricted by the integration domain to low mo- 
menta, one can take (A;^ — k''^)/{k'^ — k^) ~ 1. We shall come back to this approximation 
later on. 

The coordinate space integral for the /s component has the form: 



h{k) 



3g 



Inn 



Sirm^ Jo 



exp(— /ir) 



(fir + l)ji{kr) 



dr 



(7.41) 



where u and w are the S and D state deuteron wave functions in coordinate space. The 
extra component (72 of the final state wave function can be calculated simply from isospin 
arguments, leading to the replacement —?>g'^^^ — > — (3 -|- l)g1^^ = —Agl^^ in the am- 
phtude. 

The amplitude A has thus the following form, in coordinate space: 

,rA, 



A 



Gl9l 



+ 



NN 



mA 



exp(— /ir) 



ifxr + 1) 



1 rA, 

—=w(r) joi — 
^ V JJ2K 2 ' 



rdr 



u{r) H — -=w{r) 
V2 



(7.42) 



The structure of the second term is that of a two-body contribution of order g'^/An. 
It is similar to the one given by the pair term in the non-relativistic MEC analyses 
P, |112|, |113[. However, its strength is smaller by a factor of 2. 



Contribution from the contact term 

The one to one correspondence with the non-relativistic phenomenology in terms of MEC 
is indeed recovered if one collects all the contributions of order glp^j^/in to the light- 
front amplitude. The missing contribution originates in this particular example from the 
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contact term which appears once we have singled out the one-pion exchange mechanism 
in generating the f^{g2) component of the deuteron (^5*0) wave function. It is indicated 
in fig. |2^. There are four diagrams corresponding to the contact terms. Two of them, 



29|(a) and (b), correspond to the contact interaction at the left of the photon line. Two 
other diagrams, p5|(c) and (d), correspond to the right contact term with the two possible 
time orderings for the lower vertex. Each diagram can be calculated in terms of the 



kinematical invariants given in section U.2 of appendix ^ In leading order in 1/m, the 



contact term contribution is exactly equal to the contribution from the f^{g2) components 
and therefore the total amplitude of the deuteron electrodisintegration near threshold is 
strictly equivalent in LFD and in the non-relativistic phenomenology including meson- 
exchange currents (more precisely the pair term ) [|72| . 

This analysis can even be extended to the whole domain of momenta in eq. ( |7.40 ), 



i.e. without neglecting the factor [k"^ — k''^)/{k'^ + k^) in the integrand. This term gives 
in fact a wrong asymptotical behaviour of the wave function. As shown above in sect. 
4.3.2| and in ref . |Q , incorporation of the complete contact interaction which arises with 



the other components of the kernel (NN potential) like a and ui exchange and of a higher 
order iterations, gives a correction which turns the factor (/c^ — + k^) into 1. 

For the deuteron electrodisintegration, the dominant component (f^ of the deuteron 
wave function does not contribute to the amplitude corresponding to the contact term. 
For the elastic ed scattering with the representation ( |6.16| ) for the nucleon electromagnetic 
vertex Tp, the contact terms contribute to the deuteron electromagnetic form factor Qi 
only and do not contribute to J-'i and 

The relevance of vr exchange. Numerical results 

The 1/m analysis of the relativistic contributions to the deuteron electrodisintegration 
amplitude is based on the vr-exchange interaction with the pseudo-scalar nNN vertex. It 
is also well known that, in leading 1/m order, the pair contribution associated with the 
pseudo- vector ttNN coupling is zero. After minimal substitution in the ttNN vertex, gen- 
erating a genuine NNirj current, the deuteron electrodisintegration amplitude is however 
equivalent to the one given in the pseudo-scalar representation. In LFD, the equivalence 
between the two representations is realized in the following way. In the pseudo-vector 
representation with vr-exchange only, the /s component in the deuteron wave function is 
strictly zero in leading order. This is due to the off-shell condition at the nNN vertex. 
The contribution from the light-front contact term is also of higher 1/m order for this 
representation. On the other hand, the current originating from the direct coupling of 
the photon to the pseudo- vector irNN vertex has its analogue in LFD, providing the 
equivalence between relativistic and non-relativistic formulations in leading 1/m order. 

The accuracy of the 1/m approximation and the contribution of the pion exchange 
relative to all mesons are shown in fig. pOl. The solid curve in fig. is the ratio of /s 



for the pion exchange only, calculated in the 1/m approximation by eq. ( |7.4CI| ), to the /s. 



also for the pion exchange only, but calculated beyond the 1/m approximation. It shows 
the accuracy of the 1/m approximation, which is very good in the non-relativistic region 
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(a) 



(b) 





k (GeV/c) 

Fig. 30. The solid curve is the ratio of /5's in 1/m approximation and beyond, both for the 
pion exchange only. The dashed curve is the ratio of /s for pion only to /s incorporating 
all the meson exchanges, both beyond 1/m approximation. 
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of /c, but changes the result by a large factor for higher k. The dashed curve is the ratio 
of /s for pion only to /s incorporating all the meson exchanges, both beyond the 1/m 
approximation. The pion exchange contributes about 60 per cent in the region k ^ OA 
GeV/c, where /s dominates. 



The cross section of deuteron electrodisintegration is shown in fig.PT|. As it was 
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Fig. 31. The cross section „/ of deuteron electrodisintegration cross-section in 

cm^/MeV/sr for 9^ = 155°. The dotted curve is the non-relativistic impulse approxima- 
tion with the Bonn deuteron wave function and with the Paris scattering state wave 
function [Q. The dashed curve includes the dominant contribution to meson exchange 
currents (the standard pair term), calculated with vr-exchange only. The dot-dashed 
curve is the relativistic impulse approximation with the light-front deuteron wave func- 
tion (/i,/2,/5) and the scattering state one {gi,g2) and includes the contact term. The 
extra components /s, g2 and the contact term have been calculated with pion exchange 
only. It does not include the product /s x g2 and extra components of the wave functions 
in the calculation of the contact term. The solid curve differs from the dot-dashed one by 
the contribution of all mesons both in /s, g2 and in the contact term. 

explained above, there is equivalence, in 1/m order, between the relativistic impulse 
approximation (including all components of the wave functions, together with the con- 
tact interaction), and non-relativistic calculations including the dominant contribution 
to meson-exchange current (the so-called pair term). The comparison between meson 



exchange contributions (pair term) and the light-front approach is shown in fig. 31. The 
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dotted line is the non-relativistic impulse approximation without meson exchange cur- 
rents. The dashed line is the sum of the non-relativistic impulse approximation and of 
the contribution from the pair term for vr-exchange only. The dot-dashed line is the rela- 
tivistic impulse approximation (with /s and together with the contribution from the 
contact term. The relativistic components /s, g2 and the contact term are calculated with 
vr-exchange only and with the pseudo-scalar ttNN vertex, according to the iterative pro- 
cedure indicated in sect.[5.4|. To be consistent with this iterative procedure, we omit in the 



amplitude terms of higher order than g^j^^- The solid curve differs from the dot-dashed 
one by the contribution of all mesons both in f5,g2 and in the contact term. For the 
nucleon electromagnetic form factors the dipole fit ||114|| was used. 



The deviation between the dotted and dashed lines shows the well known influence of 
meson-exchange currents due to the pion pair term. The deviation between the dashed and 
dot-dashed lines indicates that the influence of relativistic effects to both the deuteron and 
the scattering state wave function, and to the electromagnetic operator (contact term) 
are important. The comparison between dash-dotted and solid curves shows that the 
contribution of other mesons, in addition to the pion, considerably influence the cross 
section at large Q^. 



The relativistic calculations in fig. 31 are higher than non-relativistics ones. This 



difference is related to a deviation of relativistic scattering state wave function gi from non- 
relativistic one (see fig. |16]). As seen in fig. [16|, the non-relativistic scattering state wave 



function in momentum space, in contrast to the non-relativistic one, changes sign. On the 
other hand, one can expect that the recoil-type diagram, indicated on fig. partially 
cancels the difference between both calculations. These calculations should therefore be 
considered as an estimate of relativistic effects and are not compared with experimental 
data. 

The difference between the non-relativistic impulse approximation (dotted curve) 
and the relativistic impulse approximation, calculated with the same wave functions (ne- 



glecting /s and (72)5 is negligible. This was already noticed in ref. ||105|, |106 



7.3.2 The elastic deuteron form factors 

The deuteron form factors can be calculated according to eqs.( |6.54j )-( |6.56| ) which express 
them through the tensor Jj^". This tensor is expressed through the deuteron wave func- 
tion by means of the rules of the graph technique given in chapter 0. In the impulse 
approximation this expression is given by eq. ( |6.49| ). It includes the full deuteron wave 
function 0^, with its six components. We shall take into account below the dominant 
ones, /i, /2 and f^. The contribution of the contact term to the tensor J^^ is obtained by 



calculating the amplitude corresponding to the graph of fig. |29|. All relevant kinematical 
invariants are given in appendix 0. 

Like in the case of the electrodisintegration cross-section, the component /s of the 
deuteron wave function was calculated by an iterative procedure, keeping the first degree 
of the relativistic kernel (of the order of g"^). Relativistic contributions of order g^, like 
terms in /|, have been removed from the calculation of the form factor. 
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Fig. 32. The structure function ^(Q^) for the deuteron. The dashed hne is non-relativistic 
impulse approximation calculated with the Bonn wave function. The dot-dashed curve is 
the relativistic impulse approximation with relativistic deuteron components /i and f2- 
The solid curve is the same as the dot-dashed one, but including the /s component (in 
first degree only). The dotted curve incorporates, in addition, the contact term. The 
nucleon form factor were taken in the dipole form ||114 . 
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The influence of relativistic effects on the structure function A(Q^) is shown in 
fig. 0. The dashed curve corresponds to the non-relativistic impulse approximation 
115| , |116|| with the S- and D-waves of the Bonn-QA wave function [0. The dot-dashed 



line is calculated in the light-front formalism with relativistic deuteron components /i 
and /2 (solid lines in figs. |T^ and |T^), but without /s. The solid curve is similar to the 
dot-dashed one, but includes the /s component (in first degree only). The dotted curve 
incorporates, in addition, the contact term (where the contributions of order are taken 
into account only). 




Figure BBI the B{Q'^) deuteron form factor. One can see that the function B{Q^) 
is more sensitive to the different approximations and contributions than A^Q"^). The 
comparison between the dot-dashed and solid curves shows that the influence of the extra 
component /s of the deuteron wave function is important. The minimum of -B(Q^) in this 
figure is the consequence of the /s component. 

Figure ^ shows the influence of different approximations to T20. The effect of 
incorporating the component /s is qualitatively similar to that obtained when adding the 
contribution of the pair current in non-relativistic calculations 



TT7|. 



Figure 35 illustrates the influence of B{Q^) of the higher order terms than g"^. Their 
importance indicates the need to carry out the calculations beyond the iterative procedure. 
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Fig. 34. The deuteron tensor polarization T20 at 6* = 70°. calculated by relativistic formu- 



las for form factors. The designations of the curves are the same as in figs. p2| and B3 



128 




0.0 2.0 4.0 6.0 8.0 10.0 



Q**2 (GeV/c)**2 

Fig. 35. The structure function B{Q'^). The sohd curve is the same as in fig. ^ The 
dashed curve incorporates the quadratic terms in f^. 



129 



The function ^(Q^) is less sensitive to this contribution. 

The results presented in this section should be considered only as indicative. They 
can be improved in many respects. An exact solution of the equation for the deuteron 
and scattering state wave functions is indeed feasible. It implies the redefinition of the 
parameters of the NN kernel. However, we believe a realistic estimate of relativistic cor- 
rections to the deuteron can be performed along the lines presented in this section. It 
necessitates a careful choice of the NN potential to generate the kernel (energy-dependent 
OBE potential for instance), and a good starting point for the non-relativistic wave func- 
tions (with kinematical relativistic corrections indicated in sect. [4.2| ). This may imply 
also the use of a PV nNN coupling to be consistent with chiral symmetry constraints in 
leading order. This necessitates also a careful treatment of second order contributions, 
as discussed in sect. |4.3.4 Finally, the calculation of electromagnetic observables should 
include retardation effects as indicated in fig.|22. 
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Chapter 8 

Concluding remarks 



We have reviewed the general properties of the covariant formulation of LFD, and its 
application to relativistic few-body systems. We have first detailed the general transfor- 
mation properties of the state vector, both kinematical and dynamical, and emphasized 
the most important features brought about by the covariant formulation of the state 
vector on the plane uj-x = 0. The explicit covariance provides similar advantages over 
the standard formulation, as the Feynman graph technique in comparison to old fash- 
ioned perturbation theory. We have presented a covariant light-front graph technique to 
calculate physical amplitudes. 

The covariance enables us also to construct bound and scattering states of definite 
angular momentum, separating kinematical and dynamical aspects of this problem. The 
angular momentum of a system is a property of the wave function relative to transfor- 
mation of the coordinate system and therefore has in our formalism a purely kinematical 
nature. Dynamics is involved in the composition of this angular momentum from the 
spins of its constituents. 

We illustrate our formalism with the most simple and adequate system, the two 
nucleon system, and construct explicitly the deuteron wave function and the ^5*0 scattering 
state. In both cases, the dominant relativistic components of the wave function originating 
from its dependence on the light front position (the variable n) has been found. We also 
extend our formalism to the tt and nucleon wave function in terms of valence quarks and 
antiquarks. 

This formulation is essential in order to extract the physical amplitude for all elec- 
tromagnetic transitions. This is mandatory in all practical calculations in order to avoid 
non-physical contributions coming from the definite choice of u. We show how this applies 
to the electromagnetic form factors of spin 0,1/2 and 1 systems, as well as for the deuteron 
electrodisintegration. We show also in this last example how one can recover easily, in a 
1/m expansion, the well known results obtained in the non-relativistic phenomenology in 
terms of meson-exchange currents (the so-called pair term). The dominant n-dependent 
components of the wave functions are essential in order to reproduce these results. This 
shows that the calculations which omit this dependence cannot be considered as realistic. 

These developments have been made particularly simple, and intuitive, by the three- 
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dimensional nature of the formalism, and the absence of vacuum fluctuations. This enables 
a transparent link with the non-relativistic framework developed over the last twenty years 
in the microscopic structure of few-nucleon systems. 

We expect many new developments in the near future along the lines presented in 
this review, both in nuclear and particle physics. One of them concerns the particular 
way to implement the independence of various results on the light-front orientation. Some 
restricted examples (determination of the A'^A'^-inter action, normalization of the wave 
function, . . .) have been considered in this review to show how this canbe done, but 
more extensive work is required. Another related aspect, which is related, concerns the 
construction of the generators satisfying, in the truncated Fock space, the commutation 
relation of the Poincare group algebra, that also provides independence of observables on 
the light-front orientation. 

This feature of LFD - the dependence of approximate on-shell amphtudes on u, 
which has to be excluded from observables, - has, however, a positive aspect. It gives a 
quantitative measure of the incompleteness of a given approximation. The model calcula- 
tions in other approaches where this feature is absent are often still incomplete, but this 
does not manifest itself explicitly. 

In nuclear physics, the exact calculation of the two-body wave function is worth to 
be done. This however necessitates to reparametrize the NN potential starting from the 
same relativistic formalism. This is essential in order to have quantitative predictions for 
electromagnetic observables in the few GeV range. In particle physics, the structure of 
the pion and nucleon can be investigated easily in this framework, and electromagnetic 
form factors calculated in the high momentum region. Application to virtual compton 
scattering can also be considered, as well as the extension to axial transitions. One 
other particularly interesting application is the structure of heavy quarkonium, like the 
J/ijj, in order to investigate the importance of relativistic dynamical corrections, i.e. the 
importance of higher Fock states in the J /ip wave function. 

Last but not least, the application of the covariant formulation of LFD to field 
theory should be investigated. 
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Appendix A 
Notations 



Wc use the standard covariant /contravariant notation, with imphcit summation over re- 
peated indices, according to: 



/ 1 













-1 














-1 





V 








-1 / 



The Dirac matrices satisfy the anticommutation relation: 

The standard representation for the Dirac matrices and spinors is used: 



(A.l) 



7o 



1 
-1 



7 



a 
-a 



75 = -«7o7i7273 



-1 
-1 



(A.2) 



(A.3) 



where B are the usual Pauli matrices. Note the minus sign for the definition of 75. 
The tensor is: 

Throughout the review the following convention is taken: 



(A.4) 



(A.5) 



The Dirac spinors u"{k) and v''{k) for spin 1/2 fermion and antifermion respectively 
satisfy the following Dirac equations: 

{k-m)u''{k) = %f{k){k-m) = {k + m)v''{k) = v''{k){k + m) = 0, (A.6) 

with the normalization: 

u''{k)u''\k) = 2m6aa' = -v'' {k)v''\k) . (A.7) 
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They are given by: 



/ 



1 

— * 

a-k 



\ (Ek + m) 



( 



\ 



a-k 



(A.8) 



where w'^ are the two-component spinors and e^ = \Jk'^-\- m?. The summation over spin 
indices writes in that case: 



<{k)u%k) = {k + m)ap, E <ikH{k) = {k- m)a/3. 



(A.9) 



The charge conjugation matrix is defined as: Uc = 7270- 

We use the following properties: ayO*^ — —aay, Uc^j^ = —JnUc, where the superscript 
t means transposition. 

The polarization vector e^^^ (p) of a spin 1 particle of mass m satisfies the following 
general relations: 



m 



2 ■ 



(A.IO) 
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Appendix B 

Relation to other techniques 



B.l Relation to the Weinberg rules 

As we already mentioned, the approach considered in this report is a generahzation of 
the usual LFD obtained by the replacement of the particular value of the light-front 
vector oj = (1,0,0,-1) by the most general one u = (cuq, cJ) with the only constraint 
tu^ = 0. Hence, the graph technique developed in section ^]2| is the direct generalization 
of the Weinberg rules . These rules have been obtained from old fashioned perturbation 
theory in the infinite momentum frame. We shall show below for the case of a spinless 
particle how to transform the amplitudes calculated in the covariant light-front graph 
technique to the Weinberg amplitudes. This explicit transformation was found in [Q. 

Any diagram of the covariant light-front graph technique can be represented as a 
sequence of intermediate states, each of them containing, besides the particles, one spurion 
line and only one. An example of such contribution is indicated in fig. 




Fig. 36. Amplitude with one intermediate state shown explicitly. Other intermediate 
states, if any, are absorbed by the amplitudes Mi and M2. 

The expression for the corresponding amplitude has the form: 
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dr " 

X— n e{^-k,) 5{k] - m]) d% . (B.l) 

ZTT^r — le) 

If a particle passes through a vertex without interaction, (hke in fig. |l^ the particle in 
the lower part of the diagram propagates without interaction in the upper vertex), the 
intermediate states can still be factorized. This is ensured by the following representation 
of the propagator: 



(• • ■) e{uj-k) 6{k' - m') d^k = y (• ■ ■) {2uj-ki) e{uj-ki) 6{ki - m') d% 

x6^^\ki + UJT- fca) dr e{u-k2) 6{kl - m^) d^'k-i , (B.2) 

which reduces the amplitude with the non-interacting particle to the form of eq.( p.l| ). 
From the conservation equation I] Pi + = in |B.1| , it follows that: 

Integrating ov6r (ih()j we get: 



j 



In order to transform the expression ( [B.3| ) to a form analogous to old fashioned pertur- 
bation theory in the infinite momentum frame, we introduce the variables: 

Ri = ki - Xi'^Pj, Xi = u-ki/{uj-^pj) , (B.4) 

satisfying Y^Vi = = ^ and Ri-u = R^^^-uj = with R = {Rq, Rj_, i^y), where Rj_-uj = 
and R\\ is parallel to cu. Taking into account that R^ = —R\, we can write the denominator 
of eq. (|B.3|) into the form: 



(B.5) 



This form coincides with the one given by the Weinberg rules 0. The factor d^kj/2ej in 
these variables turns into d'^Rjdxj/2xj, the domain of integration over dxj being from 
to 1. The expression 

{2uj-J2p,)6^^\J2p^ - E + 

turns into 

By this way, one can transform the expression of any intermediate state. The equivalence 
between the expression for any amplitude in the covariant formulation of light front dy- 
namics and the one given by the Weinberg rules is thus realized. The vectors R± play here 



136 



the role of the momenta transverse to the infinite momentum, and the variables Xj, are 
analogous to the fractions of the particle momenta with respect to the infinite momentum. 

Note that the factorization formula ([B.2D is generalized for the propagator of a spin 
1/2 particle as follows: 

(. ..){k + m)5ik^ - m^)d^k = [{.. .)^e(iu-ki)5(kl - m^)(ki + m)d^ki 

J m 

xS^^\ki + UJT- k2)dTe{uj-k2)5{kl - m'^){k2 - ujt + m)d^k2 . (B.6) 



The part —Cot in ( [B.6|) corresponds to the instantaneous interaction which appears in 
light front dynamics, since the fermion line with the momentum k2 now "extends over a 
single time interval" (see section |2.2|). 



B.2 Relation between the Feynman amplitudes and 
the Weinberg rules 

It is also instructive to derive the Weinberg rules directly from the Feynman diagrams. 
This allows one to avoid some subtle problems of light-front quantization. The direct 
relation between Feynman and the light-front amplitudes was shown for a few cases in 



ref. ||118|| and investigated in detail in ref. ||119|]. The problem of renormalization was 



investigated in ref. |120|. We illustrate this relation by the simplest case of the triangle 



Feynman graph shown in fig. This example has been considered also in ref. | |121| . Up 



to a multiplicative factor, the corresponding Feynman amplitude for scalar particles of 
mass m has the form: 

F^=f^ I . (B.7) 

J (27r)^ [k"^ — m? + ie){{k — — m? + ie){{k — p'Y — m? + ie) 

Introducing the new light front variables: 

K = ^ , fc- = ^ , ki_ = (fci, fca) , (B.8) 

we rewrite eq.( [B.71 ) in the form: 

dk4.d'^k I , , , r dk_ 1 



(2vr)3 '^'^ = I ^,{k.-H,){k.-H2){k.-H,y ^^-'^ 



where 03 = 'Sk^ik^ — p+){k^ — p'j^) and 

k±'^ + rrP' — ie 



Hi 



H2 = P- + 



2k+ 

{kj_ — + — ie 



n - W I {k± - p'^y + - ie 
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Let = p''^ = rn? and = {p — p'Y < 0. For simplicity we put g+ = 0, i.e., p+ = p\. At 
/c+ < and A;+ > p+ all the poles in k_ of the integrand in eg. ( P.9|) are above or below 
the real axis. Hence, the integral differs from zero < A;+ < p+ only. Application of the 
residue theorem gives: 

Ma = rJL . fB.ll^ 



{H^-H2){Hi-H; 



Introducing x = k^/p+, we get: 
FAit) - 



(27r) 



d'^k I dx 



k'j_ + 
x(l — x) 



m 



{kj_ — xKy + 
x(l — x) 



m 



2x(l — x)' 



[B.12) 



where t = = — A^. The form ( [B.12| ) reproduces the amplitude given by the Weinberg 
rules. 




Fig. 37. Box Feynman graph. 



The next step is the box diagram of fig. ^ This example is considered in ref . ||119|| . 
The amplitude has the form: 



D{ki, m)D{k2, n)D{k3, m)D{k4, n) 



d^k 



:b.i3) 



with -D(fc, m) = i/{k'^ — + ie) and ki = k, k2 = k — I, k^ = p + q — k, k^ = k — p. 
Particles are supposed to have mass m or fi depending on their appearance in the s or t 
channel. This amplitude can be represented analogously to ( |B.9| ): 



Fu 



where 04 = 16/c+(A;+ — l+){k+ — p+ — q+){k+ — p+) and -ffi_4 are analogous to (|B.10|) : 

^ + — ie 



Mr 



dk. 



1 



B.14) 



Hi 



2k. 



138 



Hi = p- + q- + 



2(fc+ - /+) 
{k± — Pj_ — q±)^ + rn? — ie 

'^{K -p+- (1+) 



ft = p_+(^^l!±4^. (B.15) 

^K-P+) 

Let p+ > Z+. Then at /c+ < and at /c+ > p+ + g+ the imaginary parts ImHi_4 have 
the same signs and, hence, the integral for Mn is zero. At < A;+ < Z+ < p+ < p+ + g+ 
there is one pole in the variable /c_ at the down half plane relative to the real axis and 
three poles are at upper half plane. The residue reproduces the expression given by the 
Weinberg rules for the diagram analogous to fig. p3, as in the case of the triangle graph 



indicated in fig. 21 




At A;+ > > /+, but k+ < p+ + q+, the amplitude corresponds to a diagram similar 
to fig. ^ but up side down. Note that if the poles have imaginary parts of the same signs, 
the corresponding lines of the time-ordered graph have the same (clockwise or counter 
clockwise) directions. For the poles with opposite imaginary parts, the directions of lines 
are opposite. 

New feature appears at /+ < A;+ < p+, when two pairs of poles are at the opposite 
sides of the real axis. In this case, any given residue does not reproduce any time-ordered 
diagram. However, the sum of them can be transformed ||119|| to a form consisting of two 
other parts, each of them coinciding with the amplitudes given by the Weinberg rules 



for the diagrams analogous to those shown below in figs. ^ We explain this below 



in more detail in the case of the covariant formulation of LFD. The general algorithm is 
developed in ref. [|119|| . 
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B.3 Relation between the Feynman amplitudes and 
the covariant light-front graph technique 

B.3.1 Spin 

Like in the case of the Weinberg rules, the Feynman approach and the covariant hght-front 
graph technique give different ways of calculating one and the same S'-matrix and, hence, 
give one and the same on-shell amplitude. However, these amplitudes are represented in 
different forms. It is therefore instructive to transform explicitly the Feynman amplitudes 
to the form of covariant light-front amplitudes. 

We will explain this transformation in the simplest case of spinless particles for the 
box diagram indicated in fig. |3^. The analytical expression for this diagram is given 
by ( [B.13| ). For this transformation, we first represent the Feynman propagator in the 



following form ||122 

iD{k,m) = 



^2 ^2 _|_ 



T — It 



+ 



T — le 

Substituting ( [Bl6D in ( [BlsD , we get: 



{uj-p) 5{p^ — m^) d'^p dr . 



:b.i6) 



d^k 



5(4 



+ 



5(4 



5(4 



+ 



5(4 



5(4 



+ 



X 



+ 



5(4 



5(4 



5(4 



[ujTi + k-pi) 



Ti — le 
{ujTi - k-pi] 



Ti — It 
{UJT2 + k-l - P2) 



T2 - le 
{1JJT2 - k + l -P2) 



0{ijj-pi) 5{p\ — m^) 
6{ijj-pi) 5{p\ — m^) 
9{uj-p2) S{p 



--(la) 



--{lb} 



d'^pidri 



--(2a) 



T2 - le 
{UJT3, + k-p- q-p-i) 



e{uj-p2) 5{pI - /i^ 



Ts - le 
(cjTg -k + p + q-p^) 



T3 - It 

{uOTi + k ~ p - P4) 



Oi^j-ps) 5{pI - rr?) 



d'^P2dT2 



--{3a) 



--{3b) 



d'^PsdTs 



T4 — le 
{ujTi ~ k + p- P4) 



9{uj-pi) 5{pI - 1/ 



9{u-pi) 5{p\ - /i^ 



= (4a) 
= (46) 



^d^PidTi . (B.17) 



r4 — le 

Let us first consider the product of the first items in the four braces in eq.( p.l7D . Using 
schematic notations, we have: 



(la) (2a) (3a) (4a)... 



d^k 5^^\uJTi + k-pi) S^^\UJT2 + k-l-p2) 

{2tiY Ti — T2 — ie 
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5^^^ {cur^ + k — p — q — p-i) 6^^^ {iOT^ + k — p — p^) 



— le 



■)dTidT2dT^dTi . (B.18) 



Tg - le 

The functions 6{uj-pi), 5{pf — yU^) and all other factors in ( p.l8|) are absorbed in 
exclude ri from the numerator, we make the following change of variables 



To 



k + UJTi = k' , T2 — Ti = t!- 



21 



r-i - n 



'3' 



ta - n 



r, 



4' 



after which we find: 



(la) (2a) (3a) (4a) 



■)dTi 



(ri - ie){T2 + Ti - ze)(r3 + ri - ie){T'^ + ri - ie) 



(B.19) 

The numerator (■ ■ ■) in ( |lj.l9| ) does not depend on ti. All the poles of the integrand in 
( p.l9| ) relative to ti are above the real axis. Therefore the integral ( |B.19| ) equals to zero. 
Similarly we get /(16) (26) (36) (46) = 0. 

We suppose that all the external momenta are on the mass shells and, for con- 
venience, uj-p > uj-l. The expression (|B.17| ) contains 16 items. However, only three of 
them differ from zero, namely, /(la) (26) (36) (46), /(la) (2a) (36) (4a) and /(la) (2a) (36) (46). 
In all other items, except for the ones considered above /(la) (2a) (3a) (4a) and 
/(16) (26) (36) (46), we get the products of the 6'-functions with mutually incompatible 
restrictions on the scalar products of uj with four-momenta. Let assign to the items (a) 
the arrow counter clockwise in the loop and to the items (6) - the clockwise arrow. Put 
on the external lines with the momenta p and q the arrows incoming to the box and to 
the lines with momenta I and p + q ~ I - the outcoming arrows. Then the items which 
disappears due to the ^-functions correspond to the diagrams containing at least one ver- 
tex with all incoming or outcoming lines, which are interpreted as vacuum fluctuations. 
One of these graphs, corresponding to the product /(16) (2a) (36) (46) is shown in fig. R^. 



(3b) 



(4b) > ^ 



(2a) 



1 



(1b) 



Fig. 39. Time-ordered box graph with the vacuum vertex No. 2. 



The vertex 2 in fig. ^ corresponds to the creation of particles from the vacuum. The 
items /(la) (2a) (3a) (4a) and /(16) (26) (36) (46), which are zero, correspond to the loops 
with all the lines directed clockwise or counter clockwise. Now consider the item 



(la) (26) (36) (46)... 



5(4) {ton + k-pi) 5(4) {iUT2 + 1 - k ~ P2) 



T2 



It 
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X : : [■■■)dTi. 

T3 — le r4 — le 

(B.20) 

Again after replacing k + ojti = k', r2,3,4 + n = r2 3 4 we get: 

;ia)(26)(36)(46) . . . = / — ^ , (B.21) 

where (■ ■ ■) does not depend on ri. Calculating the residue at ri = ie and deleting the 
prime over r's we find: 

J (la) (26) (36) (46) ... = ij e{uj-k)6{k^ - m^) 

xe{uj-{l - k)) 6{{l -k + uT^f - y?) 
x6{uj-{p + q- k)) 5{{p + q-k + un)'^ - vr?) 
xe{ujip - k)) 5{{p -k + unf - fi^) 
dT2 dr^ dr^ d^k 
(ra - ie) (rg - ie) {T4 - ie) (27r)3 ' 

(B.22) 

The expression (|B.22|) , after separating the factor i (due to its presence in the defi- 
nition ( p.54|) of the amplitude), exactly corresponds to the light-front amplitude indicated 



by the diagram in fig. ^ The integral /(la) (2a) (36) (4a) is calculated analogously (by 
integrating over the variable T3 instead of T4) and corresponds to the diagram having the 
form of the overturned trapezium (relative to fig. pS]) . 




(0X4 )lf 



Fig. 40. Time-ordered box graph. 



There are two other diagrams in the light-front dynamics shown in figs. ^ and ^ 
We can find both of them in the item /(la)(2a)(36)(46). After similar transformations 
we get: 

/ (la) (2a) (36) (46) . . . = / — r^^-^ — , (B.23) 

J J {n-ie){T2 + Ti-te){Ts-Ti-ie)in-n-ie) ' 
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Fig. 41. Time-ordered box graph. 



Taking the sum of the two residues at ri 



la) (2a) (36) (46). 



1 



+ie and ri = — r2 + ie, we obtain two items: 
1 



(B.24) 



T2T3T4 T2{t2 + Ts){t2 + T4) 

Each of them does not correspond to any diagram. However, ( [B.23| ) indeed gives the sum 



of the diagrams in figs. HQ and 41 



To transform it to the appropriate form, we shall use the simple formula: 

1 



1 

a/3 



1 ^ 
( — \ — 

a + (3^a (3' 



:B.25) 



This formula should be applied recursively to pairs of factors in the integrand having the 
poles in the integration variables with the imaginary parts of opposite signs. For example. 



applying (|B.25|) to l/(/?2/?3/54), first to 1/(32(^3 and then again to l//?2/54, we get: 
11 1 1 



+ 



+ 



(B.26) 



We apply ( p.26|) to the integrand of (p.23| ), i.e., substitute {32 = T2+Ti—ie, = r^—Ti—ie 
and /54 = r4 — Ti — ie. The factors 1//32 and /Ss in 1/^2^3 have the poles relative to the 
integration variable ri with the imaginary parts of opposite signs. The same is true for 



the factors P2 and in 1//32/34- By this way, we transform the integrand of (|B.23| ) to the 
form which after integration over dri and standard replacement of variables leads to the 
expression: 



la) (2a) (36) (46)... 



m 



e{uj-k)6{k'^ - m^) 

x9{uj-{k - l))6{{k - l + ujT2~ tOTsf - 
xe{iu-{p + q- k))6{{p + q-k + uJT^f - 
xO{uj-{j) - k))5{{p -k + uTif - jj,'^) 
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dT')dT-idTA 



X 



+ 



{t2 - ie)(r3 - ie)(r4 - ie) {2tc)^ 

e{uj-k)6{k'^ - m^) 

xOiioik - l))6{{k - / + tuTg - ujn)y - jJ) 
x9(uj-(p + q — k))5{{p + q — k + ujt2 + lut^ — lutsY — 
x9{Ljip - k))5{{p -k + unf - /i^) 
dT2dT3dT4 d'^k 



X ■ 



B.27) 



{t2 - ie)(r3 - ie)(r4 - ie) (27r)3 ' 
Two items in eq. (|B.27|) exactly correspond to the amplitudes for the diagrams of figs. ^ 



and The method can be generalized to more complicated diagrams. The general case 
is considered in ref. ||119|| and the renormalization in ref. [|120| . 



B.3.2 Spins 1/2 and 1 



For spin 1/2, one should use the following representation of the propagator pS 

k + m 



iS{k) 



J^2 _ ^2 _|_ 

6^^\u;T + k-p) 



+ 



T 


- ie 






- k 


-p) 


T 


- ie 






+ k 


-p) 


T 


- ie 




5(4) (cur 


-k 


-p) 


T - 


- ie 





(p — Cjt + m)9{uj-p) 5{p^ — m^) d'^p dr 
ip — Cjt — m)6{uj-p) dijP' — m^) d'^p dr 
{p + m)6{(jj-p) dijP' — m^) d'^p dr 
(m — p)6{ijj-p) dijp^ — m?) d^p dr 



2{u-k) ■ 



^B.28) 



In the representation given by eq.( |B.28D , the propagator contains the spin parts {p + m) 
and (m — p) with the on-shell four-momentum = and the extra contact term 
—uj/2{LJ-k), like in the rules of the graph technique, sect. |2.2.^ . At u = (1,0,0, —1) this 
contact term coincides with the one given in refs. [^, [123| , |124|| . 
For the propagator with spin 1, we have the representation: 



iD^y{k) 



k^ky 
m 



_ ^2 _|_ 

5^^\uT + k-p) 



T — le 



-9t^u + 



{p - ujt)^{p - u;t)u 



6{uj-p) 5{p^ — m^) d'^p dr 
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+ / : -g^iu H 5 6'(t^-j9) 5(p ~m ) d pdT 



T 


- ie 






+ k 


-P) 


r 


— ie 




5(4) (cjr 


- k 


-P) 



(B.29) 



-9^lu + 



^^^^ I 6iuj-p) 6(p^ — m?) d^p dr 
J 



+ / ^ -9i,u + yi^^-p) d{p -m) d p dr 

2[uj-k)m^ 4:[ijj-kYm'^ 
The item in the last hne of ( |B.3CI| ) is the vector contact term. 



Let now the particle 1 in fig. 37 has spin 1/2. Then the contact term for the particle 
1 does not contribute to the amplitude indicated in fig. Indeed, this amplitude was 
calculated in eq.( p.21|) , as the residue at ti = ie in the pole of the propagator of particle 



1. Replacing this propagator by the contact term, we see that the only pole above the 
real axis disappears and the amplitude turns to zero. However, the contact terms for the 
particles 2,3 and 4 contribute to the diagram of fig. |3^. According to (|B.28| ) they can be 
incorporated by replacing the spin part of the propagators {p ± m) by {p — ut ± m) . 

These facts agree with the general rule given in section |2.2| : the contact terms have 
to be associated with these particles which lines are directly connected by a spurion 
line. This also coincides with the rule given in ref. |3^: contact terms modify only the 
propagators corresponding to the "lines extending over a single time interval" . The lines 



of the particles 2,3,4 in fig. ^ extend over single time intervals and are associated with 
contact terms, while the line for the particle 1 extends over three time intervals and 
is not associated with any contact term. Another example is given by the diagram of 
fig. In this case, the contact term contributes for the exchanged particles (if they 
have spin) and does not contribute for the particles shown by the horizontal lines. On 



the contrary, in fig. |41| the contact term contributes for the horizontal internal lines and 



does not contribute for the exchanged particles. 
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Appendix C 

Relation between the deuteron 
components 



We give here for completeness the relations between the six invariant functions coming in 



the decomposition of the wave function (|5.1|) , (|5.2| ) and (|5.6| ). These two representations 



are connected by eq. (|5.3| ) and coincide with each other in the system of reference where 
V = 0. In this system, we have ki = —k2 = k. 

The deuteron polarization vector e^^\p) in the representation with X = j = x,y, z 



in the rest system has the form: e'f\o) = 6j, eo"'''(0) = 0. In an arbitrary system where 
p 7^ 0, it can be obtained by the Lorentz transformation (cf. ( ^.11[ )): 



Ki) 



ip) 



Ki) 



(0) + 



(e(j)(0)-p) 



P 



M{M + po) 



^^\P) 



Ki) 



(0)-p 



M 



(C.l) 



In the system of reference where "P = 0, we have: 



Po 



4£i 



P 



-n- 



With these expressions for p we get from eq. ( |C.1| ) : 



Ki) 



(p) =si + 



{2ek - M f 
AMsk 



TljTlj , 



e^Ap) 
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AMsk 



-rii 



(C.2) 



The polarization vector ( |C2| ) satisfies the general relations ( |A.10 ). 

Substituting in eqs. (|5 .11 , 15.21 ) the expressions for the nucleon spinors and Dirac matri- 
ces given in appendix ^ and the deuteron polarization vector given above, we obtain the 
structure of the wave function given in eq.( |5.5| ). Comparing the coefficients of identical 
spin structures, we can find the expressions of in terms of /j. In the approximation 
where M ?a 2m, they have the form: 



m^{2€k + rn) 
4ekk^ 



f2 + 



m 



(C.3) 
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^2 = ^(\/2/i-/2-/3-2^/4) , (C.4) 



V2k^ {2ek + m)k^ {el + AmSk + w?)k 
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Appendix D 

Two-body kinematical relations 



In the calculation of the electromagnetic amplitude we have to express the scalar products 
of the four vectors a;, ki, ^2, k[, k'2-, p, p', q with each other in terms of the integration 
variables and the four-momentum transfer squared. The variables and, hence, the kine- 
matical relations for the one loop diagram of fig. |18| (impulse approximation) and for the 
two loop diagram of fig. are different. 



D.l One loop diagram 



The one loop diagram is shown in fig. |T8[ All four-momenta are on the corresponding 
mass shells: 

kl = m\, kl = k'^ = ml, / = p'^ = M\ = -A^ = -Q^ = t, 00^ = 0. 

We consider the general case of different masses, M being the mass of the initial state, 
while Ai is the mass of the final state. 

We start from the four- vector Ri already defined in eq.( p.l6D : 

Ri = kl — xp, X = uj-ki/uj-p. 

Since Ryuj = 0, it can be represented as Ri = (i?o, -Ry) with R±-uj = and 
-R|l is parallel to uj. Similarly we represent q = (go,A,gj|) with A-u = 0. With these 
definitions, we immediately get (with u-q = 0): 

u-ki = xu-p, u!-k2 = (1 — x)u!-p, uj-k'2 = (1 — x)uj-p, uj-p = u-p' . (D.l) 

We get for the various scalar products the following relations. 
Scalar products with ki. 

krk2 = {s — m\ — m\) /2, 

kyk'2 = {s' — m\ — rn\) /2, 

krp = xM^/2 + {l-x)s/2 + {ml - ml) /2, 

krp' = krp + ki-q, 

kyq = -R_L-A + xp-q , (D.2) 
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with s and s' given by: 



X 1 — X 

where i?Y = R± — xA. 

Scalar products with 

k2-k'^ = ml + k2-q + {l-x){s'-M^-s + M^)/2, 

k2-p = (1 - x)M^/2 + XS/2 + {ml - ml) /2, 

krp' = krp + k^q, 

krq = Rl-A + {1 - x)p-q. (D.4) 
Scalar products with k2- 

K-p = K-p' - K-Q^ 

k'^-p' = {l-x)M^/2 + xs'/2 + {ml~ml)/2, 
k'^-q = /c2-g-Al (D.5) 



Scalar products with p: 



Scalar product with p': 



p-p' = {A^ + M^ + M^)/2, 

p-q = {K'^ - + M^)/2. (D.6) 
p'.q = {M^ - - A2)/2. (D.7) 



D.2 Two-loop diagram 

The two-loop diagrams indicated on fig. corresponds to the contribution of the contact 
term. One of them is reproduced in fig.^ where all the momenta are explicitly indicated. 

Some of the relations coincide with the ones given above. We repeat them here for 
completeness. For the scalar products of uj we immediately get: 

u-ki = xuj-p, u!-k2 = (1 — x)uj-p, u-k[ = x'u-p, uj-k'2 = (1 — x')uj-p, u-p = u-p, oj-q = 0. 

^ (D.8) 

All the following invariants can be transformed to the variables k, k' and n by the 
substitution: 

— * 

— * 7 

= k - (fi.k)n, X = ^^{l —- ); (D.9) 

El + 62 varepsilorii 

R'± = k'-{n-k')n,x' = -^{l~^), (D.IO) 
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p ?^ k2 P 

Fig. 42. Contact term of fig. p9|(a) contributing to the deuteron electromagnetic form 
factors. 



taking into account that A-n = 0. We denote in this equation Si = \Jk? + mf, and 
similarly for quantities with prime. We thus get the following scalar products. 
Scalar products with ki. 

^i'^2 = {s — m\ — m\) 

kyk[ = -R^-R'^ + xAr'^ - x'A-Rs_ + x(l - x')s'/2 + x'{l- x)s/2 + xx'A V2, 

krk'2 = krp' + x{s' -M^)/2-krk[, krp = {I - x)s/2 + xM^ /2, 

ki-p' = kyq + ki-p, kyq = —R±-A + xp-q. (D.H) 

Scalar products with k2: 

k2-k[ = k[-p + x'{s-M'^)/2~krk[, 
krk'^ = krp' -krk[ + {l-x){s' -M^)/2, 
k2-p = XS/2 + (1- x)M^/2, 
krp' = k'^-q + k^p, 

krq = R±-A + {l-x)p-q. (D.12) 
Scalar products with k'^. 

k\-k'2 = {s' — m\ — m^) /2^ 
k[-p = k[-p' — k[-q, 
k[-p' = (l-x')s' /2 + x'M'^/2, 

k[-q = -R'^-A + x'p'-q . (D.13) 



Scalar products with k2. 
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k'^.p' = x's'/2 + {l-x')M^/2, 
k'2-q = R'^-A + (1 — x')p-q. 



(D.14) 



The scalar products with p and p' are the same as ( p.6|) , ( p.7|) . 

The kinematics for the final plane wave contribution (neglecting the final state en- 
ergy) is obtained from the above formulas by the substitution Ai = 2m, R'^ = 0, x' = 1/2 
(or k' = 0). The kinematics for the elastic deuteron form factors is obtained by the 
replacement A4 ^ M. 
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Appendix E 

Three-body kinematical relations 



The scalar products appearing in the calculation of the amplitude indicated in fig. ^ are 
expressed through the integration variables Rij_,Xi and A. We give them below in the 
general case of different quark masses. For p-p' we have: 

p-p' = M^ + Ay2. (E.l) 

In order to find p-ki {i = 1, 2, 3), we calculate (fcj — Xip)"^ = mf — 2xip-ki + x'^M'^ = Rj = 
—Ri±- From here we find: 

p-k = ^^^^ + \m'x, ■ (E.2) 
Other scalar products are found analogously. We have: 

p'-h = p-ki - Ri^-A + AV2 , (E.3) 

P'-K = ^^^^^ + \m'x. , (E.4) 

p-k[ = p'-k[ + R!^^-A + XiAV2 , (E.5) 

ki-k'j = —Ri-R'j + Xip-k'j + Xjp'-ki — XiXjp-p for i, j = 1, 2, 3. (E.6) 

The scalar product ki-kj is obtained from eq. (Pi).6| ) by deleting primes at all variables, 
the scalar product k[-k'- is obtained by the following replacement in eq. ( [E.6|) : Ri R'^, 
ki ^ k[, p ^ p' . The variables Xi are the same for initial and final states. 

Let us now express R[^ through A and the integration variables Ri_[_ and Xj. The 
expressions are not symmetrical relative to the quark 1, which interacts with the photon, 
and to the spectator quarks 2 and 3. We have R!2 = k'2 — X2p' , R'^ = k'^ — x^p' . With 
k'2 = ^2, ^3 = ^3 we find R'2 = R2 — X2q, R'^ = R3 — x^q. With R[ = k[ — xip' = 
i?i + (1 - Xi)q + uj{t' - r), we get: 

R[^ = Rl± + (1 - Xi) A , R'2,^= R2,3± - X2,3^ . (E.7) 
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Note that R'^^^ can be also found from the relation R'^_^ = — i?2± ~ -^3± • The arguments 
g| (i = 1, 2, 3) of the initial wave function are given by 

-»2 2 2 

where e^. is the energy of the particle in the reference system where V = ki + k2 + = 0. 
It can be represented in the invariant form as £g. = V-ki/ Ai with V — p -\- cut, t — 
{M"^ - M^)/{2u;-p), that gives: 

= \p-h + - M^)\IM , (E.8) 

where is the invariant mass squared: 

= (fci + ^2 + k^f = i: ^^±^ . (E.9) 

i=l 

The arguments gf of the final wave function arc constructed similarly from the final 
momenta. 
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